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p.  12  Equation  (2-1)  should  be 

l  T 

f(x)  ■  lim  —  /  f(jc,  t)  dt 
T+-  T  o 

p,  16  1 st  line,  t  should  be  1/t 

p.  127  8th  line  “In  is  clear  . .."  should  be  "It  is  clear.*." 

p.  258  Lin,  J.  and  M.  Gad-et-Hale, . . . ,  J.  Phys.  Ocean  should  be 

...  J.  Geo .  Res . 
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ABSTRACT 


Because  the  global  weather  and  ocean  circulations  are  driven  by 
transfers  across  the  atmosphere-ocean  interface,  studies  of  it  remain  at 
the  forefront  of  research  efforts  in  atmospheric  sciences  and  oceanogra¬ 
phy.  Many  questions  concerning  the  momentum  and  energy  transfers  in  the 
water  surface  layer  beneath  the  air-water  Interface  remain  unanswered. 
The  present  experimental  program  was  undertaken  with  the  goal  of  learn¬ 
ing  how  the  energy,  which  is  transferred  from  the  wind,  is  distributed 
among  the  mean,  wave,  and  turbulent  flow  fields  in  the  water.  The 
interactions  between  these  three  flow  fields  were  also  examined. 

The  experiments  were  performed  in  the  Stanford  Ulnd,  Water-Wave 
Research  Facility.  A  two-component,  frequency-shifted,  laser  Doppler 
anemometer  system  was  used  with  frequency  trackers  to  measure  the  water 
velocity.  Three  major  sets  of  experiments  were  done:  1)  wind-generated 
waves  in  Eulerian  coordinates  with  wind  speeds  ranging  from  1.5  m/s  to 
13.1  m/s,  11)  wind-ruffled  mechanically  generated  waves  (22  mm  in  ampli¬ 
tude  and  at  a  frequency  of  1  Hz)  in  Eulerian  coordinates  with  wind 
speeds  ranging  from  0.0  to  6.2  m/s,  and  ill)  wind-ruffled  mechanically 
generated  waves  in  wave-following  Eulerian  coordinates  with  wind  speeds 
ranging  from  0.0  to  6.1  m/s.  A  novel  optical  wave-follower,  which 
allowed  velocity  measurements  very  close  to  the  mobile  surface,  was  con¬ 
structed  and  deployed  in  the  wave-following  measurements.  The  mean, 
wave,  and  turbulent  fields  were  isolated,  and  from  them,  other  quanti¬ 
ties  of  interest  were  computed. 

The  results  show  that  the  waves  affect  the  mean  flow,  but  the  mean 
velocity  profiles  show  logarithmic  behavior.  The  wave  field  generally 
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agrees  with  the  prediction  of  a  linear  water-wave  theory.  In  the  wind- 
wave  experiments,  the  turbulent  quantities  behave  similar  to  those  in 
flows  over  flat  plates.  In  the  mechanical-wave  experiments,  the  turbu¬ 
lent  quantities  have  different  trends  from  those  of  the  wind-wave  cases, 
which  suggests  possible  weak  wave-related  turbulence.  The  mean  wave- 
induced  shear  stress  is  negative,  which  implies  that  the  waves  augment 
the  mean  flow.  Thus,  the  waves  can  also  transfer  energy  indirectly  to 
turbulence  via  the  mean  flow.  The  turbulent  field  draws  its  energy  from 
the  mean  flow  through  the  mean  turbulent  shear  stress  and  from  the  wave 
field  via  the  wave-induced  turbulent  stresses.  Experimental  observa¬ 
tions  agree  with  simple  analysis  of  the  momentum  and  energy  equations. 
This  study  contributes  to  the  basic  understanding  of  the  fundamental 
principles  which  govern  transfers  at  an  air-water  Interface. 
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ug  Eulerian  surface  drift  velocity. 


xiv 


Lagrangian  surface  drift  velocity. 
Stokes  surface  drift  velocity 


USL 

uStokes 

uB  mean  air  free-stream  velocity. 

u+  (us-u)/u*;  mean  horizontal  velocity  defect  in  law-of-the-wall 

coordinates. 

v  instantaneous  vertical  velocity  component. 

Vin  input  voltage  to  optical  wave-follower. 

Vpot  position  output  voltage  of  the  optical  wave-follower. 

w  instantaneous  spanwise  velocity  component. 

x  streamwise  distance  or  fetch  measured  with  respect  to  the  air 

inlet. 

_x  position  vector  in  three-dimensional  space;  (x^^.xg)  or 

(x,y,z). 

x^  same  as  _x* 

x*  streamwise  distance  or  fetch  in  wave-following  coordinates  mea¬ 

sured  with  respect  to  the  air  inlet. 

y  vertical  distance  (positive  upwards)  measured  with  respect  to 

the  still  water  level. 

yn  wave  decay  depth. 

y+  yu*/v;  vertical  distance  in  law-of-the-wall  coordinates. 

y*  vertical  distance  in  wave-following  coordinates. 

(y*)+  y*u*/v;  vertical  distance  in  wave-following  law-of-the-wall 

coordinates. 

z  spanwise  distance  measured  with  respect  to  the  centerline  of 

the  channel. 

z*  spanwise  distance  in  wave-following  coordinates  measured  with 

respect  to  the  centerline  of  the  channel. 


Greek  Symbols 


au/ug;  ratio  of  velocity  amplitude  of  the  wave  motion  to  the 
surface  drift  velocity. 
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turbulent  boundary-layer  thickness. 


e 

n 

*s 

e 


s 


K 

V 


P 

T 


U! 


denotes  e  change  when  used  as  a  prefix. 

multiplicative  constant  used  in  or der-of -magnitude  analysis, 
surface  displacement  with  respect  to  the  still  water  level, 
amplitude  of  the  1  Hz  component  of  n* 

phase  lag  of  g  with  respect  to  the  wave  height  signal  n* 

von  Karman  constant. 

kinematic  viscosity. 

density. 

wave  period. 

2irf;  circular  frequency. 


Superscripts 

'  turbulent  component. 

wave  or  periodic  component. 
—  mean  component . 

amplitude. 


Special  Symbols  and  Usage 

-  when  used  after  a  variable  denotes  its  typical  value  in  an 

order-of-magnitude  estimate. 

<  > 


phase  average. 


CHAPTER  1 


INTRODUCTION 

1.1  BACKGROUND  AND  MOTIVATION 

Gas-liquid  boundary-layer  flows  occur  in  many  situations  of  practi¬ 
cal  importance.  Many  industrial  and  geophysical,  as  well  as  meteorolo¬ 
gical,  processes  involve  such  flows;  in  particular,  the  air-water 
interaction  problem  seems  to  be  the  most  common.  The  transfers  of 
momentum,  mass,  and  energy  across  an  air-water  interface  have  been  the 
subject  of  study  for  decades.  However,  because  turbulence  plays  a  major 
role  in  both  the  air  and  water  surface  layers  about  the  interface, 
progress  in  developing  a  deep  understanding  of  the  processes  governing 
the  transfers  has  been  slow.  Many  questions  remain,  especially  about 
transfers  in  the  water  surface  layer  beneath  the  interface. 

For  example,  when  wind  blows  over  water  waves,  how  is  the  energy 
which  is  transferred  from  the  wind  distributed  in  the  water  near  the 
Interface,  say,  among  the  mean  drift  current,  turbulence,  and  wave 
motion?  How  does  this  distribution  change  with  depth  and  by  what 
mechanism?  How  do  the  turbulent  and  wave  quantities  behave  near  the 
Interface?  Questions  such  as  these  are  perhaps  best  answered  by  per¬ 
forming  a  series  of  well-planned  laboratory  experiments. 

The  work  described  here  is  part  of  an  ongoing  program  at  Stanford 
to  study  and  gain  a  better  understanding  of  the  air-water  Interaction 
problem.  The  focus  of  this  particular  effort  is  on  the  water  surface 
layer  beneath  the  interface.  A  brief  review  of  previous  research  (at 
Stanford  and  elsewhere)  which  is  believed  to  be  most  relevant  is  given 
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in  the  next  section. 


1.2  PREVIOUS  RELATED  STUDIES 


1.2.1  Laboratory  and  Field  Studies 

Shemdln  (1972)  obtained  mean  velocity  profiles  In  the  wind-driven 
water  surface  layer  Inside  a  wind-wave  facility.  Neutrally  buoyant 
particles  were  used  as  velocity  tracers  under  reference  wind  speeds  of 
3.1,  5.7,  and  9.6  m/s.  He  showed  that  the  drift  profiles  were  logarith¬ 
mic.  Similar  studies  were  made  by  Wu  (1975),  who  used  floats  to  measure 
surface  and  near-surface  drift  current  and  a  pitot-static  tube  to  mea¬ 
sure  the  subsurface  drift  current.  The  reference  wind  speeds  used  by  Wu 
were  from  2.8  to  13.3  m/s.  He  found  that  the  current  immediately  below 
the  water  surface  varied  linearly  with  depth  and  the  current  near,  but 
not  immediately  below,  the  water  surface  followed  a  logarithmic  velocity 
distribution.  McLeish  and  Putland  (1975)  measured  wind-driven  flow 
profiles  in  the  top  millimeter  of  water  in  both  laminar  and  turbulent 
laboratory  flow  using  microscopic  bubble  traces.  They  confirmed  the 
existence  of  a  viscous  sublayer  (linear  profile),  but  it  was  consider¬ 
ably  thinner  than  that  at  a  solid  boundary.  Deeper  in  the  water  layer, 
they  found  that  the  mean  profile  became  logarithmic. 

Dobroklonsky  and  Lesnikov  (1975)  obtained  kinematic  characteristics 
of  wind-wave  flow  by  taking  motion  pictures  of  near-neutrally  buoyant 
polystyrene  indicator  balls  at  various  depths  in  the  water.  They 
reported  logarithmic  mean  drift  velocity  profiles  and  that  the  total 
Reynolds  stress  originated  from  two  sources:  a  wave  source  and  a  turbu¬ 
lent  source.  Later,  Anisimova  et  al.  (1978)  measured  drift  velocity 
profiles  over  various  sectors  of  wind-generated  waves  using  the  same 
wind-wave  tank  and  measurement  technique. 
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Howe  et  al .  (1982)  reported  laboratory  measurements  of  velocity  and 
temperature  fields  on  both  sides  of  a  wind-driven  air-water  interface  for 
wind  speeds  of  1.6  to  13.1  m/s.  Mean  velocities  and  turbulent  intensi¬ 
ties  in  water  were  obtained  using  a  two-component  laser  Doppler  anemome¬ 
ter.  Their  mean  velocity  profiles  when  plotted  in  law-of-the-wall 
coordinates  showed  that,  both  in  the  air  and  in  the  water,  the  profiles 
were  logarithmic.  In  addition,  the  mechanism  of  momentum  transfer  was 
affected  by  surface  roughness  changes.  In  the  water  surface  layer  the 
velocity  fluctuations  due  to  the  wave-related  motion  were  of  the  same 
order  as  the  purely  turbulent  motions.  Goossens  et  al.  (1982)  used  a 
LDA  to  measure  turbulent  intensities  and  Reynolds  stresses  in  air-driven 
water  flows.  More  recently,  similar  experiments  were  performed  by  Lin 
and  Gad-el-Hak  (1984).  However,  they  (Goossens  et  al.,  Lin  and  Gad-el- 
Hak)  did  not  separate  out  the  wave  effect  as  Howe  et  al.  (1982)  did. 

Bllven  et  al.  (1984)  conducted  a  series  of  laboratory  experiments 
to  measure  the  velocity  field  below  surface  gravity  waves.  The  series 
consisted  of  cases  in  wh'ch  waves  were  generated  by  a  paddle  actuated 
with  an  analog  record  of  wind-wave  surface  elevation  but  without  wind 
stress,  and  cases  with  the  paddle-generated  waves  and  wind  stress. 
Their  technique  permits  wave-generated  turbulence  to  be  investigated 
with  minimal  turbulence  from  other  sources,  such  as  shear  currents  and 
direct  wind  stress.  Their  study  showed  the  existence  of  a  Reynolds 
stress  which  increased  as  the  wave  steepness  increased  and  decayed 
exponentially  with  depth.  Turbulent  energy  increased  as  wave  steepness 
increased  and  can  be  modeled  by  an  exponential  decay.  The  ratio  of  the 
exponential  decay  rate  for  turbulent  energy  to  wave  particle  velocity 
decay  rate  was  less  than  one,  which  indicates  that  the  turbulent  energy 
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penetration  depth  is  on  the  order  of  the  vave  length  rather  than  the 
wave  height  as  suggested  by  Donelan  (1978). 

Field  measurements  of  water  velocity  under  wind  waves  were  made  by 
Bowden  and  White  (1966),  Shonting  (1964,  1967,  1968,  1970),  Simpson 
(1969),  Yefimov  and  Khristoforov  (1969,  1971),  Taira  (1971),  Thornton 
and  Krapohl  (1974),  and  Cavaleri  et  al.  (1978).  Bowden  and  White,  and 
Simpson  found  that  the  horizontal  velocities  were  generally  greater  than 
those  predicted  by  linear  theory.  According  to  Thornton  and  Krapohl, 
the  measured  wave-induced  velocities  were  2-4%  greater  than  those  calcu¬ 
lated  by  using  linear  wave  theory.  Furthermore,  the  phase  relationships 
among  the  vertical  velocity,  horizontal  velocity  and  the  surface  eleva¬ 
tion  were  in  good  agreement  with  linear  wave  theory.  Cavaleri  et  al. 
observed  that  the  amplitudes  of  the  vertical  and  horizontal  velocities 
were  about  10%  less  than  expected  on  the  basis  of  linear  wave  theory. 
However,  Shonting,  Yefimov  and  Kristoforov,  and  Cavaleri  et  al.  showed 
that  the  vertical  and  horizontal  velocities  were  not  in  quadrature,  thus 
indicating  the  existence  of  a  vertically  downward  momentum  flux. 

1.2.2  Theoretical  and  Empirical  Development 

Kondo  (1976)  proposed  a  hypothesis  of  hydrodynamic  similarity  for 
both  the  air  and  water  boundary  layers  adjoining  an  air-water  Interface. 
Using  the  high  frequency  components  of  the  surface  wave  spectrum  as  a 
measure  of  roughness,  he  treated  the  Interface  as  though  it  were  a  solid 
boundary  for  both  the  air  and  the  water  surface-layer  flows.  Although 
he  did  not  assert  that  the  instantaneous  orbital  motion  of  significant 
waves  had  no  effect  on  turbulent  transport,  he  proposed  that  such  a  wave 
effect  would  make  a  secondary  contribution  to  the  current. 


Jones  and  Kenney  (1977)  examined  a  number  of  one-dimensional  velo¬ 
city  spectra  obtained  In  aquatic  mixed  layers  (all  field  wind-wave 
data).  They  scaled  the  spectra  with  the  surface  stress  and  the  distance 
from  the  free  surface  and  compared  with  turbulent  boundary-layer  mea¬ 
surements  (Bradshaw,  1967).  The  good  agreement  between  the  spectra  of 
the  aquatic  mixed  layer  and  a  turbulent  boundary  layer  suggests  that  the 
surface  wave  orbital  velocities  act  merely  as  "inactive"  motions  and  do 
not  interfere  with  the  lower-frequency  stress-carrying  eddies. 

Csanady  (1978)  extended  the  law  of  the  wall  of  turbulent  boundary- 
layer  theory  to  turbulent  flow  along  a  sharp  density  interface  which  he 
termed  the  "law  of  the  interface.”  He  further  modified  this  law  and 
made  it  applicable  to  an  air-sea  Interface.  One  important  difference 
between  the  law  of  the  wall  and  the  law  of  the  interface,  according  to 
Csanady,  was  that  the  velocity  gradient  at  an  interface  was  typically 
much  less  than  the  stress  divided  by  the  molecular  viscosity.  He 
explained  this  finding  in  terms  of  the  presence  of  viscous  wavelets  on 
the  interface  which  contributed  to  momentum  transfer  through  a  combina¬ 
tion  of  "sheltering"  and  the  generation  of  an  interfacial  viscous  bound¬ 
ary  layer.  He  postulated  a  simple  model  for  the  effective  viscosity  for 
miscible  fluids  which  he  found  to  be  related  to  the  Keulegan  parameter 
and  a  wave  constant.  A  modified  model  was  also  postulated  for  the 
effective  viscosity  at  an  air-sea  interface  due  to  the  influence  of 
capillary  waves.  However,  for  flows  with  high  friction  velocity,  these 
two  models  were  shown  to  be  equivalent. 

According  to  classical  wave  theory,  wave-induced  orbital  motions 
beneath  the  interface  cannot  support  a  shear  stress.  However,  various 
experiments,  such  as  those  of  Shonting  (1967)  and  Howe  et  al.  (1981) 


indicate  lit*  exleleiue  of  weva-itidtu  cd  shear  ulruuu.  Vcsrnti.iv  an.t 
Grishin  (1972)  examined  the  two-dimensional  nonlinear  problem  of  wave 
perturbation  In  the  surface  layer  of  the  ocean.  They  introduced  an  eddy 
viscosity  coefficient  to  take  into  account  the  generation  of  turbulence 
by  wave  motion  and  computed  Reynolds  wave  stresses  by  the  method  of 
matched  asymptotic  expansions.  They  found  that  the  wave  stresses  led  to 
a  stationary  Eulerian  velocity  and  a  mass-transport  velocity.  This 
interaction  of  the  wave  field  with  the  turbulence  field  and  the  concept 
of  the  wave  eddy  viscosity  were  also  discussed  by,  for  instance, 
Navrotskii  (1967),  Kitaigorodskii  and  Miropolskii  (1968),  and  Phillips 
(1977).  A  recent  review  on  the  subject  of  wave  turbulence  was  given  by 
Katsaros  (1976). 

Liu  et  al.  (1979)  developed  a  model  for  the  marine  atmospheric 
layer  including  the  interfacial  sublayers  on  both  sides  of  the  air-water 
interface  where  molecular  constraints  on  transport  are  important.  Flux- 
profile  relations  which  were  based  on  the  postulation  of  intermittent 
renewal  of  the  surface  fluid  were  matched  to  the  logarithmic  profiles 
and  compared  with  both  field  and  laboratory  measurements.  By  using 
these  relations,  one  can  make  numerical  determinations  of  air-water 
exchanges  of  momentum,  heat,  and  water  vapor  by  employing  the  bulk 
parameters  of  mean  wind  speed,  temperature,  and  humidity  at  a  certain 
height  in  the  atmospheric  surface  layer  and  the  water  temperature. 

Street  (1979)  presented  a  simple  theory  for  turbulent  heat  and  mass 
transfers  across  a  rough,  air-water  Interface.  The  theoretical  predic¬ 
tions  and  results  obtained  from  laboratory  data  were  in  quantitative 
agreement;  however,  the  assumptions  and  data  needed  to  allow  a  predic¬ 
tion  were  somewhat  restrictive  (cf.,  Liu  et  al.,  1979).  Klotz  and 


Street  (1981a)  presented  a  solution  for  the  coupled  laminar  boundary-* 
layer  equations  of  heat,  mass,  and  momentum  at  a  smooth  Interface;  the 
numerical  solution  applies  to  flows  with  and  without  pressure  gradient 
and  uses  the  Keller  box  method.  Klotz  and  Street  (1981b)  described  the 
numerical  simulation  of  the  turbulent  boundary  layer  at  a  rough,  air- 
water  interface.  A  key  problem  lies  in  the  parameterization  of  the 
roughness  effects  at  the  interface  and  their  influence  on  the  eddy 
coefficients  in  the  near-interface  regions  (see,  for  example,  their 
Fig.  4  where  it  is  clear  that  the  heat  transfer  has  been  successfully 
handled,  but  the  momentum  transfer  has  not). 

Donelan  (1978)  measured  velocity  spectra  beneath  wind  waves  on  Lake 

Ontario  using  a  miniature  drag  sphere.  Linear  statistical  methods  were 

used  to  separate  the  velocity  into  wave  and  turbulent  parts  (cf., 

Benilov  et  al. ,  1974).  The  most  prominent  feature  of  these  spectra  is 

the  enhancement  about  the  dominant  wave  frequency.  Sufficiently  far 

-5 /3 

from  the  peak,  the  frequency  dependence  tends  to  f  ,  with  the  high 
frequencies  showing  an  increase  in  spectral  density.  Lumley  and  Terray 
(1983)  proposed  a  new  model  to  explain  the  shape  of  these  velocity  spec¬ 
tra.  Their  model  shows  that  these  features  can  be  understood  without 
recourse  to  dynamics,  but  arise  instead  from  an  essentially  kinematic 
process  in  which  frozen,  isotropic  turbulence  is  convected  bodily  by  the 
orbital  velocity  field  associated  with  the  surface  wave  motion. 

1.2.3  Conclusions 

It  is  widely  accepted  that  a  viscous  sublayer  and  a  layer  in  which 
there  is  a  logarithmic  variation  of  the  mean  velocity  profile  exist  in 
the  water  layer  at  an  air-water  interface  (Shemdin,  1972;  McLeish  and 
Putland,  1975;  Wu,  1975;  Howe  et  al.,  1980).  However,  there  seems  to  be 
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so mo  dispute  about  the  viscous  sublayer  thickness  as  compared  to  that  at 
a  solid  boundary.  McLeish  and  Putland  concluded  from  their  experiment 
that  the  viscous  sublayer  at  a  free  surface  was  considerably  thinner 
than  that  at  a  solid  boundary.  However,  according  to  Csanady  (1978), 
the  presence  of  viscous  wavelets  on  the  interface  implies  a  thicker 
viscous  sublayer  due  to  higher  energy  dissipation,  and  an  augmented  eddy 
viscosity  is  needed  to  model  the  experimental  results. 

It  has  been  suggested  also  that  the  orbital  motion  of  waves  in  the 
water  layer  has  little  effect  on  turbulent  transport  (Kondo ,  1976;  Jones 
and  Kenney,  1977).  However,  the  experiments  of  Shonting  (1967)  and 
Yefimov  and  Kristoforov  (1969)  showed  otherwise.  Howe  et  al.  (1982) 
reported  that  near,  but  not  immediately  below,  the  interface  the  wave- 
related  motions  were  of  the  same  order  as  the  purely  turbulent  motions, 
and  deviation  of  mean  velocity  profiles  from  a  1/k  slope  (k  =  0.4)  was 
apparent  when  waves  were  present  compared  to  the  wave-absent  case.  It 
is  quite  possible  that  right  below  the  interface  there  exists  a  zone 
where  the  organized  wave  motion  dominates,  hence  becoming  important  in 
transport . 

It  is  difficult  to  control  conditions  in  field  experiments.  Labor¬ 
atory  studies  such  as  those  of  Howe  were  done  using  fixed  probes,  but 
led  to  significant  results.  Experiments  with  a  progressive,  mechani¬ 
cally  generated  wave  train  under  the  shearing  action  of  a  turbulent  wind 
may  also  prove  to  be  interesting  as  shown  by  Hsu  et  al.  (1981).  How¬ 
ever,  the  interface  cannot  be  approached  as  closely  as  one  desires  due 
to  the  larger  amplitude  of  the  mechanically  generated  waves.  Further¬ 
more,  Hsu  et  al.  (1981)  showed  that  a  wave-following  coordinate  system, 
which  allowed  measurements  in  the  air  between  the  wave  crests  when 


mechanically  generated  waves  were  present,  was  crucial  for  understanding 
the  physical  processes  of  energy  and  momentum  transfer  to  the  waves.  In 
view  of  all  the  aforementioned  questions  and  arguments,  it  is  only  natu¬ 
ral  to  conclude  that  further  experimental  investigations  are  necessary. 

1.3  OBJECTIVES  OF  THE  STUDY 

A.s  a  first  attempt  to  answer  some  of  the  questions  stated  in  the 
previous  section,  a  series  of  experiments  were  undertaken  with  the  fol¬ 
lowing  objectives: 

1.  Obtain  velocity  data  (mean  and  turbulent)  in  the  water 

layer  under  wind-generated  waves  in  Eulerian  coordinates. 

2.  Obtain  velocity  data  (mean  and  turbulent)  in  the  water 

layer  under  wind-ruffled  mechanically  generated  waves  in 
Eulerian  coordinates. 

3.  Design,  construct,  and  qualify  an  optical  wave-following 
mechanism  which  enables  the  velocity  probe  to  sample  very 
close  to  the  mobile  Interface. 

4.  Obtain  velocity  data  (mean  and  turbulent)  in  the  water 

layer  under  wind-ruffled  mechanically  generated  waves  in 
wave-following  Eulerian  coordinates. 

5.  Identify  similarities  and  differences  between  results 
obtained  in  Eulerian  coordinates  and  wave-following 
Eulerian  coordinates. 

6.  Compare  measured  quantities  with  existing  data  and  theo¬ 
retical  models. 

7.  Develop  a  qualitative  model  of  the  flow  which  explains 

its  observed  characteristics. 
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1.4  SELECTION  OF  EXPERIMENTAL  CONDITIONS 

The  range  of  variables  for  Che  experimental  program  is  summarized 
in  Table  1.  Two  specific  surface  conditions  were  considered:  1)  wind¬ 
generated  waves  only,  and  2)  wind-ruffled  mechanically  generated  waves. 
Also,  two  reference  frames  were  used  for  the  second  flow  condition. 

The  air  free-stream  velocity  ranged  from  1.5  m/s  to  13.1  m/s  for 
wind-generated  wave  runs.  For  the  wind-ruffled  mechanically  generated 
wave  condition,  the  air  velocity  ranged  from  0.0  m/s  to  6.0  m/s.  The 
mechanically  generated  waves  had  an  amplitude  of  22  mm  and  a  frequency 
of  1.00  Hz.  These  values  were  chosen  so  that  the  waves  were  of  deep 
water  type  and  the  small  amplitude  wave  theory  was  applicable.  For 
almost  all  Che  cases,  the  temperature  of  the  air  was  slightly  warmer 
than  that  of  the  water. 


Table  1 


Summary  of  Range  of 

Variables  for 

the  Experimental  Program 

Mechanically 

Amplitude  of 

Case 

u«  (m/s) 

Generated 

Wave 

1  Hz  Component 

Wave 

Following 

of  the  Wave  (mm) 

1 

1.5 

no 

no 

— 

2 

2.6 

no 

no 

- 

3 

3.2 

no 

no 

- 

4 

4.7 

no 

no 

- 

5 

6.7 

no 

no 

- 

6 

9.9 

no 

no 

- 

7 

13.1 

no 

no 

- 

8 

0.0 

yes 

no 

21.9 

9 

1.7 

yes 

no 

21.1 

10 

2.5 

yes 

no 

21.9 

11 

4.1 

yes 

no 

22.1 

12 

6.2 

yes 

no 

20.9 

13 

0.0 

yes 

yes 

21.3 

14 

1.7 

yes 

yes 

20.9 

15 

2.5 

yes 

yes 

21.9 

16 

4.1 

yes 

yes 

21.5 

17 

6.1 

yes 

yes 

20.5 

•  •>  1%  VW  irvl 


CHAPTER  2 


THEORETICAL  FOUNDATIONS 


2 . 1  OVERVIEW 

The  present  study  consists  of  three  sets  of  distinct,  but  closely 
related  experiments.  Each  set  of  experiments  has  its  own  theoretical 
framework.  In  the  following  sections,  the  theoretical  background  for 
each  set  is  outlined;  then  a  comparison  of  the  similarities  and  differ¬ 
ences  among  them  is  given. 


2.2  CASE  I:  WIND  WAVE  IN  EULERIAN  FRAME 
2.2.1  Triple  Decomposition 

Following  Benilov  et  al.  (1974),  one  can  decompose  any  flow  vari¬ 
able  f(x,t),  such  as  velocity,  into  three  components: 

f(x,t)  -  f 00  +  f(x,t)  +  f'(x,t) 


where 


the  mean  value  of  f(x,t) 


f(:c,t)  *  the  wave-induced  component 
f'(jc, t)  *  the  fluctuating  component 

x  -  the  Cartesian  coordinate  vector  (xj,  X2»  X3,  or  x,  y,  z), 
*1  or  x  ■  the  streamwise  distance  measured  with  respect  to  the  air 
inlet 

X2  or  y  *  the  vertical  distance  (positive  upwards)  measured  with 
respect  to  the  still-water  level,  and 
X3  or  z  ■  the  spanwise  distance  measured  with  respect  to  the  cen¬ 
terline  of  the  test  facility. 


The  time  average  is  defined  as: 


T(x)  -  lim  /  f (x,t)  dt 

T-md  o  * 


For  a  sampled  data  system,  the  time  average  takes  the  form: 


(2-1) 


1  N 

f(x)  -  lim  =  l  f(x,t  ) 

M  _ _  W  1  1 


N-m»  i-1 

It  is  also  assumed  that  both  f  and  f*  have  zero  means,  i.e.,  f  =  0  and 
f  '  *  0.  The  statistical  properties  of  f  and  f '  can  be  found  by  the 
wave-separation  method  outlined  in  the  next  section. 

2.2.2  Wave-Separation  Method 

In  order  to  isolate  the  wave  component  of  a  flow  variable,  a  coher¬ 
ent  relationship  is  assumed  between  the  wave-induced  component  and  the 
wave  motion  at  the  interface.  That  is,  the  wave-induced  component 
fOc.t)  is  assumed  to  be  statistically  coherent  with  the  surface  dis¬ 


placement  n(r,t): 


f(x,t)  -  Ln(r,t) 


where  L  is  some  linear  operator  and  r  is  a  horizontal  vector  at  the 
interface.  A  necessary  second  assumption  is  mutual  orthogonality  of  the 
three  components: 

ff  -  0,  ff '  -  0,  ff '  -  0 

Naturally,  n  is  composed  of  wave  components  only,  i.e.,  n  ■  n  and  n  s 
n'  =  0.  In  Appendix  A,  it  is  shown  that  information  about  all  second 
moments  of  f  and  f*  can  be  obtained  without  reconstructing  realizations 
of  f  and  f'.  For  this  study,  only  second  moments  are  of  interest; 
hence,  no  reconstruction  of  f  or  f*  is  required. 


2.2.3  Governing  Equations 


The  governing  equations  in  this  case  are  those  for  an  incompres¬ 


sible  Newtonian  fluid,  viz., 


continuity: 


(2-2) 


2 

3u  3u,  .  .  3  u. 

„  i  ,  i  1  3p  ,  i  ,  _ 

momentum:  -sr-—  +  u  .  -» —  - +  v  -s — * —  +  F . 

3t  j  3Xj  p  3x,3xj  i 

where  *  F-j  -  0,  F2  *  “g*  Let 


(2-3) 


Ui  "  ui  +  ui  +  ui 


(2-4a) 


P  “  P  +  P  +  P' 


(2-4b) 


Substituting  Eqs .  2-4a  and  2-4b  into  Eq.  2-2  and  then  taking  the  time 


average  gives: 


(2-5a) 


<5i  +  UP  '  ° 


(2-5b) 


Substituting  Eqs.  2-Aa  and  2-Ab  into  Eq.  2-3  and  using  the  continuity 
equation  yields: 


3 


+  ut  +  u')  +  [(Uj  +  u.  +  u^KUi  +  UjL  +  up] 


‘  5  1^  +  '  +  p'>  +  v  (ui  +  “1  +  “p  +  Fi 


(2-6) 


It  <U1 +  “1 +  UP  +  fej  <“luj +  ui“j +  uluj 


+  uluj +  Vj +  Vi +  UP) +  “pj +  “pp 


I  fe  +  f  +  »'>  +  u  HIT.  <ui  +  “i  +  UP  +  ri 

i  J  j 


(2-7) 


.*•**-  .  “*  /  •  _  - 


■•>‘vvv':‘va-: 


Time  averaging  Eq.  2-7  produces: 


wi  (Vj +  Vj +  ui“? 


p  3xt  3xj3xj  1 


W~  <ui“j>  ■  ‘  \  1^  +  v  '  1^  (uiuj  +  “i"j>  +  F1 


(2-8) 


Since  no  simple  averaging  process  is  available  for  the  random  wave- 
induced  components,  further  simplification  of  equations  is  not  possible. 
A  different  decomposition  was  given  by  Benilov  and  Zaslavskii  (1974), 
but  such  decomposition  complicates  the  equations  and  makes  comparison 
with  phase-averaging  results  (see  Sec.  2.3.1)  impossible. 

If  it  is  assumed  that  the  mean  flow  and  the  wave-induced  motion  are 
two-dimensional,  then  the  time-averaged  equations  (2- 5a,  2-8)  become: 


continuity: 


momentum: 


“  (u  u  )  +  4—  (u  v  ) 
3x  7  3y 


3u  .  3v 

—  +  -r—  ••  0 

3x  3y 


-  I  iP  +  vf®J«  + 

p  3x  v^.  2  .  2J 

3x  3y 


(2-9) 


(u’u’) 


h  (u’v') ' I;  (“ " >  'I?  > 


(2-10) 


—  2-  2— 

3  , —  —  .  ,  3  , —  —  .  1  3p  ,  /3  v  ,  3  Vi  3  . — ; — f-. 

3x<vu)  '5y<vv)  'p1y  +  ^T)  '  Si  <v  u  > 


3  /  v  t  \  3 

-  -5^  (v’v’)  -  - 


fe  (v  u  )  -  ^  (V  v  )  -  g  (2-11) 


2.2.4  Boundary  Conditions 


The  no-slip  boundary  condition  applies  at  the  channel  floor: 

u^  ■  0  at  y  ”  -H 

where  H  is  the  water  depth.  A  similar  condition  also  applies  at  the 


-  v  .  »  O  -  -  •  „ 


.  •  „*v  „V  \  . 

■  •  •  .  •  .  7  .  - 


'WKv;-. 


.  »  *  *  -  ■  ■  §  «  •* 
t*  '•  »  »  k  *  *.  "  •  ,*  »  "  w  * 

V-V.’.'.'.vV.VAAv 
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air-water  interface  provided  no  (or  negligible)  wave  breaking  occurs: 


u1  *  u,  at  y  “  n 

The  labels  "air”  and  "water"  refer  to  the  respective  velocity  components 
on  each  side  of  the  interface-  The  no-wave-breaking  assumption  is  valid 
for  this  study  with  the  possible  exception  of  some  of  the  high  wind- 
speed  cases.  Other  dynamical  conditions  to  be  satisfied  at  an  interface 
for  the  general  case  of  a  viscous  fluid  with  surface  tension  are: 

1.  continuity  of  the  tangential  stress  as  one  passes  through 
the  interface. 

2.  discontinuity  in  the  normal  stress  proportional  to  the 
mean  curvature  of  the  boundary  surface  due  to  surface 
tension  as  one  passes  through  the  interface. 

Mathematical  formulation  of  these  conditions  can  be  found  in 
Wehausen  and  Lai  tone  (1960). 

2.3  CASE  11:  WIND-RUFFLED  MECHANICALLY  GENERATED  WAVE  IN  EULERIAN  FRAME 
2.3.1  Phase  Average 

One  can  decompose  each  flow  variable  into  three  components  in  a 
manner  similar  to  that  used  for  the  wind-wave  case: 

f(x,t)  -  f(x)  +  f(x,t)  +  f'(*,t) 

Again,  the  time  average  of  f  is  as  defined  in  Eq.  2-1. 

In  order  to  separate  out  the  periodic  component  f,  we  define  a 
phase  average: 

,  N-l 

<f(x,t)>  ■■  lim  ^  £  f(x,t+m)  (2-12) 

N+<*>  n*0 

where  t  is  the  time  period  of  the  imposed  periodic  excitation.  The 
periodic  excitation  in  this  case  is  the  mechanically  generated  wave, 


and  t  ia  the  frequency  of  the  wave.  Of  course,  in  practice,  N  is  always 


finite  but  large. 

Because  f'  is  uncorrelated  to  the  periodic  excitation  over  time, 
the  phase  average  of  f  contains  only  the  mean  and  the  periodic  compo¬ 
nents,  l.e., 

<f(x,t)>  -  f(x)  +  f(x,t) 


The  time  average  of  the  phase  average  is,  therefore,  equal  to  the  mean, 
viz . , 

<f(x,t)>  -  f (x) 


Thus, 

f(jc,t)  -  <f(£,t)>  -  f(x) 
f'(_x,t)  -  f(jc,t)  “  <f(jc,t)> 

The  following  relations  can  be  developed  from  the  above  definitions: 


<f'>  -  0 

fg  -  ?g 
f  -  0 

7’  -  0 


<fg>  -  f<g> 

<7>  -  7 

<fg>  -  f<g> 

Ig7  -  o 


<fg’>  -  o 


2.3.2  Governing  Equations 

Equations  2-5a  and  2-8  still  hold  for  this  case.  By  use  of  the 
definition  of  phase  average,  further  decomposition  of  the  equations  is 
possible.  Taking  the  phase  average  of  Eq.  2-2  and  then  substractlng 
Eq.  2~5a  gives: 

3u.  9u! 

T-i  -  0  ,  ~  •  0  (2-13a,b) 


Again,  phase  averaging  Eq.  2-6  produces: 


¥T  +  lij  (ui°j  +  Vj  +  ”i“j  +  Vj  +  <UM>> 


'  ‘  5  %  <»  +  »>  +  v  3^  <ui  *  “i>  +  ri 


(2-14) 


Subtracting  Eq.  2-8  from  Eq.  2-14  yields: 


sr  *  <  Vj +  v*j +  Vj +  uiu? 


.1*  +v^ 
P  V  8Xj3Xj 


or,  upon  rearrangement, 


3  _  _  _ 

IT  +  a5EJ  (uiuj  +  uiuj> 


1  3?  ,  ..  a2“l 

p  3x^  3Xj3Xj 


^  (vj +  ui“i>  <2"15> 


It  is  again  assumed  that  the  mean  and  the  periodic  flow  fields  are 
two-dimensional,  that  is: 


(u,  v,  w) 


(2-16a) 


(u,  v,  0) 


(2-16b) 


-  (u,  v,  0) 


(2-16c) 


(*,  y,  *) 


(2-16d) 


Using  Eq.  2-16,  one  finds  that  the  time-  and  phase-averaged  continuity 
and  momentum  equations  become: 


3u  3v 

■T—  +  -jr-  *  0 

3x  ay 


(2-17a) 


3u  ,  3v  n 

—  +  -r—  =  0 

3x  3y 


(2-17b) 
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<;  s  >  +  (a  v  )  - 


3  /TT\  3  ,x-,r\  3  >'B~s~s 

-  Tj  <u  v  >  '  3i  (u  u  >  ’  37  (u  v  > 


!i  (v  a  )  +  (v  v  ) 


-  2-  2- 

1  3p  .  /■  3  v  ,  3  v  \  3  . — j — T\ 


■  I7  <v’v’>  ■  h  (v’ G )  ~  h  ^  * ) ' 8  (2 


Ft  +  2  lx  (“  u  }  +  I7  (u  V  +  U  v) 


1  3p  ,  (32u  32u 


T~2  +  7*1  J  ”  lx 

3x  3y 


)  -  -r-  (U  U  +  U'u'  )  “ 


(U  V  +  U'V') 


^  +  (v  G  +  V  U)  +  2  (V  7) 


*  -  —  +  vf— — =  +  — x)  -  4—  (u  V  +  u'v')  ~  -5—  (v  v  +  v'v')  (2 

p  3y  3x2  3y  3x  3y 

2.3.3  Boundary  Conditions 

The  boundary  conditions  stated  in  2.2.4  also  hold  for  this  case 


2.4  CASE  III:  WIND-RUFFLED  MECHANICALLY  GENERATED  WAVE  IN  EULERIAN 
WAVE-FOLLOWING  FRAME 

2.4.1  Transformed  Coordinate  System 


The  mobile  interface  poses  the  biggest  challenge  to  a  satisfactory 
theoretical  as  well  as  experimental  treatment  of  the  problem.  In  par¬ 
ticular,  the  random  nature  of  the  interface  in  the  wind-wave  case  makes 
the  problem  virtually  impossible  to  solve.  However,  for  the  mechani¬ 
cally  generated  wave  cases,  some  further  theoretical  development  is 
feasible  due  to  the  well-defined  periodicity  of  the  main  wave.  One  way 


*.  V  V  ' 

'  -\  V 

■l  .uju  *.  m  r.  *7 . 


to  resolve  the  Interface  problem  for  mechanical-wave  cases  is  to  trans¬ 
form  the  entire  problem  to  a  new  coordinate  system  such  that  both  the 
interface  and  the  lower  boundary  of  the  flow  are  fixed  in  the  new  coor¬ 
dinates  . 

One  can  transform  the  Cartesian  coordinate  system  and  the  flow 
variables  into  a  new  coordinate  system  such  that  the  velocity  components 
are  aligned  with  the  axes  of  the  new  curvilinear  orthogonal  system 
(Benjamin,  1959;  Miles,  1959;  Gent  and  Taylor,  1976).  In  addition,  a 
pair  of  constant  spatial  values  in  this  new  coordinate  coincide  with  the 
channel  floor  and  the  interface.  This  transformation  simplifies  theo¬ 
retical  analysis,  but  leads  to  great  difficulties  in  experimental  real¬ 
ization  because  one  needs  to  constantly  align  the  probe  direction  with 
the  axes  in  order  to  measure  the  correct  velocity  components,  while  at 
the  same  time  maintaining  the  probe  at  a  constant  spatial  location  in 
the  transformed  system.  Therefore,  many  theoretical  works  using  this 
type  of  transformation  remain  experimentally  untested.  Another  approach 
is  to  transform  only  the  coordinate  system,  but  not  the  flow.  That  is, 
the  Cartesian  flow  field  is  interpreted  in  a  curvilinear  coordinate 
system  (Norris  and  Reynolds,  1975;  Hsu  et  al.,  1981).  Experimental 
realizations  can  be  achieved  quite  easily  in  this  system. 

One  form  of  such  a  coordinate  system  is  given  as  follows: 


t  -  t*  ;  x  -  x* 

y  ■  y*  +  f (y*)n  ;  z  *  z* 


n  ■  a  cos(kx*  -  u>t*) 


f(y*) 


slnh(kH  +  ky*) 
s i nh(kH) 


(2-22) 


19 


where  (x,  y,  z)  denotes  spatial  points  in  the  Cartesian  coordinate 
system  and  (x*,  y*,  z*)  denotes  spatial  points  in  the  new  coordinate 
system.  The  amplitude  of  the  imposed  mechanical  wave  is  a.  The  trans¬ 
formation  between  (x,  y,  z)  and  (x*,  y*,  z*)  is  one-to-one,  but  non- 
orthogonal;  however,  it  has  the  following  desirable  properties: 

y*  ■  0  <«>  y  *  n 
y*  -  -H  <->  y  =*  -H 


Experimental  realization  of  this  coordinate  system  can  be  achieved  by 
oscillating  the  probe  vertically  at  a  fixed  fetch  so  the  probe  lies  at 
selected  (x*,  y*,  z*)  *  constant  points.  The  particular  function 

(Eq.  2-22)  is  chosen  because  y*  ■  constant  surfaces  are  streamlines  of 
an  inviscid  flow  below  a  monochromatic,  progressive  wavy  surface. 

2.4.2  Governing  Equations  in  the  Transformed  Coordinate  System 

The  governing  equations  in  (x,  y,  z,  t)  coordinates  can  be  trans¬ 
formed  into  (x*,  y*,  z*,  t*)  coordinates  using  the  chain  rule  for  par¬ 
tial  derivatives.  This  transformation  is  similar  to  that  of  Hsu  et  al. 
(1981),  but  all  the  higher  order  terms  are  retained.  Then  substituting 
the  decomposition  (Eq.  2-4)  into  the  transformed  equations  and  taking 
time  averages,  we  obtain  the  two-dimensional  mean  flow  equations  in  the 
(x*,  y*,  z*,  t*)  coordinates,  for  example, 
continuity: 


[1/(1  -  i  ^  iij  +  »4/(i  ♦  iij  S)  .  o  (2-23) 

ax*  ay*l/v  ay*  n  ax*  ay*  ay  /  ay 


If  the  phase  averages  are  taken  before  the  time  averages,  the  difference 
between  the  phase-  and  time-averaged  equations  yields  the  wave- 
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perturbation  equations,  for  example, 


continuity: 


3u  3q  3u  3n  3u 

3x*  3x*  3y*  3x*  3y 

Because  y*  ■  0  corresponds  to  y  -  q,  so  provided  n*  «  q,  any 
position  (x*,  y*,  2*)  with  y*  <  0  is  always  beneath  the  interface. 
Thus,  these  equations  are  valid  for  the  entire  water  regime. 

2.4.3  Boundary  Conditions  in  the  Transformed  Coordinate  System 

In  the  transformed  coordinates,  the  no-slip  boundary  conditions 
become : 


3(v  +  v 

IT 


* 


1  + 


3f 

a? 


n) 


(2-24) 


and 


at  y*  =  -H 


u.  =  u(  at  y*  »  0  provided  q'  «  q 

air  1water 

The  conditions  for  the  stresses  at  the  interface  given  in 
Sec.  2.2.4  are  also  valid  for  this  case. 


2.5  COMPARISON  OF  THE  THEORIES 

According  to  the  wave-separation  method  (Sec.  2.2.2),  the  surface 
elevation  q  i3  defined  to  be  composed  entirely  from  wave  motions  and  the 
fluctuating  component  q'  =  0.  For  the  other  two  cases  (II  and  III), 
the  deviation  of  the  water  surface  from  the  periodic  part  n  is  the 
fluctuating  component  n'.  This  is  the  fundamental  difference  in  the 
definition  of  wave  and  fluctuating  components  for  the  surface  elevation 
between  Case  I  and  Cases  II  and  III.  It  is  possible  to  use  both  the 
phase-averaging  and  wave-separation  methods  in  Case  II.  For  no  wind  and 
low  wind-speed  cases,  q'  <<  q,  so  the  results  given  by  the  wave- 
separation  method  should  be  good  approximations  to  those  obtained  by 
phase  averaging. 
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We  believe  that  Che  phase-averaging  results  are  more  accurate 
because  linearity  is  not  assumed  in  the  separation  process.  Hence,  an 
error  estimate  can  be  made  on  the  wave-separation  method  by  comparing 
the  results  given  by  these  two  different  schemes  (see  Sec.  3. 5. A). 
Although  it  can  be  shown  that  the  theoretical  error  is  small  for  the 
linear  wave-separation  method  (see  Benilov  et  al.,  197A;  Benilov,  1978), 
an  error  higher  than  the  theoretical  limit  is  likely  due  to  finite  reso¬ 
lution  of  spectra  and  numerical  integration  of  the  spectra  in  the  wave- 
separation  method. 

The  main  difference  between  Cases  II  and  III  is  the  coordinate 
system.  In  the  (x,  y,  z,  t)  system,  u(x,  y,  t)  and  v(x,  y,  t)  are 
Cartesian  velocities  that  describe  actual  water  motion  about  the  mean 
velocities  u(x,  y)  and  v(x,  y)  at  fixed  x  and  elevation  y.  On  the  other 
hand,  u(x*,  y*,  t*)  and  v(x*,  y*,  t*)  are  also  Cartesian  velocities 
which  represent  the  periodic  flow  field  about  the  mean  velocities  u(x*, 
y*)  and  v(x*,  y*)  at  constant  x*  and  y*.  Also  in  Cases  I  and  II,  quan¬ 
tities  such  as  u,  u,  etc.,  are  not  well  defined  in  the  domain  |y|  <  |n|. 
In  this  region,  a  given  Eulerian  space  point  is  sometimes  in  the  air 
and  sometimes  beneath  the  water.  However,  all  statistical  quantities 
are  well  defined  in  the  transformed  coordinate  system  in  Case  III  for 
y*  <  0,  provided  n'  «  n. 

The  time-averaged  governing  equations  are  the  same  for  Cases  I  and 
II.  Then,  using  the  phase  average,  the  periodic  flow  field  can  be  sepa¬ 
rated  from  the  mean  and  turbulent  fields  in  Cases  II  and  III.  However, 
the  time-averaged  and  phase-averaged  equations  in  Case  III  are  quite 
different  from  those  of  Case  II.  Although  many  of  the  terms  in  the 
fixed-frame  equations  can  readily  be  identified  in  the  wave-following 


frame  equations — this  can  be  seen  by  comparing  Eqs.  2-17a,b  with 
Eqs.  2-23  and  2-24 — there  are  many  extra  inhomogeneous  terms  added  to 
the  equations  in  Case  III.  Hsu  et  al.  (1981)  discarded  many  of  these 
inhomogeneous  terms  because  they  are  higher-order  wave-perturbation 
quantities  and  are  small  compared  to  the  air  velocity.  In  this  study, 
the  wave-perturbation  quantities  are  large  compared  to  the  mean  flow,  so 
the  higher-order  terms  cannot  be  neglected  without  further  justifica¬ 
tions  . 


CHAPTER  3 


EXPERIMENTAL  FACILITIES  AND  INSTRUMENTATION 

The  experimental  facilities  and  some  of  the  specially  designed 
instruments  are  described  in  this  section.  The  results  of  the  qualify¬ 
ing  procedures  for  the  instruments  are  also  included  here.  Special 
methods  or  techniques  used  during  data  acquisition  and  reduction  are 
outlined  in  Secs.  3.4  and  3.3. 

3.1  THE  WIND,  WATER-WAVE  CHANNEL 

The  channel  used  for  the  experiments  was  described  in  detail  by  Hsu 
(1965).  This  facility  has  cross-sectional  dimensions  of  1.93  m  high  by 
0.91  m  wide.  The  overall  length  is  35  m  with  a  22-m  long  glass  test 
section.  The  cross-sectional  area  varies  by  about  IX  along  the  test 
section.  For  the  current  set  of  experiments,  the  water  depth  in  the 
channel  was  0.965  m  leaving  0.965  m  for  the  vertical  depth  of  the  air 
flow.  A  schematic  view  of  the  facility  is  shown  in  Fig.  1. 

Situated  at  the  upstream  end  of  the  channel  is  a  horizontal- 
displacement-type  mechanical  wave-generating  plate  which  is  actuated  by 
a  hydraulic  cylinder.  The  motion  of  the  hydraulic  cylinder  is  con¬ 
trolled  by  a  closed-loop  servo-system.  The  signal  from  a  function 
generator,  with  a  frequency  resolution  of  i  0.05Z,  forms  the  input  to 
the  servo-system  for  the  generation  of  a  progressive  wave  train.  At  the 
downstream  end  of  the  channel  la  a  beach  which  consists  of  stainless 
steel  turnings  packed  in  expanded-metal  baskets  mounted  upon  a  1  on  5 
sloping  frame.  The  reflection  coefficient  for  waves  at  a  frequency  of 
1  Hz  is  found  to  be  less  than  5X  for  all  the  test  cases.  The  reflection 
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coefficient  was  expected  to  be  less  for  wind  waves  because  of  smaller 
amplitude  and  higher  frequency. 

Just  above  and  downstream  of  the  beach  Is  a  centrifugal,  alrfoil- 
bladed  fan  which  draws  air  through  the  test  section.  An  air¬ 
straightening  vane  consisting  of  a  100-mm  thick  stainless  steel  honey¬ 
comb  section  prevents  the  vortex  motion  created  by  the  fan  from  extend¬ 
ing  upstream  into  the  channel  test  section. 

The  air  inlet  is  located  4.92  m  downstream  of  the  mean  position  of 
the  wave  generating  plate,  thus  allowing  the  mechanically  generated 
water  waves  to  develop  before  being  acted  on  by  wind.  Fiberglass  fur¬ 
nace  filters  are  mounted  at  the  entrance  to  the  air  inlet  to  prevent 
large  particles  of  dust  from  getting  into  the  channel;  a  100-mm  thick 
honeycomb  section  and  fine-meshed  screens  are  also  installed  at  the  exit 
of  the  air  inlet  to  produce  a  uniform  turbulent  velocity  profile  in  the 
test  section.  However,  there  was  no  transition  plate  and  the  basic  air 
inlet  conditions  were  similar  to  those  employed  by  Hsu  (1977)  and  Chen 
(1981),  except  for  small  differences  in  the  water  depth. 

3.2  PROBE  SUPPORT  MECHANISM 

A  motorized  traversing  and  position  indicator  (hereafter  called  the 
traverse)  permitted  control  of  probe  location  in  the  vertical  and  span- 
wise  planes  within  the  wind-wave  facility.  A  detailed  description  of 
the  traverse  is  given  in  a  report  by  Chambers  et  al.  (1970). 

The  traverse  supports  an  air  free-stream  temperature  probe,  a  bulk 
water  temperature  probe,  a  pitot-static  tube  and  wave-height  gauge(s) 
(WHG).  A  typical  configuration  of  the  probe  setup  is  shown  in  Fig.  2. 

For  the  wind-wave  experiments,  the  wave-height  gauge  was  mounted 
right  next  to  the  laser  Doppler  anemometer  (LDA)  probe  volume.  However, 
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for  the  mechanically  generated  wave  experiments,  the  plane  of  the  two 
wave-height  gauges  was  offset  0.21  m  to  one  side  of  the  channel  center 
line  (about  midway  between  the  side  wall  and  the  center  line  of  the 
channel),  but  the  sensing  wire  of  the  first  wave-height  gauge  remained 
at  the  same  fetch  as  the  LDA  probe  volume*  The  LDA  probe  volume,  the 
wave-height  gauge  sensing  wire,  the  pitot-static  tube,  and  the  air  free- 
stream  temperature  probe  were  at  a  fetch  of  13.0  m  from  the  air-inlet. 
The  bulk,  water  temperature  probe  was  mounted  0.34  m  behind  the  LDA  probe 
volume  and  0.35  m  below  the  mean  water  level  surface. 

The  pitot-static  tube  was  mounted  0.70  m  ahead  of  the  main  support 
of  the  traverse.  This  distance  was  the  shortest  distance  ahead  of  the 
traverse  at  which  the  probes  could  be  mounted  and  still  experience 
minimal  interference  effects  due  to  the  traverse  (Chambers,  1970).  For 
wind-wave  experiments,  the  flow  interference  effect  of  the  sensing  wire 
and  support  of  the  WHG  on  the  water  velocity  was  tested  by  measuring 
velocity  with  and  without  the  gauge  next  to  the  LDA  probe  volume.  No 
noticeable  differences  were  found  for  the  mean  and  fluctuating  velocity 
components.  Indeed,  the  wire  diameter  was  only  0.15  mm  which  was  much 
smaller  than  the  length  of  the  LDA  probe  volume  (~  3  mm),  and  the  sup¬ 
port  was  150  mm  behind  the  measuring  point.  However,  for  mechanically 
generated  waves,  the  lower  support  of  the  WHG  frame  was  found  to  produce 
some  interference  effect  on  velocity  measurements  due  to  the  slower 
decay  with  depth  of  the  longer  wave.  Hence,  the  wave-height  gauges  were 
moved  off  the  center  line  for  all  the  mechanically  generated  wave  cases. 
The  bulk  water  temperature  probe  was  found  to  have  no  influence  on 
velocity  measurement  for  it  was  quite  a  distance  behind  and  below  the 
LDA  probe  volume. 
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3.3  INSTRUMENTATION  AND  QUALIFICATION  PROCEDURES 

3.3.1  Wave  Height 

Capacitance  wave-height  gauges  were  used  to  measure  the  water 
surface  elevation.  Two  gauges  were  used  for  all  mechanically  generated 
wave  cases  in  order  to  resolve  various  components  of  the  water-wave 
train.  However,  only  one  gauge  was  used  for  wind-wave  runs.  Each  gauge 
is  constructed  with  #36  Nylclad  insulated  copper  conductor  which  is 
stretched  between  two  ends  of  a  U-shaped  frame  and  connected  In  parallel 
to  one  arm  of  a  Sanborn  958-1100  capacitance  bridge  amplifier.  The 
gauge  (or  gauges)  was  aligned  with  the  wire  normal  to  the  mean  water 
level  and  with  half  its  length  Immersed  in  the  water.  A  thorough  dis¬ 
cussion  of  the  characteristics  of  the  capacitance-wire  gauge  and  its 
associated  electronics  was  given  by  Colonell  (1966). 

The  capacitance  gauge  was  calibrated  before  and  after  each  experi¬ 
ment  for  the  wind  wave  runs  because  different  gain  settings  were  used 
for  the  various  flow  cases.  For  the  mechanically  generated  wave  cases, 
the  gauges  were  calibrated  before  and  after  a  series  of  experiments. 
When  the  probes  were  cleaned  periodically,  the  change  in  calibration  was 
found  to  be  less  than  0.5%  over  a  period  of  2  months. 

Static  calibration  was  performed  by  lowering  and  raising  the  probe 
support  mechanism  so  as  to  change  the  immersion  depth  of  the  wire.  The 
immersion  depths  were  read  off  the  probe  support  mechanism  (readable  to 
±  0.25  mm)  with  the  corresponding  voltages  recorded  by  the  data  acquisi¬ 
tion  system  (Sec.  3.3.9). 

3.3.2  Air  Free-Stream  Velocity 

The  mean  air  free-stream  velocity  or  the  reference  wind  speed  was 
obtained  from  the  pressure  differential  of  a  pitot-static  tube.  The 
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pitot-static  tube,  which  has  a  3.0  mm  diataeter,  is  a  Model  PAC-12-KL 
manufactured  by  United  Sensors  and  Control  Corporation. 

A  Pace  differential-pressure  transducer  (Model  P90D)  with  a  full 
range  of  ±  21.2  mm  of  water  was  used  to  measure  the  pressure  differen¬ 
tial  from  the  pitot-static  tube.  The  Pace  transducer  was  connected  as 
one  leg  of  an  Inductive  bridge  of  a  Sanborn  656-1100  carrier  amplifier. 
The  transducer  amplifier  system  was  calibrated  for  several  attenuation 
scales  using  a  Combust  Institute  micromanometer  (with  resolution  of 
±  0.006  mm  of  fluid  of  specific  gravity  0.82).  The  calibration  curves 
were  checked  from  time  to  time  during  experimental  runs  and  the  maximum 
deviation  was  2%. 

3.3.3  Temperature 

Mean  air  free-stream  temperature  and  mean  bulk  water  temperature 
were  measured  using  Victory  Engineering  43A38105  bead  in  glass  thermis¬ 
tor  probes  in  conjunction  with  bridges.  Calibration  was  performed  in  a 
Rosemont  constant  temperature  bath;  the  temperature  standard  was  a 
Hewlett  Packard  2801A  quartz  thermometer.  Parabolas  were  least-squares 
fitted  to  the  calibration  points  with  typical  standard  deviations  of 
±  0.1°C  over  a  temperature  range  of  8°C-28°C  and  a  time  span  of  10  months. 

3.3.4  Surface  Drift  Velocity 

The  Lagranglan  surface  drift  was  determined  by  measuring  the  time 
required  for  paper  punchings  to  travel  a  fixed  distance  (usually  0.5  m 
or  1.0  m)  immediately  upstream  of  the  measuring  station.  To  obtain  a 
single  value  of  Lagranglan  surface  drift  for  one  day's  experiment,  approx¬ 
imately  five  measurements  were  made  before  the  experiment  was  started  and 
five  during  the  experiment.  The  mean  time  of  travel  for  all  the  trials 
was  computed  to  obtain  the  surface  drift  velocity  for  the  experiment. 
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For  wind-wave  experiments,  the  standard  deviations  for  the  mean 
Lagrangian  surface  drifts  were  about  10Z  for  wind  speed  less  than  5  m/s. 
At  higher  wind  speeds,  standard  deviations  were  as  large  as  13Z.  This 
was  due  mainly  to  actual  fluctuations  and  timing  errors  (travel  times  at 
the  high  wind  speeds  were  on  the  order  of  two  to  three  seconds). 

For  mechanically  generated  wave  cases,  the  main  wave-induced  peri¬ 
odic  motion  to  the  floats  and  timing  errors  could  possibly  be  on  the 
order  of  one  second  (the  period  of  the  main  wave).  Fortunately,  the 
range  of  wind  speeds  used  for  mechanical  wave  experiments  was  low  and 
the  percentage  error  introduced  in  the  float  transit  timing  would  not  be 
too  high.  The  standard  deviations  of  the  mean  Lagrangian  surface  velo¬ 
city  were  found  to  be  about  71  of  the  mean  values. 

3.3.5  Distance  from  the  Interface 

The  LDA  probe  volume  position  was  read  from  a  dial  height  gauge 
(readable  to  ±  0.01  mm)  or  from  a  reference  voltage,  which  is  linearly 
related  to  the  displacement,  on  the  LDA  support  table.  The  gauge  was 
zeroed  at  still-water  level  by  moving  the  table  so  that  the  plane  of  the 
blue  and  the  green  laser  beams  of  the  LDA  coincide  with  the  still-water 
level  surface.  Some  water  level  change  over  the  course  of  an  experiment 
was  unavoidable.  A  linear  trend  with  time  was  assumed  and  taken  into 
account  in  establishing  the  actual  probe  volume  location.  In  addition, 
the  water  set-up  due  to  wind  was  monitored  by  four  manometers  along  the 
channel.  The  set-up  at  the  measuring  location  was  negligible  for  wind 
speeds  lower  than  6  m/s.  For  higher  wind  speeds,  depth  corrections  were 
necessary  due  to  set-up.  The  overall  accuracy  of  the  corrections  was 


estimated  to  within  ±  0.5  nm. 


3.3.6  Water  Velocit 


Instantaneous  values  of  the  horizontal  (u)  and  vertical  (v)  velo¬ 
city  components  were  measured  with  a  two-component  DISA  53X  modular 
optical  laser  Doppler  anemometer  (LDA)  system.  The  system  Is  equipped 
with  Bragg  cell  frequency  shift  and  Is  capable  of  operating  in  both  for¬ 
ward-scatter  and  backward-scatter  modes.  Forward-scatter  mode  was  used 
for  all  fixed-frame  experiments,  while  backward-scatter  mode  was  used  In 
wave-following  experiments. 

The  laser  used  for  the  LDA  was  a  Spectra  Physics  Model  164-06 
argon-ion  type.  A  600  mm  focal  length  achromatic  focusing  lens  was 
used,  giving  half-intersection  angles  of  2.47°  and  2.53°  for  the  LDA 
beams  (see  Fig.  3).  The  beams  were  oriented  such  that  the  plane  con¬ 
taining  the  green  and  the  blue  beams  was  always  parallel  to  the  still- 
water  level  with  the  optical  axis  of  the  lens  perpendicular  to  the 
channel  sidewall.  The  forward-scatter  optics  were  normally  placed  at  a 
10°  angle  to  the  optical  axis  of  the  system  on  the  opposite  side  of  the 
channel.  The  probe  volume  in  the  forward-scatter  mode  was  about  0.15  mm 
in  diameter  and  3.0  mm  long.  The  probe  volume  in  the  backward-scatter 
mode  was  about  0.15  mm  in  diameter  and  3.5  mm  long. 

Signal  processing  was  accomplished  using  DISA  type  55N20  frequency 
trackers  and  DISA  type  55L72  range  translators.  The  down  shifted 
Doppler  signals  from  the  range  translators  were  band-pass  filtered  using 
Krohn-Hite  3100  filters.  The  signal  fed  into  the  trackers  was  thus 
enhanced  as  the  level  of  extraneous  noise  was  reduced.  However,  care 
was  taken  so  as  not  to  filter  out  relevant  frequencies.  For  wind-wave 
experiments,  the  filter  bandwidth  ranged  from  40  kHz  (centered  at  the 
zero  velocity  frequency)  for  the  lowest  wind  speed  to  140  kHz  for  the 
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highest  wind  speed.  During  mechanical-wave  experiments,  the  filters 
were  set  at  bandwidths  of  120  kHz. 


The  linear  frequency  to  voltage  relationships  of  the  frequency 
trackers  were  found  by  two  point  calibrations.  This  was  done  periodic¬ 
ally  during  each  experiment  by  using  the  range  translators  for  known 
frequencies  and  recording  the  output  voltages  (after  the  signal  condi¬ 
tioning  electronics)  on  the  data  acquisition  system  (see  Sec.  3.3.9). 
Hence,  we  accounted  for  all  the  electronic  components  (filters,  amplifi¬ 
ers)  in  the  calibrations.  The  relationships  were  linear  and  drifts  were 
less  than  ±  0.5  mm/s  in  each  experiment. 

Given  the  frequency,  one  can  compute  the  velocity  from  the  well- 
known  relationship: 

Vi 
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blue  or  green  LDA  system 
Doppler  frequency  of  the  i-th  system 
laser  wavelength  of  the  i-th  system 
half-intersection  angle  of  the  i-th  system 
velocity  component  of  the  i-th  system. 


The  velocity  components  in  the  Cartesian  coordinate  system  can  be  found 
by  combining  vectorially  the  velocity  components  of  the  blue  and  green 
systems. 

In  addition  to  the  qualification  tests  for  the  scanner-LDA  system 
mentioned  in  the  following  section,  another  qualification  test  was  run 
to  estimate  the  noise  level  when  the  LDA  was  operated  in  a  simultaneous 
forward-scatter  and  backward-scatter  mode  with  the  optical  scanner 


(Sec.  3.3.7).  The  set-up  for  the  test  was  similar  to  that  of  Cheung  and 
Koseff  (1983),  except  that  the  scanner  was  mounted  In  place  and  the  test 
flow  was  generated  by  the  wind-wave  facility.  The  scanning  mirror  was 
not  moving  during  this  test  because  the  forward-scatter  optics  was 
always  operated  in  the  fixed-frame  mode.  However,  the  back-scattered 
light  was  collected  through  the  mirror  system.  Two  wind  speeds  were 
used  for  the  test:  0.0  m/s  and  4.1  m/s. 

The  results  for  the  test  are  listed  in  Tables  2  and  3.  In  the 
cases  of  the  mean  and  of  the  peak  amplitude  of  the  phase-averaged  velo¬ 
cities,  the  forward-scatter  and  backward-scatter  results  differed  by 
about  1  mm/s.  This  difference  is  most  likely  due  to  small  errors  in  the 
calibration  curves  of  the  trackers.  The  root-mean-squared  (RMS)  values 
were  higher  for  the  backward-scatter  mode  due  to  a  larger  probe  volume 
and  higher  background  noise  level  (see  Cheung  and  Koseff,  1983).  If  the 
true  RMS  values  were  obtained  by  cross-correlating  the  two  simultaneous 
signals,  then  the  noise  equivalent  RMS  velocities  are  2-3  mm/s  for  the 
backward-scatter  mode  and  1  mm/s  for  the  forward-scatter  mode.  A  bonus 
result  of  this  test  was  the  establishment  of  the  stability  of  the  scan¬ 
ner  when  there  was  a  constant  input  voltage.  The  forward-scatter  signal 
did  not  deteriorate  during  each  run  (10  min),  indicating  that  the  probe 
volume  stayed  within  the  field  of  view  of  the  forward-scatter  receiving 
optics  (about  0.5  mm  at  the  focal  point).  Therefore,  the  scanner  is 
able  to  hold  the  probe  volume  at  a  fixed  point  to  at  least  ±  0.5  mm. 
3.3.7  Wave  Follower 

A  sketch  of  the  LDA  beams  steering  mechanism  is  shown  in  Fig.  4. 
It  consists  mainly  of  a  fixed  front-surfaced  mirror,  an  optical  scanner 
(General  Scanning  G-300  POT)  with  a  front-surfaced  mirror,  and  an 


electronic  control  unit  for  the  scanner  (General  Scanning  CCX-101-T).  A 


close-up  picture  of  the  scanner  and  a  perspective  view  of  the  LDA- 
scanner  system  in  relation  to  the  water  channel  are  shown  in  Figs.  5  and 
6. 

The  angle  of  rotation  of  the  scanner  mirror,  and  hence  the  actual 
position  of  the  LDA  probe  volume,  can  be  found  from  the  position  trans¬ 
ducer  of  the  scanner.  The  probe  volume  actually  traces  out  an  arc 
rather  than  a  vertical  line;  however,  the  maximum  scan  angle  remains 
small  (6.2°).  Thus,  the  arc  is  an  excellent  approximation  to  the 
assumed  vertical  motion.  Deviation  of  the  scan  trajectory  from  the 
vertical  was  1.5  mm  at  the  most.  This  distance  is  of  the  order  of,  but 
less  than,  the  length  of  the  LDA  probe  volume. 

The  input  to  the  scanner  control  unit  was  a  narrowly  band-pass- 
filtered  signal  from  a  wave-height  gauge  situated  at  the  same  fetch  as 
the  LDA  probe  volume.  The  band-pass  filter  was  set  at  1.00  Hz  with 
24  db/octave  attenuation  slope  and  the  phase  lag  was  zero  at  1.00  Hz.  A 
variable  attenuator  was  used  to  adjuBt  the  overall  gain  of  the  WHG- 
scanner  system  so  that  the  LDA  probe  volume  moved  in  synchronism  with 
the  mechanically  generated  waves. 

Because  the  laser  beams  were  transmitted  through  a  series  of  media 
(air,  glass,  water),  a  computer  program  was  written  to  simulate  the 
actual  motion  of  the  beams  and  to  compute  the  change  in  intersecting 
angles  among  the  beams.  The  effects  on  the  intersecting  angles  were 
found  to  be  minimal,  and  the  moving  mirror  did  not  induce  any  Doppler 
shift  to  the  outgoing  or  Incoming  light.  A  detailed  account  of  the 
analysis  can  be  found  In  Appendix  B. 
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Two  qualification  tests  were  run  to  verify  the  analysis  given  in 
Appendix  B  and  to  establish  the  relationship  between  the  position  trans¬ 
ducer  voltage  output  and  the  actual  distance  scanned  by  the  probe  volume. 
In  the  qualification  tests,  the  wave  follower  and  LDA  system  were  made 
operational  as  if  they  were  being  used  in  an  actual  experiment.  However, 
the  driving  signal  was  a  1.00  Hz  sine  wave  from  a  frequency  generator. 

In  the  first  qualification  test,  the  probe  volume  was  first  made  to 
scan  on  a  piece  of  graph  paper  (graduated  in  mm).  The  LDA  measured  the 
velocities  sensed  by  the  probe  volume  scanning  on  the  graph  paper.  The 
LDA  laser  beams  were  in  the  air  for  this  solid-surface  scan  test.  The 
approximate  scan  distance  was  also  read  from  the  graph  paper  (±  0.5  mm). 
Suppose  the  position  transducer  output  voltage  was  VpQt  *  VQ*sin  2nft 
which  corresponded  to  some  distance  D  ■  kVQ»sin  2nft  *  D0*sin  2nft 
scanned  by  the  LDA  probe  volume.  Then  the  LDA  measured  D  -  kV0(2nf)* 
cos  2nf t  ■  2xfDQ.cos  2itft.  From  a  data  analysis  program,  maxD  (=  2nfDQ) 
and  VQ  were  calculated  via  a  cosine  fitting  procedure  (see  Appendix  C) . 
Since  kVQ  ■  maxD/2xf,  k  can  be  found. 

The  second  qualification  test  was  done  with  the  probe  volume  inside 
the  channel  and  with  the  laser  beams  going  through  the  air-glass-water 
Interfaces  to  test  the  equations  of  the  model  given  in  Appendix  B.  Here, 
the  water  was  assumed  to  be  still  as  the  probe  volume  oscillated.  The 
constant  k  established  in  the  first  test  was  4.269  mm/V  and  4.254  mm/V 
in  the  second  test,  the  difference  being  well  within  experimental  uncer¬ 
tainties. 

In  actual  experimental  runs,  the  measured  velocity  (yM)  is  the 
velocity  of  the  fluid  (y)  relative  to  the  moving  probe  volume,  i.e., 
Xm  *  Jt  “  Xyp  where  Hyp  represents  the  wave-follower  scan  velocity.  The 
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scan  velocity  jjyp  can  be  found  by  differentiating  the  scanner  position 
voltage  and  then  multiplying  the  result  by  the  constant  k.  Next,  the 
scanning  test  data  were  reanalyzed  with  velocity  correction  so  that  the 
actual  Eulerian  velocity  (y)  could  be  found.  Obviously  in  these  tests, 
the  Eulerian  velocity  was  identically  zero  and  all  the  averages  should 
vanish.  Thus,  this  exercise  told  us  how  good  the  velocity  compensation 
scheme  was;  furthermore ,  any  non-zero  averages  would  indicate  the  noise 
level  and  the  limits  of  the  wave-following  LDA  system  in  measuring 
turbulent  quantities. 

For  the  solid  surface  scan  test,  the  mean  u  and  v  were  at  most 
1  mm/s  (see  Table  4).  The  mean  velocities  were  slightly  larger  for  the 
scan  test  inside  the  channel  (Table  5);  this  probably  was  due  to  some 
small  non-zero  mean  velocities  in  such  a  large  body  of  water,  i.e.,  the 
water  was  not  truly  "still."  The  phase-average  results  for  u  and  v, 
after  compensation  of  the  vertical  velocity,  were  found  to  be  at  most 
±  1  mm/s  for  u  and  ±  2  mm/s  for  v;  typical  figures  were  about  ±  0.5  mm/s 
or  less  for  both  tests.  Other  "turbulent"  quantities  are  also  shown  in 
Tables  4  and  5.  Note  that  the  "turbulent"  quantities  have  the  same 
numerical  values  either  with  or  without  vertical  correction  because  of 
the  sinusoidal  driving  signal.  In  the  second  test  (Table  5),  the  noise 
equivalent  values  for  the  mean  velocities  and  the  phase-averaged  veloci¬ 
ties  were  of  the  order  of  1  to  2  mm/s,  about  2  to  5  mm/s  for  the  root- 
mean-square  turbulent  quantities,  and  about  2  to  20  (mm/s)^  for  the  tur¬ 
bulent  stresses.  The  noise  level  was  higher  when  the  wave-follower  scan 
amplitude  (max  D)  exceeded  22  mm  (see  Table  5)  because  the  laser  beams 
were  near  the  point  where  they  failed  to  intersect  inside  the  channel . 
In  actual  experiments,  the  scan  amplitude  was  always  less  than  22  mm. 
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Therefore,  the  wave-following  LDA  system  was  shown  to  be  capable  of  mea-  j 

< 

suring  mean  and  turbulent  quantities  in  the  test  facility.  < 

* 

The  system  transfer  function  of  the  optical  scanner  and  its  control  j 

electronics  was  found  by  monitoring  the  input  and  output  voltages.  The  j 

phase  lag  of  the  optical  scanner  system  at  1.00  Hz  was  found  to  be 
essentially  zero  (less  than  0.5°). 

In  the  light  of  the  above,  the  two  constants  used  in  reducing  data  ! 

t 

from  the  experiments  reported  in  Sec.  4.3  were:  < 

1)  scan  distance  versus  position  output  voltage:  4.254  mm/V;  and 

2)  scan  distance  versus  input  voltage:  129.7  mm/V  at  1.00  Hz. 

3.3.8  Signal  Conditioning  Electronics 

The  low-level  analog  outputs  of  the  wave-height  gauge(s)  and  Pace 
crystal  transducer  for  the  pitot-static  tube  were  amplified  by  Hewlett 
Packard  2470B  data  amplifiers  at  gains  of  twenty.  This  was  done  to 
minimize  quantization  errors  inherent  to  any  analog-to-digital  data 

I 

acquisition  system.  I 

The  LDA  tracker  outputs,  the  wave-height  gauge  output(s),  the  Pace 

I 

crystal  transducer  output,  and  the  position  output  of  the  optical  scan¬ 
ner  were  low-pass  filtered  at  40  Hz  (-3db  point)  by  third-order  Butter- 
worth  anti-aliasing  active  filters.  The  maximum  phase  difference  among 
different  channels  of  the  filters  was  found  to  be  less  than  one  degree 
over  the  0  to  40  Hz  frequency  range.  Therefore,  the  phase  characteris¬ 
tics  were  virtually  identical  among  all  the  filter  channels.  This 
property  is  crucial  In  comparing  phases  among  different  signals  from  the 
filters.  Small  adjustments  in  the  gains  of  the  filters  were  made  to 
allow  use  of  the  full  dynamic  range  of  the  data  acquisition  system. 
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3.3.9  Data  Acquisition/Rcduction  System 

The  data  acquisition  system  consists  of  four  major  components:  a 
minicomputer  (Hewlett  Packard  2100a),  a  digital  magnetic  tape  drive 
(Hewlett  Packard  7970B),  a  hard  disc  drive  (Hewlett  Packard  7900),  and 
an  analog-to-digital  subsystem  (a  modified  Hewlett  Packard  2313  with 
simultaneous  sample-and-hold  capability  on  32  analog  input  channels). 
The  analog  voltages  from  the  instruments  were  sampled  and  digitized  by 
the  analog-to-digital  subsystem,  and  then  recorded  on  digital  magnetic 
tapes  for  future  analysis.  Each  analog  input  has  a  dynamic  range  of 
±  10.230  V  and  a  resolution  of  i  5  mV.  This  system  was  documented  by 
Takeuc'ni  and  Mogel  (1975). 

Data  reduction  was  performed  on  a  Hewlett  Packard  1000  Model  45 
computer  system.  In  addition  to  the  fast  F-series  computer,  the  system 
has  full  graphics  capability  and  a  library  of  data  reduction  programs 
tailored  to  the  experiments. 

3.4  DATA  ACQUISITION  PROCEDURES  AND  TECHNIQUES 

A  flow  chart  depicting  signal  paths  is  shown  in  Fig.  7.  The  chart 
is  for  the  case  of  wind-ruffled  mechanically  generated  waves  with  Eule- 
rian  wave-following.  The  second  wave-height  gauge,  the  wave-follower, 
and  the  sine  wave  reference  signal  were  not  present  for  the  wind-wave 
case.  Also,  the  wave-follower  was  not  used  for  the  mechanical-wave 
fixed-frame  experiments. 

The  sampling  frequency  was  100  Hz  for  the  wind-wave  case  and  200  Hz 
for  the  others.  A  higher  sampling  rate  was  used  for  mechanical-wave 
runs  to  give  better  resolution  of  the  main  wave  in  phase  averages.  The 
sampling  period  for  each  data  point  was  kept  at  10.24  minutes  despite 
the  difference  in  sampling  rate.  For  mechanical-wave  cases,  this  record 


length  was  sufficient  to  resolve  the  anticipated  lowest  frequency  compo¬ 
nent  and  accommodate  more  than  614  periods  of  the  1  Hz  generated  wave. 
A  longer  time  record  would  be  desirable  for  a  few  experiments  as  indi¬ 
cated  by  Howe  et  al.  (1978).  However,  the  10.24-minute  record  length 
was  chosen  as  a  compromise  which  allowed  a  vertical  profile  to  be  com¬ 
pleted  in  a  reasonable  time,  thereby  minimizing  the  effects  of  changes 
in  ambient  and  bulk-water  conditions. 

The  channel  was  filled  with  tap  water  at  least  three  days  before 
each  set  of  experiments  to  allow  time  for  de-aeration.  Chemical  treat¬ 
ment  of  the  water  to  prevent  algae  growth  was  unnecessary.  Little  or  no 
seeding  of  the  water  was  needed  for  the  LDA  in  fixed-frame  experiments 
where  forward-scatter  mode  was  used.  However,  it  was  essential  to  seed 
the  water  with  titanium  dioxide  to  obtain  sufficient  numbers  of  good 
Doppler  bursts  in  wave- following  runs  where  the  backward-scatter  mode 
was  used. 

There  are  17  experiments  which  comprise  the  three  main  cases  in 
this  study.  The  general  experimental  conditions  for  each  experiment  are 
tabulated  in  Table  1.  The  typical  sequence  of  steps  taken  for  each 
experiment  is  as  follows: 

1.  The  LDA  probe  volume  was  adjusted  to  coincide  with  the  initial 
still  water  level;  this  was  the  Initial  depth  reference  level. 
Manometer  levels  were  also  recorded  under  the  no-wind  condition. 

2.  The  fan  was  turned  on  and  the  desired  wind  speed  set.  After  the 
wind  speed  reached  a  steady  state,  the  manometer  levels  were 
recorded  to  note  the  water  set-up.  The  fan  was  left  on  overnight 
(about  16  hr  between  successive  experiments)  to  assure  that  the  air 
and  water  flow  were  in  equilibrium. 


3. *  The  wave  plate  was  turned  on  at  least  an  hour  before  the  first  data 

run  to  ensure  full  development  of  the  water  and  air  wave  fields. 
The  wave  generator  was  set  to  produce  a  wave  amplitude  of  about 
22  mm  and  a  frequency  of  1.00  Hz. 

4.  The  manometer  levels  were  recorded  just  prior  to  the  first  run  to 
find  the  initial  water  loss. 

5.  The  LDA  table  was  lowered  for  the  first  data  point.  The  LDA  probe 
volume  was  brought  as  close  to  the  interface  as  possible,  without 
letting  it  leave  the  water. 

6.  The  data  acquisition  program  was  activated  and  all  necessary  in 
line  calibrations  were  made,  such  as  for  the  trackers  and  the  pres¬ 
sure  transducer. 

7. ^  The  scanning  amplitude  was  adjusted  by  varying  the  calibrated 

attenuator  of  the  scanner  control  unit  so  that  the  LDA  probe  volume 
followed  the  transformed  coordinate  system  (y*  ■  constant). 

8.  The  computer  was  instructed  to  commence  the  data  taking  routine. 

9.  The  LDA  table  was  lowered  for  the  next  point  and  Steps  7  through  9 
were  repeated  until  sufficient  number  of  data  points  (at  about  20 
depths)  were  obtained. 

10.  Once  again,  the  manometer  levels  were  recorded  to  estimate  the 
final  water  loss. 

11.  The  fan  was  turned  off  and  when  quiescence  was  reached,  the  manome¬ 
ter  levels  were  recorded  to  give  a  second  estimate  of  the  total 
water  loss. 

*For  mechanically  generated  wave  cases  only. 

^For  Eulerlan  wave-following  case  only. 


12.  Water  was  added  to  the  channel  if  necessary  and  the  whole  procedure 
was  repeated  for  the  next  wind  speed. 

3.5  DATA  REDUCTION  PROCEDURES 

In  this  section,  the  data  reduction  procedures  for  each  of  the  three 
sets  of  experiments  in  this  study  are  discussed.  Then  in  Sec.  3. 5. A,  a  com¬ 
parison  between  the  wave-separation  and  phase-averaging  methods  is  given. 

3.5.1  Wind-Generated  Wave  Case 

The  wave-separation  method  (Sec.  2.2.2)  was  used  to  deduce  the 
means  and  second  moments  of  all  quantities  of  interest.  Three  major 
steps  were  necessary  in  analyzing  the  data,  viz., 

1 .  The  time  mean  for  each  channel  of  data  in  an  experiment 
was  calculated  and  the  results  stored  in  a  disc  file. 

2.  The  wave-separation  program  was  run  to  find  the  wave- 
induced  and  turbulent  quantities.  The  results  were 
stored  in  a  disc  file. 

3.  The  mean,  wave,  and  turbulence  data  files  were  combined 
to  form  a  summary  file  for  an  experiment. 

The  autospectra  and  cross-spectra  among  the  velocity  components 
(u,  v)  and  the  water  surface  elevation  (n)  were  tabulated  as  well  as 
plotted  in  Step  2.  However,  the  spectra  were  not  stored  on  disc  due  to 
the  large  amount  of  disc  space  needed. 

3.5.2  Mechanical-Wave  Eulerian  Frame  Case 

As  indicated  in  Sec.  2.5,  it  is  possible  to  analyze  this  case  using 
the  wave-separation  method  and  the  phase-averaging  scheme.  The  wave- 
separation  method  used  was  identical  to  that  of  the  wind-wave  case.  The 
steps  involved  for  the  phase-averaging  data  analysis  procedure  are  as 
follows : 
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1.  The  time  mean  for  each  channel  of  data  was  calculated. 

2.  The  wave  period  of  the  reference  wave-height  signal  was 
found  by  counting  the  number  of  sampled  points  between 
two  successive  zero  crossings  with  positive  slope. 

3.  The  phase  averages  were  calculated  as  defined  by  Eq.  2-12. 

4.  The  mean  turbulent  stresses  (u^u^),  mean  wave-induced 
stresses  (u^uj),  and  other  wave  related  quantities  (ujuj) 
were  calculated. 

5.  The  amplitude  and  phase  of  u,  v,  u'u',  and  n  were  found 
by  the  cross-spectral  method  and  by  the  cosine-fitting 
techniques  (see  Appendices  A  and  C) . 

6.  The  mean,  phase-averaged,  and  turbulence  data  were  stored 
in  a  disc  data  file. 

3.5.3  Mechanical-Wave  Wave-Following  Euler ian  Frame  Case 


The  time-  and  phase-averaging  data  reduction  procedures  in 
Sec.  3.5.2  were  also  used  for  this  case.  However,  the  vertical  velocity 
measured  has  to  be  corrected  because  of  the  motion  of  the  wave-follower. 
The  time  derivative  of  the  wave-follower  position  output  voltage  is 
related  to  the  vertical  velocity  component  measured  by  the  LDA.  A 
central  differencing  scheme  was  used  to  differentiate  this  voltage.  The 
result  was  converted  into  a  velocity  by  use  of  the  constants  given  in 
Sec.  3.3.7.  Therefore,  the  actual  vertical  velocity  is  given  by 

v  "  VM  +  *WF 

where  vM  and  v^F  are  the  measured  and  the  wave-follower  scan  velocities, 
respectively.  Because  the  driving  signal  of  the  wave-follower  was  band¬ 
pass  filtered  and  periodic  with  zero  mean,  v'  ■  v  ■  0. 


Ocher  quantities  of  interest  were  found  using  the  same  procedures 
as  used  for  the  fixed-frame  case  (Sec.  3.5.2). 

3.5.4  A  Quantitative  Comparison  Between  the  Wave-Separation  and  Phase- 
Averaging  Methods 

As  noted  earlier  in  Sec.  2.5,  it  is  possible  to  use  both  the  wave- 
separation  and  phase-averaging  methods  in  analyzing  mechanical-wave 
fixed-frame  measurements.  Because  for  no-wind  and  low-wind-speed  cases, 
n'  <<  n,  the  results  obtained  by  the  wave-separation  method  should  be 
good  approximations  to  those  given  via  the  phase-averaging  process. 

Two  experiments  were  run  with  the  following  objectives: 

1.  to  show  that  phase  averaging  is  a  better  method  in  ana¬ 
lyzing  mechanical  wave  data,  and 

2.  to  infer  from  (1)  the  errors  made  by  the  wave-separation 
method  in  reducing  wind -wave  data. 

The  experiments  were  performed  at  zero  wind  speed  and  under  mechan¬ 
ically  generated  waves.  A  short  vertical  velocity  profile  was  taken  in 
each  experiment.  The  sampling  rate  was  200  Hz  for  the  first  experiment 
and  100  Hz  (i.e.  the  same  as  that  used  in  wind-wave  experiments)  for  the 
second  one.  Because  there  was  no  wind  and  the  waves  were  essentially 
irrotational ,  except  for  the  very  thin  Stokes  boundary  layers  along  the 
channel  sides  and  at  the  interface  due  to  viscous  effects,  all  the 
turbulent  quantities  should  be  very  small  for  both  test  cases. 

The  results  of  the  experiments  are  presented  in  Table  6.  For  the 
first  experiment,  it  is  clear  that  the  turbulent  quantities  obtained  by 
the  phase-averaging  method  are  more  credible  than  those  given  by  the 
wave-separation  scheme.  The  turbulent  quantities  are  very  small,  which 
is  representative  for  the  test  flow,  and  are  within  the  uncertainty 
intervals  of  the  LDA.  The  wave-separation  method  gives  a  much  higher 


turbulence  level;  however,  the  separation-method  wave  quantities  are  in 

excellent  agreement  with  the  phase-averaging  results. 

The  high  turbulence  level  given  by  the  wave-separation  method  can 
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be  explained  as  follows.  The  mean  square  wave  quantities,  such  as  u  , 

are  computed  by  integrating  the  spectra  (Appendix  A).  Next,  the  mean 

square  turbulent  component  is  obtained  by  subtracting  the  mean  square 

— 2  ~7  77 

wave  component  from  the  total  mean  square  value,  i.e.  u*  *  u  -  u  . 

Because  the  total  mean  square  value  is  the  sum  of  the  wave  and  turbulent 

7? 

parts,  errors  introduced  in  computing  u  ,  say  by  numerical  integration, 

— I 

lead  to  errors  in  u'  .  If  we  assume  the  "true"  results  are  those  given 

by  the  phase-averaging  method,  for  example  the  first  data  point  in 

7 2  -7  — jr  -y 

Table  6,  we  have  for  phase  averaging:  u  /u  =  0.99,  u'  / u  *  0.01  and 

77  ~2  — 2 

for  wave  separation:  u  /u  a  0.96,  u*  /u  *  0.04.  This  means  that  the 

77  — 2 

wave-separation  method  has  a  4%  error  in  computing  u  ,  hence  in  u' 
also.  Although  this  is  a  very  small  percentage  error  in  terms  of  mean 
square  values,  the  percentage  error  in  root  mean  square  values  is  much 

— j  1/2  — ^  1/2 

larger:  (u'  )  /(u  )  *  20% !  This  is  the  reason  for  the  seemingly 

high  turbulence  level  in  the  wave-separation  scheme.  So,  based  on  the 
results  of  this  test,  the  phase-averaging  process  is  shown  to  be  the 
better  method. 

As  pointed  out  earlier,  the  phase-averaging  scheme  cannot  be  used 
for  wind-wave  experiments.  However,  as  shown  above,  the  wave-separation 
method  can  introduce  seemingly  large  errors  in  the  calculation  of  turbu¬ 
lent  quantities.  Therefore,  it  is  necessary  to  establish  some  error 
bounds  on  the  data  obtained  by  this  method. 

A  second  experiment  was  run  and  analyzed  with  all  the  sampling 
parameters  identical  to  those  for  the  wind-wave  experiments.  The 
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following  results  are  obtained,  1)  phase-averaging:  max(u  /u  )  * 

■  ■  j  — * 

0.99,  max(u'  /u  )  »  0.01  and  11)  wave  separation:  max(u  /u  )  *  0.98, 

— I  T 

max(u'  /u  )  «  0.02.  So  the  error  Is  reduced  for  the  wave-separation 
method;  furthermore,  the  percentage  error  In  terms  of  the  root  mean 
square  value  Is  at  most  13%.  The  improvement  is  due  mainly  to  a  lower 
sampling  rate  that  results  in  a  narrower  signal  bandwidth.  Also  note 
that  the  turbulent  quantities  given  by  the  phase-averaging  method  are 
slightly  higher  than  those  in  the  first  experiment.  This  is  due  to  a 
coarser  resolution  of  the  main  wave  when  it  is  sampled  at  a  lower  fre¬ 
quency.  The  result  is  a  larger  phase  jitter  which  leads  to  a  higher 
turbulence  level.  Although  the  test  was  for  mechanically  generated 
waves,  we  believe  that  the  percentage  error  in  root  mean  square  values 
(about  13%)  is  also  applicable  for  wind-wave  experiments.  This  error 
bound  is  assumed  and  used  in  subsequent  sections  for  wind-wave  data. 

3.6  UNCERTAINTY  ESTIMATES 

The  N^-order  uncertainty  intervals  as  defined  by  Moffat  (1981) 
were  estimated  for  all  important  variables.  The  Nth-order  uncertainty 
of  a  measured  quantity  is  an  estimate  of  its  absolute  accuracy  and 
includes  Instrument  calibration  uncertainty  as  well  as  the  effects  of 
unsteadiness  and  interpolation.  The  use  of  such  uncertainty  estimates 
is  essential  for  reporting  results  and  assessing  the  significance  of 
differences  among  results  from  different  experiments. 

Table  7  contains  the  estimated  uncertainties.  They  were  calculated 
by  using  the  procedure  described  by  Kline  and  McClintock  (1953).  As 
indicated,  many  of  the  estimates  are  dependent  upon  actual  flow  condi- 
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tiona. 


Table  2 


Simultaneous  Forward-Scatter  Backward-Scatter  Runs 
(Tracker  Range  333  kHz,  *  0.0  m/s) 


UB  (BS) 

uB  (FS) 

uB  (BS) 

uB  (FS) 

(BS) 

/  Ug2  (FS) 

/  Ug2  (CC) 

(mm/s) 

(mm/s) 

(mm/s) 

(mm/s) 

(mm/s) 

(mm/s) 

(om/s) 

92.0 

93.3 

-1.31 

-2.27 

5.78 

3.63 

3.12 

83.2 

84.9 

-1.62 

-2.85 

6.02 

2.60 

1.73 

75.3 

76.5 

-2.05 

-3.20 

5.16 

2.72 

1.91 

68.2 

69.1 

-1.90 

-2.79 

3.92 

2.52 

1.68 

56.4 

57.7 

-0.89 

-2.24 

5.07 

2.24 

1.11 

50.5 

51.3 

-0.75 

-1.80 

4.04 

2.32 

1.31 

43.7 

44.6 

-0.38 

-1.46 

4.08 

2.08 

0.81 

35.3 

36.0 

0.69 

-1.12 

3.72 

2.11 

0.86 

Notation:  BS,  backward-scatter  mode;  FS,  forward-scatter  mode;  CC, 

cross-correlated.  The  subscript  “B"  denotes  the  LDA  blue  channel. 


Table  3 

Simultaneous  Forward-Scatter  Backward-Scatter  Runs 
(Tracker  Range  333  kHz,  u.  »  4.1  m/s) 


I 

uB  (BS) 
(mm/s) 

uB  (FS) 
(mm/s) 

uB  (BS) 
(mm/s) 

uB  (FS) 
(mm/s) 

/  ^  (BS)  » 
(mm/s) 

'  u'|  (FS) 
(mm/s) 

'  Ug2  (CC) 
(mm/s ) 

89.6 

90.6 

3.23 

2.70 

6.30 

5.57 

5.21 

81.1 

81.9 

1.86 

1.02 

5.86 

4.75 

4.37 

* 

74.5 

75.1 

-0.48 

-1.29 

5.16 

4.14 

3.72 

,  •  .. 

67.1 

68.3 

1.52 

-1.02 

6.39 

4.08 

3.58 

f-  . 

r-  - 

57.4 

59.8 

0.66 

-1.46 

8.38 

3.74 

3.21 

k 

52.1 

52.9 

-0.68 

-1.65 

4.83 

3.29 

2.68 

9 

r.  -  ■ 

44.5 

45.5 

0.02 

-1.15 

5.49 

3.18 

3.17 

'v 

c.  • 

36.6 

37.1 

-0.69 

» 

H-* 

• 

wn 

o 

4.50 

3.01 

2.40 

Notation:  BS,  backward-scatter  mode;  FS,  forward-scatter  mode;  CC, 

cross-correlated.  The  subscript  "B"  denotes  the  LDA  blue  channel. 
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Table  6 


Comparison  of  Wave-Separation  (Method  1) 
and  Phase-Averaging  (Method  2)  Methods 


(nun/s) 

FI 

/  U 

Me  thod  1 

Method  2 

(mm/s) 

/7 

Method  1 

Method  2 

/? 

f? 

/=f 

i  U 

/7 

/? 

/T* 

ft 

r? 

Experiment  1  (f 

sample 

-  200 

Hz) 

73.3 

72.0 

13.4 

73.3 

1.7 

77.8 

76.2 

15.7 

77.7 

1.7 

68.2 

67.2 

11.6 

68.0 

2.6 

72.3 

70.8 

14.6 

72.3 

2.6 

64.3 

63.4 

10.8 

64.5 

1.4 

68.1 

66.7 

13.6 

67.9 

1.6 

60.0 

59.2 

9.9 

59.8 

1.8 

63.5 

62.3 

12.7 

63.6 

1.9 

56.4 

55.6 

9.1 

56.4 

1.4 

59.8 

58.6 

11.8 

59.7 

1.5 

52.2 

51.6 

8.2 

52.3 

1.6 

55.3 

54.2 

10.9 

55.1 

1.6 

47.9 

47.4 

7.4 

47.9 

1.4 

50.7 

49.7 

10.0 

50.7 

1.6 

44.0 

43.5 

6.6 

44.0 

1.3 

46.4 

45.5 

9.2 

46.4 

1.4 

40.1 

39.6 

5.9 

40.0 

1.6 

42.3 

41.5 

8.3 

42.4 

1.8 

35.6 

35.3 

5.2 

35.5 

1.0 

37.5 

36.7 

7.3 

37.5 

1.1 

32.1 

31.8 

4.6 

32.0 

1.0 

33.6 

33.0 

6.6 

33.6 

1.1 

29.0 

28.7 

4.2 

29.0 

1.0 

30.2 

29.6 

5.9 

30.2 

1.1 

26.3 

26.0 

3.8 

26.2 

1.2 

27.4 

26.9 

5.4 

27.4 

1.2 

23.4 

23.1 

3.4 

23.3 

1.1 

24.2 

23.8 

4.8 

24.2 

1.1 

Experiment  2  (f 

sample 

-  100 

Hz) 

73.7 

73.0 

9.9 

73.5 

3.4 

77.7 

77.2 

8.9 

77.8 

3.3 

69.3 

68.8 

8.1 

69.2 

6.3 

73.8 

73.3 

8.7 

73.5 

6.6 

66.5 

66.0 

7.4 

66.6 

2.9 

70.3 

69.8 

8.4 

70.0 

3.1 

63.2 

62.8 

6.9 

63.2 

2.9 

65.0 

64.5 

7.9 

64.8 

3.3 

59.8 

59.4 

6.6 

59.5 

3.2 

61.0 

60.6 

7.5 

61.1 

3.5 

55.6 

55.3 

5.9 

55.4 

2.5 

58.6 

58.2 

7.2 

58.7 

2.6 

52.2 

51.9 

5.3 

52.3 

3.0 

55.2 

54.8 

6.9 

54.9 

3.3 

48.4 

48.4 

4.6 

48.3 

2.5 

51.3 

51.0 

6.1 

51.2 

2.8 

47.6 

47.3 

5.6 

47.4 

3.4 

54.0 

53.5 

7.1 

54.0 

3.9 

35.5 

35.4 

3.1 

35.4 

2.0 

37.0 

36.7 

4.8 

37.0 

2.1 

33.3 

33.2 

2.7 

33.2 

1.4 

34.7 

34.4 

4.4 

34.7 

1.6 

30.7 

30.6 

2.5 

30.7 

1.1 

32.1 

31.8 

4.1 

32.1 

1.2 

28.3 

28.2 

2.3 

28.3 

1.6 

29.4 

29.2 

3.7 

29.4 

1.7 

26.0 

25.9 

2.2 

26.1 

1.0 

27.1 

36.8 

3.4 

26.9 

1.2 

24.2 

24.1 

2.1 

24.2 

2.0 

25.0 

24.8 

3.2 

24.8 

2.2 

22.3 

22.2 

1.9 

22.1 

1.1 

23.0 

22.8 

3.0 

23.1 

1.3 

20.7 

20.6 

1.7 

20.6 

0.9 

21.3 

21.1 

3.0 

21.3 

1.3 

19.0 

19.0 

1.7 

19.1 

1.0 

19.6 

19.4 

2.6 

19.4 

1.1 

17.9 

17.8 

2.0 

17.8 

1.5 

18.3 

18.2 

2.5 

18.4 

1.4 

- : 
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Table  7 


Order  Uncertainty  Estimates  at  20:1  Odds 


Quantity 


Typical  Uncertainty 


1,  V 

2%  or  ±  0.6  mm/s  100  kHz  Tracking  Range 

±2.0  mm/s  333  kHz  Tracking  Range 

— j  1/2 

u’  )  , 

— y  1/2 

(V  ) 

13%  of  total  RMS  value 

— y  1/2 

u  )  , 

— y  1/2 

<v  ) 

2%  of  total  RMS  value 

■u'v' 

20%  in  fixed-frame  experiments 

30%  in  wave-following  frame  experiments 

■u  V 

10%  in  fixed-frame  experiments 

20%  in  wave-following  frame  experiments 

.  A 

».  V 

2% 

'u*  0$ 

1° 

20% 

>» 

20° 

i»  n 

2%  or  ±  0.3  mm 

0.5%  or  ±  0.1°C 

^Bulk 

0.5%  or  ±  0.1°C 

»« 

2-6%,  3%  at  6  m/s 

r 

±  0.5  mm 

CHAPTER  4 


EXPERIMENTAL  RESULTS 

As  stated  In  Chap.  2,  the  present  study  consists  of  three  sets  of 
experiments: 

1.  wind-wave  in  Eulerian  frame; 

2.  wind-ruffled  mechanically  generated  wave  in  Eulerian 
frame;  and 

3.  wind-ruffled  mechanically  generated  wave  in  wave¬ 
following  Eulerian  frame. 

The  results  for  all  three  cases  are  given  in  the  following  sections. 
These  results  include  the  mean,  wave-induced,  and  turbulent  velocity 
quantities  presented  in  terms  of  various  dimensional  and  non-dimensional 
parameters.  The  general  trends  and  special  behaviors  that  follow  natu¬ 
rally  from  the  data  are  indicated  as  they  occur.  However,  a  discussion 
of  the  implications  of  these  results  is  withheld  and  presented  in 
Chap.  5  after  all  the  results  are  shown  here. 

4.1  WIND  WAVE  IN  EULERIAN  FRAME 

4.1.1  Water  Surface  Conditions 

The  water  surface  varied  from  smooth  to  rough  in  the  literal  sense 
for  the  range  of  reference  wind  speeds  used  in  this  set  of  experiments. 
Time  traces  of  typical  water  surface  conditions  are  shown  in  Fig.  8.  Each 
trace  is  four  seconds  in  duration.  These  time  traces  represent  the  water 
surface  elevation  at  a  point  as  the  waves  pass  by  the  wave-height  gauge. 

For  the  lowest  wind-speed  case  (1.5  m/s),  no  visible  gravity  wave 
was  present;  however,  small  capillary-type  perturbations  were  observed 
during  the  runs.  The  wave-height  gauge  is  unable  to  detect  such  small 


disturbances  (<  0.2  nm).  Tiny  wavelets  were  present  for  the  2.6  m/s 
case,  but  the  water  surface  remained  relatively  smooth.  At  higher  wind 
speeds,  the  waves  were  clearly  visible.  The  wave  amplitude  and  the  wave 
slope  increase  with  increasing  wind  speed.  Small  wave  breaking  was 
observed  when  the  wind  speed  was  greater  than  6.7  m/s  as  indicated  by 
the  pointedness  of  the  wave  crests  and  small  discontinuities  in  the 
wave-height-gauge  signal. 

Another  parameter  which  is  often  used  in  describing  the  surface 
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conditions  is  the  root-mean-square  (RMS)  wave  amplitude  (rr)  •  This 
quantity  is  plotted  versus  wind  speed  in  Fig.  9.  The  RMS  n  is  less  than 
0.2  mm  for  the  two  lowest  wind  speeds,  but  increases  almost  linearly  at 
higher  wind  speeds. 

The  dominant  frequency  of  the  wind-generated  waves  decreases  with 
increasing  wind  speed.  This  gradual  shift  of  dominant  wave  frequency 
with  wind  speed  can  be  seen  from  the  time  traces.  All  the  measured 
wave-height  power  spectra  (not  shown  here)  showed  a  local  maximum  at  a 
frequency  close  to  twice  the  peak  frequency.  These  results  are  consis¬ 
tent  with  the  findings  of  Sutherland  (1968),  Plate  et  al.  (1969),  and 
Chambers  et  al .  (1970). 

4.1.2  Surface  Drift  Velocity 

The  surface  drift  velocity  measured  by  timing  the  paper  punchings 
is  a  Lagrangian  velocity.  It  is  necessary  to  convert  this  Lagrangian 
velocity  to  a  Eulerian  velocity  in  order  to  use  it  with  the  LDA  measure¬ 
ments,  which  are  Eulerian.  A  relation  between  the  Lagrangian  velocity 
UL(x,  t)  and  the  Eulerian  velocity  u(x,  t)  induced  by  water  waves, 
correct  to  the  second  order  of  invlscid  water-wave  theory,  is  given  by 
(see  Phillips,  1977): 
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u^(_x ,  t)  “  u(£,  t)  +  (/  u(x,  t' )dt’ )*Vu(x,  t) 


-  “<*.  o  +  “St0ke,(*.  O 

The  second  term  on  the  right-hand  side  of  the  above  equation  is  equiva¬ 
lent  to  the  Stokes  mass  transport  velocity.  Assuming  a  linear  theory, 
this  term  can  be  approximated  from  a  one-dimensional  surface  elevation 
spectrum  (Srif)(f))  and  is  (see  Bye,  1967): 


(Wf 

Stokes  ‘  g  nn 
o 


(4-1) 


where  f  taken  to  be  10  Hz  here,  is  the  upper  frequency  limit  of  grav¬ 
ity  waves.  Therefore,  the  mean  Eulerian  surface  drift  velocity  (u,.)  at 
the  location  of  the  wave-height  gauge  is: 


US  “  USL  "  uStokes 

where  UgL  is  the  mean  Lagrangian  surface  drift  velocity. 

The  variation  in  the  mean  Eulerian  and  Lagrangian  surface  drift 
velocities  with  wind  speed  is  shown  in  Fig.  10.  The  Eulerian  surface 
drift  velocity  is  found  to  be  about  2%  of  the  free  stream  velocity.  If 
no  correction  to  the  Lagrangian  velocity  is  made,  then  the  Lagrangian 
surface  drift  velocity  is  about  3%  of  the  free-stream  wind  speed.  These 
results  are  consistent  with  the  findings  of  other  investigators  (Wu, 
1968;  Phillips  and  Banner,  1974). 

4.1.3  Mean  Velocity  Profiles 

Mean  velocity  profiles  are  shown  in  Fig  11.  The  velocity  defect  ug  - 
u  is  plotted  instead  of  the  mean  velocity  u  so  that  the  profiles  may 
resemble  those  for  flows  over  a  solid  surface.  This  allows  comparison 
with  other  well-established  data  sets,  such  as  those  for  flows  over 
smooth  or  rough  plates.  Full  profiles  were  not  taken  because  the  main 


51 


goal  of  this  study  is  to  investigate  the  surface  and  near-surface  flow 
phenomena.  The  maximum  depth  for  velocity  measurements  is  in  the  vicin¬ 
ity  of  300  mm  or  roughly  one  third  of  the  total  water  depth. 

It  is  evident  that  presence  of  the  moving  waveforms  makes  it  diffi¬ 
cult  to  approach  close  enough  to  the  surface  with  the  LDA  to  observe  the 
behavior  of  mean  velocity  profiles  near  the  interface.  For  high  wind- 
speed  cases,  the  measuring  point  closest  to  the  surface  is  of  the  order 
of  10  to  20  mm  from  the  mean  water  level.  The  general  trend  of  the 
profiles  is  similar  to  flows  over  solid  surfaces.  However,  there  are 
differences  especially  when  the  profiles  are  plotted  in  law-of-the-wall 
coordinates. 

In  Fig.  12,  mean  velocity  profiles  are  shown  in  wall  coordinates, 
where  u+  *  (u<.  -  u)/u*  and  y+  ■  yu*/vw-  The  flow  was  turbulent, 
as  defined  by  Tennekes  and  Lumley  (1972),  even  at  the  lowest  wind 
speed.  The  friction  velocity  u*  is  determined  by  a  linear  least-squares 
fit  (with  depth)  of  the  direct  measurements  of  the  turbulent  shear 
stress  (-uTvT)  and  extrapolation  of  the  trend  line  to  the  interface.  A 
detailed  account  of  the  turbulent  shear  stresses  is  given  in  Sec.  4.1.4. 
Also  plotted  in  the  figure  are  the  universal  law  of  the  wall  (Schlichting, 
1979)  for  a  turbulent  flow  with  zero  pressure  gradient  over  smooth  walls: 

u+  ■  i  ln(-y+)  +  5.5  (4-2) 
and  that  for  the  rough  walls: 

u+  -  I  ln(-y+)  -  2.1  (4-3) 

where  k  (■  0.4)  is  the  von  Karman  constant.  In  wall  coordinates,  the 
velocity  profile  for  the  lowest  wind  speed  clusters  around  the  smooth 
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wall  relationship  and  as  wind  speed  increases,  the  usual  shift  associ¬ 


ated  with  the  apparent  increase  In  surface  roughness  occurs.  For  wind 
speeds  higher  than  3.2  m/s,  all  the  profiles  lie  below  the  relationship 
representing  the  start  of  the  fully  rough  flow  regime.  All  the  profiles 
exhibit  logarithmic  behavior.  The  profiles  for  wind  speeds  less  than 
3.3  m/s  have  slopes  close  to  2.5  (*  I/O*  However,  for  higher  wind 
speeds,  the  profiles  have  a  slope  smaller  than  2.5.  Such  seemingly 
anomalous  behavior  is  also  reported  by  Howe  et  al .  (1982). 

The  mean  vertical  velocity  profiles  are  shown  in  Fig.  13.  The 
scales  of  this  plot  are  the  same  as  those  used  for  the  horizontal  velo¬ 
city  defect  profiles  for  easy  comparison.  The  actual  data  are  tabulated 
in  Appendix  D.  For  all  wind  speeds,  with  the  exception  of  the  highest, 
the  mean  vertical  velocity  profiles  fall  within  a  band  of  i  3  mm/s  about 
zero.  At  the  highest  wind  speed,  the  mean  vertical  velocity  is  about 
zero  near  the  surface  and  increases  with  depth.  However,  the  mean 

vertical  velocity  is  very  much  smaller  than  the  mean  horizontal  velocity 
defect  for  all  runs.  Furthermore,  the  mean  vertical  velocity  gradient 
(with  depth)  is  very  close  to  zero  near  and  away  from  the  interface. 
4.1.4  Mean  Turbulent  Quantities 

The  turbulent  shear  stress  (-u'v')  profiles  for  all  the  wind-wave 
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experiments  are  shown  in  non-dimensional  forms  in  Fig.  14.  The  u* 
values  are  obtained  by  a  linear  least-squares  fit  with  depth  of  the 
direct  measurements  of  the  shear  stress  and  extrapolation  of  the  trend 
line  to  the  interface.  The  depth  is  normalized  by  the  depth  (Delta(S)] 
at  which  the  turbulent  shear  stress  vanishes.  The  Deltas  are  determined 
from  the  same  leas  t-squares-f  it  lines  used  in  deducing  u*;  that  is, 
Delta  is  the  depth  at  which  the  line  meets  the  y-axis.  This  definition 


of  Delta  has  the  equivalent  physical  meaning  of  a  boundary- layer 
thickness  for  the  turbulent  boundary  layer  (Goossens  et  al.,  1982). 

Relatively  large  amounts  of  data  scatter  are  present  for  low  speed 

runs  or  at  depths  where  the  turbulent  shear  stress  is  expected  to  be 

very  small.  This  Is  due  to  the  small  magnitude  of  the  shear  stress  (of 

the  order  of  10  °  nr/s  )  which  approaches  the  resolution  limit  of  the 
—ft  9  9 

LDA  system  (*  10  °  nr/s'y.  However,  data  scatter  Is  very  much  reduced 
for  higher  wind-speed  runs  where  the  linear  behavior  of  the  turbulent 
shear  stress  profiles  is  obvious. 

A  decrease  of  the  shear  stress  near  the  surface  can  be  seen  for  the 
two  low  speed  runs  as  well  as  some  higher  wind-speed  flows.  However, 
the  measurements  were  made  in  a  fixed  frame.  Thus,  the  shear  stress 
values  obtained  near  the  interface  were  averaged  over  points  which  were 
very  close  to  the  surface  (near  the  wave  trough)  and  points  which  were 
some  distance  away  from  the  instantaneous  water  surface  (under  the  wave 
peak).  Accordingly,  it  is  unwise  to  conclude  that  such  a  decrease  in 
the  turbulent  shear  stress  is  physical  and  not  due  to  kinematics  of  the 
waves . 

On  the  other  hand,  the  two  low  wind-speed  experiments  are  fairly 
close  to  the  smooth  flow  regime  as  indicated  by  the  mean  profiles  in 
wall  coordinates.  Then,  the  decrease  of  turbulent  shear  stress  near  the 
surface  is  due  to  the  existence  of  a  viscous  sublayer.  Another  piece  of 
evidence  for  the  existence  of  the  viscous  sublayer  is  that  the  velocity 
profiles  adhere  close  to  the  u+  ■  -y+  line  for  these  two  low  wind- 
speed  cases. 

The  friction  velocities  (u*)  are  plotted  versus  wind  speed  (u*,)  in 
Fig.  15.  The  friction  velocity  grows  almost  linearly  with  wind  speed 


for  u.  <  5  m/s  and  more  rapidly  at  higher  wind  speeds.  This  result 
agrees  with  that  of  Howe  et  al.  (1981,  1982).  The  corresponding  turbu¬ 
lent  boundary-layer  thickness  6  can  be  found  in  Appendix  D.  Delta  is 
about  140  mm  for  the  lowest  wind  speed  and  stays  relatively  constant 
around  300  mm  for  higher  wind  speeds  where  waves  are  present. 

— *-  1/2  — k-  1/2 

The  root-mean-square  (RMS)  velocity  fluctuation  [(u'  )  ,  (v'  ) 

profiles  are  plotted  in  non-dimensional  forms  in  Figs.  16  and  17.  The 
RMS  u'  profiles  for  all  wind  speeds  collapse  within  a  narrow  band  when 
normalized  by  u*  and  6.  The  RMS  v'  profiles  for  wind  speeds  higher  than 
3  m/s  also  coalesce  into  a  band  when  normalized  by  the  same  parameters. 
Clearly,  u*  and  6  are  the  appropriate  velocity  and  length  scales  relat¬ 
ing  turbulent  velocity  fluctuations  and  turbulence  decay  with  depth. 

Also  plotted  in  the  same  figures  are  the  data  of  Klebanoff  (1955) 
for  a  turbulent  boundary  layer  with  zero  pressure  gradient.  The  RMS  u' 
profiles  follow  this  reference  data  reasonably  well,  but  generally  lie 
above  it.  The  RMS  v'  profiles  also  group  around  the  reference  data 
except  for  the  lower  wind-speed  experiments.  It  is  clear  that  there  is 
anisotropy  between  RMS  u'  and  v'  within  the  boundary  layer.  In  addi¬ 
tion,  the  degree  of  anisotropy  increases  towards  the  interface.  How¬ 
ever,  it  is  unclear  whether  this  trend  continues  into  the  inner  part  of 
the  boundary  layer  due  to  extreme  difficulty  in  obtaining  data  above  the 
wave  trough  in  a  fixed  frame.  On  the  other  hand,  turbulent  fluctuations 
decrease  and  become  more  isotropic  near  the  edge  of  the  boundary  layer. 

The  RMS  v'  in  normalized  units  for  wind  speeds  between  1.5-3. 2  m/s 
does  not  seem  to  fall  into  the  general  trend  of  the  data.  In  this  wind 
speed  range,  the  non-dimensional  RMS  v'  increases  with  increasing  wind 
speed,  but  lies  below  the  main  sequence  of  the  rest  of  the  data. 
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A. 1.5  Mean  Wave-Induced  Quantities 

Figure  18  shows  the  mean  wave-induced  shear  stress  (-u  v)  pro¬ 
files*  The  normalization  parameters  are  the  same  as  those  used  for  the 
turbulent  shear  stress  profiles,  thus  allowing  direct  comparison  between 
the  two  components  of  the  total  shear  stress.  The  wave-induced  shear 
stress  is  negative  very  close  to  the  interface  and  decays  very  rapidly 
with  depth.  One  can  define  a  penetration  depth  for  the  wave-induced 
shear  stress  as  the  depth  at  which  the  wave-induced  shear  stress  van¬ 
ishes.  Then,  the  penetration  depth  of  the  wave-induced  shear  stress 
increases  with  wind  speed  but  is  at  most  equal  to  0.26  for  the  highest 
wind  speed  and  is  considerably  less  for  lower  wind-speed  experiments. 
The  wave-induced  shear  stress  is  of  the  order  of  u£  near  the  surface  and 
is  a  larger  percentage  of  the  total  shear  stress  within  that  region. 
However,  the  turbulent  shear  stress  seems  to  be  the  dominant  component 
over  most  parts  of  the  boundary  layer. 

TJf  1/2  ^  1/2 

The  RMS  wave-induced  velocity  profiles  [(u  )  ,  (v  )  ]  are  shown 

in  Figs.  19  and  20.  The  normalizing  parameters  are  the  RMS  surface 
TJ  1/2 

orbital  velocity  (uq)  computed  from  the  wave-height  signal  and  the 

TJf 

wave-decay  depth  y^.  First,  uq  is  given  by: 


T?  ,c  2 

u‘  -  /  (2*f)\  (f)df 

o 


(4-4) 


where  S^^f)  is  the  one-sided  power  spectrum  of  the  wave-height  signal 
and  fc  (■  10  Hz)  is  the  upper  frequency  limit  for  gravity  waves.  Next, 
yn  is  the  depth  at  which  the  integrated-surface  orbital  amplitude  decays 
by  a  factor  of  e~l  (Howe  et  al.,  1982);  hence,  y  satisfies  the  follow¬ 


ing  relationship: 


(4-5) 


.fc  2 
/  C  S  (f)  exp[-2( 2uf )  y  / g]df 

o  n  -2 

f  "  e 

/  C  S  (£ )df 
_  nn 
o 

Here,  y^  Is  used  instead  of  some  other  length  scale,  such  as  the  domi¬ 
nant  wavelength,  because  y^  takes  into  account  the  amplitude  and  fre¬ 
quency  distribution  of  the  entire  wind-generated  surface  spectrum.  In 
fact,  the  definition  of  y^  is  equivalent  to  the  reciprocal  of  an  inte¬ 
grated  wave  number. 

The  data  collapse  is  excellent  using  these  parameters.  Some  of  the 

data  scatter  is  partly  due  to  the  logarithmic  scale  which  accentuates 

TT  1/2  1/2 

the  distance  for  y/y^  <  1*  The  normalized  (u  )  and  (v  )  decrease 
only  slightly  with  increasing  wind  speed  but  exponentially  with  depth. 
However,  the  vertical  component  is  always  larger  than  the  horizontal 
component  at  the  same  depth.  The  profiles  for  the  two  lowest  wind-speed 
cases  are  not  plotted  because  waves  were  not  present  or  were  insignifi- 
cant.  Therefore,  neither  uq  nor  y^  could  be  computed. 

When  the  RMS  wave-induced  velocity  profiles  are  normalized  by  the 
u*  and  6  as  shown  in  Figs.  21  and  22,  the  profiles  do  not  collapse  into 
a  band.  However,  the  figures  enable  us  to  make  a  direct  comparison 
between  the  RMS  wave-induced  velocity  components  and  the  RMS  turbulent 
velocity  components.  At  the  lowest  wind  speed  there  are  essentially  no 
wave-induced  fluctuations,  so  RMS  u  <  RMS  u'  and  RMS  v  <  RMS  v*. 
At  uw  “  2.6  m/s,  the  RMS  u  is  smaller  than  the  RMS  u',  but  in  the 
region  -y/6  <  0.03,  the  RMS  v  is  larger  than  the  RMS  v'. 

At  uw  ■  3.2  m/s,  the  RMS  u  is  of  the  same  order  of  magnitude  as 
the  RMS  u*.  For  -y/6  >  0.02,  the  RMS  u'  becomes  the  dominant  component 
as  the  RMS  u  decays  rapidly  with  depth.  Furthermore,  the  ratio  of  RMS  v 


to  RMS  v'  Is  about  five-to-one  at  -y/6  -  0.01  and  the  ratio  decreases 
to  unity  at  -y/6  *  0.06,  then  the  RMS  v'  starts  to  dominate. 

For  U",  values  greater  than  3.2  m/s,  the  maximum  RMS  u  exceeds  the 
maximum  RMS  u’  by  a  factor  of  approximately  2,  while  the  maximum  RMS  v 
exceeds  the  maximum  RMS  v*  by  a  factor  of  about  3.  Of  course,  as  -y/6 
increases,  the  wave-induced  quantities  decrease  more  rapidly  than  the 
turbulent  quantities. 

4.2  MECHANICAL  WAVE  IN  EULERIAN  FRAME 
4.2.1  Water  Surface  Conditions 

A  one-hertz,  sinusoidal,  mechanically  generated  motion  was  used  as 
the  basic  propagating  wave  train.  The  actual  measured  amplitude  of  the 
one-hertz  wave  component  is  tabulated  in  Table  1;  the  amplitude  is  about 
22  mm  for  the  mechanical-wave  experiments.  Time  traces  of  the  wave- 
height-gauge  signal  for  the  five  different  wind  speeds  are  shown  in 
Fig.  23;  each  trace  is  of  a  two-second  duration. 

At  zero  wind  speed,  there  are  no  ripples  on  the  water  surface,  and 
the  wave  has  a  slightly  longer  trough  than  the  peak.  This  is  due  to  a 
slight  nonlinearity  of  the  small,  but  finite,  amplitude  mechanically  gen¬ 
erated  waves  (ka  £  0.09).  At  u*,  -  1*7  m/s,  the  waveform  is  slightly  dis¬ 
torted  from  a  pure  sinusoid,  but  the  water  surface  is  still  relatively 
smooth.  The  slightly  clipped  wave  peak  is  due  to  meniscus  effect  of  the 
wave-height-gauge  wire  and  is  not  indicative  of  the  actual  wave  shape. 
Small  wavelets  or  ripples  can  be  seen  riding  on  the  dominant  wave  at 
u,,,  ■  2.5  m/s.  The  ripples  grow  in  size  and  steepness  with  increasing  u„, 
as  can  be  seen  from  the  time  traces  for  higher  wind-speed  runs . 

Since  phase  averaging  is  used  in  analyzing  this  set  of  data,  the 
surface  elevation  t\  can  be  decomposed  as; 


The  random  component  or  the  ripple  tj'  is  the  deviation  of  the  instanta¬ 
neous  water  surface  from  the  imposed  periodic  wave  field  n«  An  indica¬ 
tor  of  the  surface  "roughness"  on  the  periodic  wave  surface  is  given 
- K  1/2  —"A  1/2 

by  (n,z)  •  The  variation  of  (n'O  with  wind  speed  u,,,  is  plotted  in 

—  H  1/2 

Fig.  24.  It  is  evident  that  (n,z)  grows  directly  with  u^,  in  this 

— j-  1/2 

wind  speed  range.  The  non-zero  value  of  (n'  )  at  u^  *  0.0  m/s 
represents  the  noise  level  of  the  wave-height  gauge. 

The  reflection  coefficient  of  the  traveling  wave  train  was  evalu¬ 
ated  from  the  signals  of  the  two  wave-height  gauges.  The  details  of  the 
method  can  be  found  in  Papadimitrakis  (1982).  The  reflection  coeffi¬ 
cient  varies  slightly  with  wind  speed,  but  is  always  less  than  5%,  which 
agrees  with  the  results  of  Bole  and  Hsu  (1967). 

4.2.2  Surface  Drift  Velocity 

The  mean  Lagrangian  and  Eulerian  surface  drift  velocities  versus  wind 
speed  are  plotted  in  Fig.  25.  The  mean  drift  velocities,  both  Lagrangian 
and  Eulerian,  for  this  set  of  experiments  are  consistently  lower  than 
those  for  wind-wave  runs  at  equivalent  wind  speeds.  As  mentioned  in 
Sec.  3.3.4,  the  timing  error  can  be  as  much  as  one  second  in  this  case, 
but  this  alone  cannot  account  for  such  consistent  results.  It  is  believed 
that  the  mechanical  wave  provides  some  sheltering  effects  from  the  wind 
along  the  wave  trough,  thus  leading  to  a  lower  mean  drift  velocity. 

4.2.3  Mean  Velocity  Profiles 

The  mean  velocity  defect  (ug  -  u)  profiles  are  shown  in  Fig.  26. 


Since  the  mechanically  generated  wave  has  an  amplitude  of  about  22  mm, 
the  closest  measuring  point  is  at  least  22  mm  from  the  mean  water  level. 
The  velocity  defect  increases  with  wind  speed.  At  each  wind  speed,  the 


defect  increases  with  depth,  except  for  the  zero  wind-speed  case.  The 
mean  Lagrangian  surface  drift  velocity  measured  for  the  zero  wind-speed 
case  is  very  close  to  zero,  but  the  computed  mean  Stokes  surface  drift 
is  not.  Therefore,  the  mean  Eulerian  surface  drift  velocity  is  nega¬ 
tive.  It  follows,  for  this  case,  that  the  mean  velocity  u  is  negative 
near  the  surface  and  goes  to  zero  deeper  down.  Thus,  the  velocity 
defect  decreases  with  depth. 

The  velocity  defect  profiles  in  wall  coordinates  are  plotted  in 
Fig.  27.  The  definitions  for  u+  and  y+  are  given  in  Sec.  4.1.3  and  the 
same  technique  was  used  to  deduce  u*.  All  the  profiles  are  logarithmic 
and  cluster  very  close  to  each  other.  There  is  a  slight  downward  shift 
of  the  profile  with  increasing  wind  speed.  Near  the  interface  (small 
-y+),  there  is  substantial  deviation  from  the  logarithmic  trend  for  the 
two  lowest  wind-speed  experiments- 

The  well-known  relationships  for  smooth  and  rough  wall  flows  are 
also  shown  in  Fig.  27  for  comparison.  The  profiles  lie  across  the  rough 
wall  flow  regime  and  have  gentler  slopes  than  the  reference  lines. 
Furthermore,  the  slopes  of  the  profiles  are  about  the  same  as  those  in 
the  high-wind-speed  wind-wave  experiments. 

The  mean  vertical  velocity  profiles  are  plotted  in  Fig.  28.  For 
the  two  low  wind-speed  experiments,  the  mean  v  are  within  ±  2  mm/s  about 
zero,  this  increases  to  ±  6  mm/s  at  higher  wind  speeds.  Therefore,  the 
mean  vertical  velocity  gradient  with  depth  is  very  much  smaller  than  the 
mean  horizontal  velocity  gradient. 

4.2.4  Mean  Turbulent  Quantities 


The  mean  turbulent  shear  stress  (-u’v')  profiles  for  all  wind  speeds 


other  than  zero  are  shown  in  Fig.  29.  The  friction  velocity  u*  and  the 


turbulent  boundary-layer  thickness  6  are  found  using  the  same  procedure 
outlined  In  Sec.  4.1. 4.  Since  u*  and  6  cannot  be  defined  for  the  zero 
wind-speed  case,  this  data  is  not  plotted,  but  is  tabulated  in  Appen¬ 
dix  D.  The  LDA  frequency  trackers  were  set  at  a  wider  tracking  range 
for  higher  wind-speed  experiments  (2. 5-6. 2  m/s)  resulting  in  a  three¬ 
fold  Increase  in  the  uncertainty  interval  for  the  mean  turbulent  shear 
stress.  In  addition,  there  is  about  -5.0  x  10“^  m^/s^  bias  in  the  mean 
turbulent  shear  stress  values.  This  small  bias  is  insignificant  for 
wind  speeds  at  4.1  and  6.2  m/s.  However,  it  is  significant  for  the 
intermediate  wind  speed  at  2.5  m/s.  The  u*  for  this  particular  case  was 
found  by  two  different  methods.  The  first  method  was  a  direct  interpo¬ 
lation  of  the  least-squares-f it  line  through  the  origin  from  Fig.  30. 
The  second  method  was  by  taking  out  the  bias  from  the  measured  values 
before  using  the  usual  linear  least-squares  procedure.  The  u*  values 
obtained  by  these  two  methods  agree.  Accordingly,  6  was  found  from  the 
second  method.  No  correction  was  necessary  for  any  other  data;  the 
actual  measured  values  are  summarized  in  Appendix  D  for  reference. 

When  uw  *  0.0,  the  measured  mean  turbulent  shear  stress  is  at  the 
noise  level  and  is  essentially  zero  throughout.  At  higher  wind  speeds, 
the  linear  distribution  of  the  mean  turbulent  shear  stress  with  depth  is 
apparent.  Contrary  to  the  wind-wave  results,  there  is  more  data  scatter 
for  higher  wind-speed  runs.  The  increase  in  data  scatter  at  higher  wind 
speed  is  partially  due  to  a  wider  tracking  range.  Also,  at  high  wind 
speeds,  the  distortion  of  the  instantaneous  wave  shape  from  the  periodic 
component  such  as  u,  v,  and  n  is  greater.  Given  the  large  amplitude  of 
the  periodic  component,  a  very  small  error  in  the  phase-averaging  pro¬ 
cess  can  lead  to  a  relatively  large  error  in  the  mean  turbulent  shear 


stress.  All  In  all,  the  mean  turbulent  shear  stress  profiles  are 
believed  to  be  representative  for  such  flow  fields. 

The  RMS  u'  and  RMS  v'  profiles  are  presented  in  Figs.  31  and  32 
using  the  same  scales  as  those  for  wind-wave  runs.  These  quantities 
behave  very  differently  from  the  wind-wave  case.  The  normalized  RMS  u’ 
and  v’  have  a  value  of  unity  for  -y/6  >  1,  then  increase  rapidly  toward 
the  interface.  Also  shown  in  the  figures  are  the  data  of  Klebanoff 
(1955).  The  RMS  u'  and  v*  deviate  radically  from  these  reference  data. 
In  addition,  there  is  no  anisotropy  between  RMS  u'  and  v'.  The  data 
collapse  for  various  wind  speeds  is  good.  Therefore,  in  this  set  of 
experiments,  u*  and  6  remain  the  prominent  parameters  for  the  turbulent 
part  of  the  flow  field. 

4.2.5  Mean  Wave-Induced  Quantities 

The  mean  wave-induced  shear  stress  (-u  v)  profiles  are  shown  in 
Fig.  33.  Similar  to  the  wind-wave  experiments,  the  mean  wave-induced 
shear  stress  is  negative.  At  u,  -  1.7  m/s,  the  magnitude  of  the  mean 
induced  shear  stress  is  about  30  times  that  of  the  mean  turbulent  shear 
stress  C-u’v’)  near  the  interface,  then  decreases  to  about  15  times  that 
value  at  the  edge  of  the  turbulent  boundary  layer  (-y  ■*  6).  For  uffi  * 

2.5  m/s,  this  ratio  is  about  15  near  the  surface  and  decreases  to  about 
10  at  a  depth  of  6.  When  the  wind  speed  was  4.1  m/s,  the  ratio  is  about 
9  near  the  interface  and  decays  to  3  at  6.  At  the  highest  wind  speed, 
this  ratio  is  7  near  the  surface  and  diminishes  to  almost  zero  at  -y  * 
6.  Therefore,  even  at  the  highest  wind  speed,  the  mean  wave-induced 
shear  stress  remains  the  dominant  component  in  the  total  mean  shear 
stress  within  the  region  6  where  turbulence  is  active. 


Let  k  be  the  wave  number  of  a  1-Hz  monochromatic  deep-water  wave. 
If  1/k  is  used  as  the  normalizing  depth,  then  the  mean  wave-induced 
stress  decays  to  about  zero  at  -ky  ■  1  as  shown  in  Fig.  34.  In  physi- 
cal  units  (omVs  ),  as  tabulated  in  Appendix  D,  the  mean  wave-induced 
shear  stress  profiles  remain  relatively  unchanged  for  u,,,  »  0.0  and 
1.7  m/s.  There  is  a  slight  increase  in  the  shear  stress  value  at 
2.5  m/ s  and  it  has  a  slightly  slower  rate  of  decay.  Then  at  uro  = 
4.1  m/s,  there  is  a  significant  increase  in  the  mean  wave-induced  shear 
stress  deeper  in  the  water  than  the  lower  wind-speed  runs.  At  the 
highest  wind  speed,  the  mean  wave-induced  shear  stress  values  are  almost 
twice  as  those  at  uM  ■  4.1  m/s.  So  in  physical  units,  there  is  a 

definite  increase  in  mean  wave-induced  shear  stress  with  wind  speed. 

— 2  1/2  — ^  1/2 

The  RMS  wave-induced  or  periodic  velocity  components  [(u  )  ,  (v  ) 

are  plotted  in  Figs.  35  and  36  with  depth  normalized  on  1/k.  The  nor¬ 
malizing  velocity  scale  is  the  theoretical  RMS  orbital  velocity  given  by 
the  1-Hz  component  of  the  water  wave  train,  namely,  2irfDng//2,  where 

A 

fjj  ■  1  Hz  and  is  the  mean  amplitude  of  the  1-Hz  component  of  the 
water  wave  spectra.  The  1//2  factor  converts  the  peak  velocity  into  an 
RMS  velocity.  Then,  using  this  velocity  scale  and  1/k  as  the  length 
scale,  the  normalized  RMS  u  and  v  profiles  collapse  into  a  very  narrow 
band.  The  data  collapse  is  excellent  for  -ky  greater  than  0.3.  At 
lower  -ky  values,  i.e.,  closer  to  the  surface,  both  the  normalized  RMS  u 
and  v  seem  to  decrease  slightly  with  increasing  wind  speed.  In  addi¬ 
tion,  the  RMS  u  is  always  less  than  the  RMS  v  at  each  wind  speed. 

The  velocity  scale  used  here  is  different  from  that  used  in  the 
72  1/2  72  1/2 

wind-wave  case  [(uq)  ].  If  (UQ)  Is  used  as  the  normalizing  parame¬ 

ter,  the  RMS  u  and  v  profiles  do  not  collapse  into  a  band,  because  the 
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definition  of  the  wave  or  periodic  velocity  component  is  different  for 

each  case.  In  the  phase-averaging  process,  only  the  1-Hz  component  and 

its  bound  harmonics  of  the  signal  for  u,  v,  q,  etc.  are  Included  in  the 

wave  or  periodic  part.  On  the  other  hand,  the  random  ripples  riding  on 

TJT  1/2 

the  main  wave  are  included  in  the  calculation  of  (u  )  ,  which,  there- 

o 

fore,  has  a  larger  value  than  2ttfQrig//2.  The  difference  between  these 

two  velocity  scales  increases  with  wind  speed  because  more  ripples  are 

Tz  1/2 

generated  at  high  wind  speeds  which  lead  to  a  higher  (uq)  while  rig 
remains  relatively  constant  for  all  wind  speeds.  Since  the  random  rip- 
pies  are  effectively  filtered  out  by  phase  averaging,  (uq)  is  not  the 
proper  velocity  scale  for  normalizing  quantities  obtained  via  phase 
averaging,  viz.  RMS  u  and  v. 

Once  again,  when  the  RMS  u  and  v  profiles  are  normalized  by  u*  and 
6,  the  profiles  do  not  coalesce  into  a  band  as  presented  in  Figs.  37  and 
38.  In  the  figures,  the  data  bands  representing  the  normalized  RMS  u' 
and  v'  profiles  are  indicated  by  the  dashed  lines.  The  RMS  u  is  almost 
an  order  of  magnitude  bigger  than  RMS  u'  at  -  1.7  m/s.  At  u^  - 

2.5  m/s,  the  ratio  of  RMS  u  to  RMS  u'  decreases  to  about  6.  Then  at 
u^  *  4.1  m/ s  and  6.2  m/s,  the  ratios  become  4  and  2,  respectively.  The 
ratios  of  RMS  v  to  RMS  v'  are  slightly  higher  than  those  for  the  hori¬ 
zontal  component  because  RMS  v  is  always  bigger  than  RMS  u.  It  is  quite 
apparent  that  the  RMS  u  and  v  dominate  over  the  whole  region  even  though 
the  RMS  u'  and  v'  become  more  significant  at  high  wind  speeds. 

4.2.6  Wave  Perturbation  Velocities 

The  wave  perturbation  velocities  u  and  v  were  found  by  the  phase¬ 
averaging  technique.  A  typical  result  for  u  and  v  is  shown  in  Fig.  39. 
Also  shown  in  the  figure  is  the  phase  average  for  the  water  surface 


elevation  in  the  same  run.  It  is  clear  that  the  periodic  velocities  are 
nearly  sinusoidal  with  the  same  frequency  as  that  of  the  progressive 
water  waves.  The  amplitudes  (u,  v)  and  the  phase  angles  ( 9 - ,  0^),  with 
respect  to  n,  of  the  first  harmonic  of  the  periodic  velocities  were 
obtained  by  spectral  analysis.  This  technique  is  fully  documented  in 
Hsu  et  al .  (1977).  The  total  harmonic  content  of  the  phase  average 
velocity,  in  terms  of  energy  is  given  by: 


mean  square  total  harmonic  content 


where  i  *  1  or  2.  The  total  harmonic  energy  content  for  u  and  v  at 
each  depth  in  an  experiment  are  tabulated  in  Appendix  D.  The  u  usually 
has  more  harmonic  components  than  v  near  the  surface  but  the  difference 
decreases  with  depth  until  the  harmonic  content  becomes  almost  identical 
for  both  u  and  v.  For  the  worst  case,  the  total  harmonic  energy  content 
for  u  or  v  is  at  most  2.6%.  Therefore,  it  is  sufficient  to  present  the 
amplitude  and  the  phase  of  the  periodic  velocity  at  the  excitation  fre¬ 
quency.  Although  the  full  results  for  u  and  v  at  each  data  point  were 
graphed,  they  are  too  numerous  to  present  here.  Readers  are  referred  to 
a  separate  report  and  files  on  800  BPI  magnetic  tapes  available  at  the 
Stanford  University  Environmental  Fluid  Mechanics  Laboratory. 

The  plots  of  the  amplitudes  (u,  v)  and  the  phase  lags  ( 0 ~ ,  9^), 
relative  to  n,  versus  -ky  are  given  in  Fig.  40.  The  amplitudes  are 
normalized  by  the  theoretical  peak  orbital  velocity  given  by  the  1-Hz 
component  of  the  water  wave  train,  namely  Ztif^rig.  Both  the  u  and  v 


profiles  show  exponential  decay  with  depth.  Also,  the  u  is  always 


smaller  than  v  for  this  set  of  experiments.  In  the  proximity  of  the 


interface,  there  Is  a  slight  decrease  in  u  and  v  with  increasing  wind 
speed.  This  decrease  is  most  prominent  by  comparing  profiles  for  u#  = 

1.7  m/s  and  u0  *  2.5  m/s.  Some  of  these  trends  are  also  found  in 

1/2  -Zty  1/2 

plots  of  (uz)  and  (v^)  versus  -ky,  but  they  do  not  give  any  phase 
relationships.  The  numerical  values  for  the  phase  lags  (9q,  8^)  can  be 
found  in  Appendix  D.  The  following  observations  on  the  phase  lags  can 
be  made  Irrespective  of  wind  speeds  and  depth: 

1.  The  phase  lag  of  u  with  respect  to  the  phase  of  n<j  is 
very  close  to  0  but  not  exactly  so.  The  phase  lag  varies 
most  between  -3°  and  +2°  about  0  for  all  data  points. 

2.  The  phase  difference,  0-  -  0^,  is  always  slightly  less 

than  90°.  Specifically,  the  phase  difference  remains 

relatively  constant  at  about  89°  except  at  u,,,  ■*  6.2  m/s 
where  the  phase  difference  decreases  to  about  86°. 

The  velocity  amplitudes  and  their  respective  phases  were  also  found 
by  the  cosine  curve-fitting  techniques  (see  Appendix  C).  The  spectral 
and  cosine-fitting  methods  gave  virtually  identical  results:  the  ampli¬ 
tude  values  agreed  to  within  ±  0.5%  and  ±  0.1°  for  the  phase  angles.  So 
the  cosine-fitting  technique  verified  the  results  of  the  spectral  method. 

4.2.7  Wave-Induced  Stresses 


The  wave-induced  turbulent  stresses  are  defined  as  r^j  «  uiuj  and 


were  obtained  by  the  phase-averaging  method.  In  general,  the  magnitudes 
of  rjj  are  very  small,  so  the  signal-to-noise  ratio  is  much  lower  than 
that  of  u  or  v.  The  phase-averaged  results  were  used  to  calculate  the 


amplitudes  r^j  and  the  phase  angles  8^^  at  the  excitation  frequency. 


There  could  be  substantial  harmonic  content  in  r.,,  but  only  r<  ,  and 


8.  are  considered  here.  The  numerical  values  for  r.  ,  are  stored  on 
fij 

800  BPI  magnetic  tapes  and  all  the  plots  are  collected  in  a  separate 
report. 

The  rj.  and  0*  values  are  given  in  Appendix  D.  They  are  also 
J  1  j 

plotted  in  Fig.  41.  The  amplitudes  are  normalized  by  the  square  of  the 
friction  velocity  in  these  plots.  Because  the  friction  velocity  is  not 
defined  in  the  zero  wind-speed  experiment,  its  results  are  not  shown. 

A  A 

At  each  depth,  the  r^  and  ^2  are  approximately  equal  in  magnitude 
regardless  of  wind  speeds.  Furthermore,  they  decrease  quickly  with 
depth.  For  wind  speeds  less  than  3  m/s,  the  normalized  amplitudes  drop 
to  the  background  noise  level  for  -ky  greater  than  0.2.  At  higher  wind 

A 

speeds,  this  number  increases  to  0.3.  On  the  other  hand,  r^  is  much 

A  A 

smaller  than  either  r^  or  ^2  and  remains  near  the  noise  level  for 
practically  all  wind  speeds  and  depths. 

The  phase  lags  8f^  stay  very  close  to  180°  in  experiments  with 
wind  speeds  less  than  4  m/s.  At  u,,,  *  4.1  m/s,  there  is  a  fair  amount 
of  data  scatter,  but  the  overall  trend  of  the  phase  angles  stay  in  the 
neighborhood  of  180°.  At  the  highest  wind  speed,  the  phase  lags  do  not 
have  any  trend  with  depth.  The  apparent  data  scatter  at  high  wind 
speeds  may  be  due  to  the  large  amount  of  harmonics  present,  thus  making 
it  difficult  to  determine  the  phase  lags  at  the  first  harmonic. 

The  wave-induced  wave  stresses  are  defined  as  UjUj.  They  can  be 
found  by  multiplying  u^  and  Uj  over  the  period  of  the  main  wave  and 
removing  the  mean  wave-induced  shear  stresses  ujuj  from  the  products. 
It  is  shown  in  Sec.  4.2.6  that  u^  is  almost  sinusoidal,  so  we  can  char¬ 
acterize  u^  as  follows: 

*»  A 

u^  ■  u^  cos(u>t  -  8j)  +  harmonics 


Therefore,  vl^Qj  becomes, 

uiuj  “  ~2~*  (c°8(2wt  -  8^  -  9j)  +  cos(9^  -  9j)]  +  harmonics  (4-6) 

Then  u^uj  Is  just  simply  the  time-varying  part  of  Eq.  4-6.  It  i9  obvi¬ 
ous  that  u^Uj  has  a  dominant  frequency  twice  that  of  the  excitation  fre- 

A  A  w 

quency  and  has  an  amplitude  of  u^Uj/2.  The  harmonic  content  of  ut  is 
very  small,  consequently  the  harmonic  content  of  u^uj  relative  to  the 
dominant  frequency  is  also  negligible.  For  this  set  of  experiments, 

A  A  A 

u^uj/2  is  very  much  larger  than  r^j.  The  implication  of  this  can  only 
be  resolved  by  further  analysis  of  the  effects  of  measurements  made  in  a 
non-wave-following  coordinate  system  and  the  governing  equations  (see 
Chap.  5). 


4.3  MECHANICAL  WAVE  IN  WAVE-FOLLOWING  EULERIAN  FRAME 

4.3.1  Water  Surface  Conditions 

A  one-hertz,  sinusoidal,  mechanically  generated  wave  was  also  used 
in  this  set  of  experiments.  The  measured  amplitudes  of  the  one-hertz 
wave  components  are  given  in  Table  1.  Time  traces  of  the  wave-height- 
gauge  signal  under  different  wind  speeds  are  shown  in  Fig.  42,  with  each 
trace  being  two  seconds  in  duration.  The  variation  of  surface  "rough- 

— j  1/2 

ness"  (n'^)  with  wind  speed  u*,  is  plotted  in  Fig.  43.  The  wind 
speeds,  wave  shapes,  and  reflection  coefficients  of  the  traveling  wave 
train  are  virtually  identical  to  those  of  the  fixed-frame  mechanical- 
wave  experiments. 

4.3.2  Surface  Drift  Velocity 

The  mean  Lagrangian  and  Eulerian  surface  drift  velocities  versus 
wind  speeds  are  shown  in  Fig.  44.  The  results  agree  with  those  for  the 


fixed-frame  mechanical-wave  case;  thus,  the  remarks  made  for  those 
experiments  are  applicable  here  also. 

4.3.3  Mean  Velocity  Profiles 

Under  the  proposed  coordinate  transform  (Sec.  2.4.1),  the  depth  is 
now  -y*  Instead  of  -y.  However,  the  mean  position  of  the  LDA  probe 
volume  when  observed  in  the  laboratory  frame  is  y*  =  y.  Then,  as  noted 
in  Sec.  2.5,  the  interpretation  of  the  measured  quantities  in  the  wave¬ 
following  frame  is  different  from  that  of  the  fixed  frame. 

The  mean  velocity  defect  (Ug  -  u)  profiles  are  presented  in  Fig.  45, 
showing  that  the  mean  velocity  defect  increases  with  wind  speed.  Simi¬ 
lar  to  the  fixed-frame,  mechanical-wave  experiments,  the  computed  mean 
Gulerian  surface  drift  velocity  and  the  measured  mean  velocity  are 
negative  at  zero  wind  speed.  Most  importantly,  we  see  that  the  wave¬ 
following  system  enables  measurements  very  close  to  the  wavy  surface. 
For  low  wind-speed  runs,  where  the  ripples  on  the  one-hertz  dominant 

wave  were  small,  it  was  possible  to  reach  a  y*  of  -1  to  -2  mm.  Even  at 

—7 

the  highest  wind  speed,  where  (n,z)  was  about  2  mm,  the  probe  was 
still  able  to  make  measurements  at  y*  ■  -5.6  mm. 

It  is  important  to  stress  that,  although  the  measured  velocity  is 
in  a  wave-following  frame,  the  velocities  are  still  Eulerian  quantities. 
Therefore,  the  mean  Eulerian  surface  drift  is  used  in  computing  the  mean 
velocity  defect.  Of  course,  the  measured  mean  Eulerian  surface  drift 
is  perfect  for  use  in  this  wave-following  coordinate  system  since  the 
velocity  was  obtained  by  measuring  the  mean  Lagrangian  surface  drift 
exactly  at  y*  m  0,  and  then  removing  the  mean  Stokes  surface  drift. 
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The  mean  velocity  defect  profiles  in  wall  coordinates  are  shown  in 
Figs.  46  and  47.  Each  figure  uses  a  different  set  of  u*  values.  The 
uA's  obtained  from  the  fixed-frame  experiments  are  used  for  normaliza¬ 
tion  in  Fig.  46.  The  u*'s  obtained  by  stress  continuity  at  the  inter¬ 
face  from  the  experiments  of  Hsu  et  al.  (1983)  are  used  in  Fig.  47.  A 
discussion  of  the  validity  of  using  either  of  these  u*  values  is  given 
in  the  next  section. 

In  Fig.  46,  the  two  low  wind-speed  profiles  lie  close  together  and 
fall  between  the  smooth-wall  and  rough-wall  data  lines  with  a  strong 
bias  towards  the  latter  line.  Near  the  interface,  the  profiles  follow 
the  u+  ■  -y*+  line.  Then,  the  profiles  become  logarithmic  for  -y*+ 
greater  than  20.  For  the  two  high  wind-speed  cases,  the  profiles  also 
cluster  together,  but  lie  below  those  for  the  lower  speed  runs,  as  well 
as  below  the  data  line  representing  the  fully  rough  flow  regime  on  solid 
walls.  Furthermore,  the  profiles  exhibit  logarithmic  trends  for  -y*+ 
greater  than  50. 

Results  similar  to  the  above  are  observed  in  Fig.  47  where  a  dif¬ 
ferent  set  of  u*  was  used  in  normalization.  The  two  low-speed  mean 
velocity  defect  profiles  now  lie  almost  totally  between  the  smooth  and 
fully  rough  flow  regime.  However,  the  two  high  wind-speed  profiles 
remain  relatively  unchanged. 

The  mean  vertical  velocity  profiles  are  shown  in  Fig.  48.  The  mean 
velocity  lies  within  ±  5  mm/s  about  zero  for  wind  speeds  from  0  to 
4.1  m/s.  At  the  highest  wind  speed,  the  mean  vertical  velocity  is 
positive  and  is  about  2  to  9  mm/s  for  the  whole  profile.  Again,  the 
mean  vertical  velocity  gradient  (with  depth  -y*)  is  small  compared  to 
the  mean  horizontal  velocity  gradient. 


4.3.4  Mean  Turbulent  Quantities 

Due  to  data  scatter,  well-defined  friction  velocities  (u*)  cannot 
be  found  from  the  measured  mean  turbulent  shear  stress  via  the  usual 
curve-fitting  method  as  used  in  the  two  previous  sets  of  experiments. 
However,  good  estimates  of  u*  can  be  obtained  by  the  following  methods: 


Method  1.  For  the  two  low  wind-speed  experiments,  the  first  point 
(y*  ~  0)  of  the  velocity  profile  may  lie  in  the  viscous  sublayer  where 
the  velocity  gradient  is  linear.  Then,  the  friction  velocity  can  be 
estimated  as: 


+  *+ 
u  ■  -y* 


u„  -  u  -y*u 

C  *  4 


2  US  -  U  I 

u*  *  v  -=yr . 

J  '  1st  data  point 

The  u*  values  estimated  by  this  method  for  the  two  low-speed  runs  are 
given  in  Table  8. 

Method  2.  Hsu  et  al.  (1983)  found  friction  velocities  on  the  air 
side  in  a  transformed  coordinate  system  under  similar  experimental 
conditions.  Utilizing  the  u*  values  for  the  air  flow  and  assuming 
stress  continuity  across  the  interface,  we  can  find  the  interpolated  or 
extrapolated  friction  velocities  for  the  water  flow  given  in  Table  8. 

Method  3.  The  u*  values  obtained  in  the  mechanical-wave  fixed- 
frame  experiment  may  be  good  estimates  for  the  u*  values  in  a  trans¬ 
formed  coordinate  system  and  these  are  also  given  in  Table  8  for  com¬ 
parison. 


It  can  be  seen  from  Table  8  that  for  the  two  low  wind-speed  runs, 
there  is  little  difference  between  the  u*'s  obtained  by  Methods  1  and 


3.  Although  the  u*'s  given  by  Method  2  are  lower  than  those  by 
Method  3,  especially  for  the  lowest  wind-speed  run,  they  are  of  compar¬ 
able  magnitude.  The  low  u*  value  for  Method  2  at  1.7  m/s  may  be  due  to 
the  curve-fitting  procedure  used  by  Hsu  et  al.  (1977).  As  pointed  out 
by  Papadimitrakis  (1982),  the  u*  values  can  vary  as  much  as  10%  by  using 
different  fitting  methods,  and  each  method  is  very  sensitive  to  the 
fitting  region.  At  higher  wind  speeds,  Method  1  is  not  applicable 
because  of  large  random  ripples  present  on  the  main  wave.  Method  2  now 
yields  higher  u*  values  than  Method  3.  If  the  difference  is  truly 
significant  and  not  a  result  of  uncertainty,  then  the  difference  might 
be  attributed  to  the  fact  that  a  portion  of  the  momentum  from  the  wind 
goes  into  wave  growth,  thus,  less  shear  stress  is  measured  in  the  water. 

The  friction  velocities  given  by  the  mechanical-wave,  fixed-frame 
experiments  are  likely  to  be  the  best  estimates  for  u*  for  the  present 
case  also.  Therefore,  they  are  used  throughout  for  all  mechanical-wave 
experiments  unless  stated  otherwise.  In  any  case,  it  is  true  that  none 
of  the  interpretations  of  the  data  set  are  changed  significantly  by  the 
estimates  for  u*. 

For  completeness,  the  measured  mean  turbulent  shear  stresses  in  the 
transformed  coordinate  system  are  tabulated  in  Appendix  D.  They  are 
also  shown  in  normalized  form  in  Fig.  49.  At  zero  wind  speed,  the  mean 
turbulent  shear  stress  is  very  small  and  is  of  the  order  of  5.0  ^ 
10  m  /s  ,  which  is  very  close  to  the  resolution  limit  of  the  LDA. 
There  is  quite  a  bit  of  data  scatter  as  can  be  seen  in  Fig.  49.  In 
addition,  the  normalized  mean  turbulent  shear  stress  does  not  go  to  one 
near  the  interface  but  is  about  15  for  u^,  *  1.7  m/s,  5  for  u,,,,  *  2.5 
and  4.1  m/s,  and  roughly  unity  at  uw  ■  6.1  m/s.  This  result  suggests 
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the  shear  stress  may  vary  over  the  main  wave  with  possibly  a  minimum  at 
the  wave  trough. 


The  RMS  u'  and  v'  are  normalized  by  the  wave  parameters  and  shown 
in  Figs.  50  and  51.  The  wave  parameters  are  used  here  because  they  have 
less  uncertainty.  The  RMS  u'  and  v'  both  increase  with  wind  speed  and 
decrease  with  depth.  There  is  little  anisotropy,  except  at  the  two 
higher  wind  speeds  where  the  RMS  v'  is  slightly  higher  near  the  inter¬ 
face  and  lower  than  RMS  u'  away  from  the  surface.  Furthermore,  there  is 
no  data  collapse  because  the  wave  number  k  and  the  RMS  orbital  velocity 
at  the  surface  are  not  the  appropriate  scales  for  turbulent  processes. 

Since  the  mechanical-wave  fixed-frame  u*'s  are  used  in  this  case, 
the  corresponding  turbulent  boundary-layer  thicknesses  (5)  may  also  be 
suitable  length  scales  to  use.  Figures  52  and  53  show  the  RMS  u'  and  v’ 
normalized  by  u*  and  6.  Data  collapse  is  fair,  thus  implying  that  u* 
and  5  may  also  be  the  appropriate  scales  for  the  turbulent  processes  in 
the  transformed  coordinates.  Comparing  Figs.  52  and  53  with  Figs.  31 
and  32  for  the  fixed-frame  case,  one  sees  that  the  wave-following 
results  appear  as  extensions  to  those  of  the  fixed-frame  experiments  to 
shallower  depths. 


4.3.5  Mean  Wave-Induced  Quantities 


of  a  lower  magnitude  in  general,  the  trends  are  similar  to  those  for  the 
fixed-frame  experiments. 

The  RMS  wave-induced  velocity  components  are  plotted  in  normalized 
coordinates  in  Figs.  55  and  56 <  Data  collapse  is  excellent  despite  a 
slightly  larger  scatter  for  the  RMS  u.  Again,  RMS  v  is  always  larger 
than  RMS  u  at  equivalent  wind  speeds  and  depth.  By  comparing  Figs.  55 
and  56,  and  Figs.  35  and  36  for  the  fixed-frame  case,  one  sees  again 
that  the  wave-following  results  appear  as  extensions  to  the  fixed-frame 
plots  to  shallower  depths.  The  normalized  RMS  u  approaches  a  value  of 
about  0.9,  and  the  normalized  RMS  v  attains  a  limiting  value  of  about 
unity  near  the  interface. 

4.3.6  Wave  Perturbation  Velocities 

Similar  to  the  fixed-frame  experiments,  the  wave  perturbation 
velocities  u  and  v  are  almost  sinusoidal.  The  normalized  amplitude  and 
the  phase  angles  of  the  periodic  velocities  at  the  excitation  frequency 
are  shown  in  Fig.  57.  In  each  experiment,  the  amplitudes  decay  exponen- 

A 

tially  with  depth.  The  vertical  amplitude  component  v  is  always  larger 

A 

than  the  horizontal  component  u,  especially  for  non-zero  wind  speed 

A  A 

runs.  Both  u  and  v  decrease  slightly  with  wind  speed  at  equivalent 

A  A 

depths;  the  wind  seems  to  have  a  lesser  effect  on  v  than  on  u.  The 

A  A 

normalized  profiles  for  u  or  v  approach  unity  as  -ky*  goes  to  zero 
regardless  of  wind  speeds.  Even  though  this  trend  is  also  observed  in 
fixed-frame  experiments,  it  is  more  vividly  shown  here  in  the  wave- 
following  coordinates. 

The  total  harmonic  energy  content  for  u  and  v  at  each  depth  are 
summarized  in  Appendix  D.  Both  u  and  v  have  higher  harmonic  content 
near  the  interface  than  at  deeper  depths.  Furthermore,  the  u  usually 
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has  more  harmonic  components  than  v,  except  at  u,,,  -  2.5  m/s  where  the 
reverse  Is  true.  At  greater  depths,  both  u  and  v  have  roughly  the  same 
amount  of  harmonic  energy  content  of  about  1.4%. 

A  A 

The  phase  lag  of  u  with  respect  to  the  phase  of  rig  varies  at  most 
between  -3°  and  +2°  In  the  neighborhood  of  zero  at  all  depths  except  at 
points  where  there  are  obvious  data  scatter.  The  data  scatter  was 
partially  due  to  the  lower  signal-to-noise  ratio  of  the  backward-scatter 
collection  method  for  the  LDA;  especially  at  deeper  depths  where  u  and  v 
were  smaller  than  near  the  surface.  These  data  points  are  not  discarded 
because  even  though  the  phase  angles  may  have  a  larger  error,  other 
statistical  quantities  are  still  valid.  The  phase  difference,  6-  -  0^ 
stays  at  90°  ±  2°.  The  decrease  of  the  phase  difference  at  u,,,,  = 
6.2  m/s  present  in  the  fixed-frame  experiment  is  not  found  here. 

The  cosine-fitting  technique  was  also  used  to  determine  the  velo¬ 
city  amplitudes  and  phase  angles.  Similar  to  the  fixed-frame  experi¬ 
ments,  the  results  agree  with  those  given  by  the  spectral  method.  The 
plots  of  all  u  and  v  are  collected  In  a  separate  report  and  their  numer¬ 
ical  values  are  stored  on  800  BPI  tapes. 

4.3.7  Wave-Induced  Stresses 

The  amplitudes  r^j  (normalized  by  u£)  and  the  phase  angles  of 

ij 

the  wave-induced  turbulent  stresses  r^j  at  the  excitation  frequency  are 
shown  in  Fig.  58.  There  is  some  data  scatter  and  a  few  extraneous 
points  are  eliminated  in  these  plots,  but  their  numerical  values  are 
retained  in  Appendix  D  for  completeness. 

A 

At  zero  wind  speed  (not  plotted),  the  r^j  are  very  small  and  most 
likely  represent  the  background  noise  level  rather  than  any  physical 

A  A 

quantities.  At  u„  ■  1.7  and  2.5  m/s,  rn  and  roo  are  roughly  of  equal 
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magnitudes  while  is  usually  smaller  than  both  of  them.  The  all 
decay  rapidly  with  depth  to  the  noise  level  for  -ky*  greater  than  0.3. 
At  u^  -  4.1  and  6.1  m/s,  the  r^j  are  quite  high  near  the  surface  but 
decrease  quickly  to  the  noise  level  for  -ky*  bigger  than  0.3.  There  is 

A 

no  systematic  increase  of  r^j  with  wind  speed  as  observed  in  the  fixed- 
frame  experiments.  The  apparent  increase  in  r^j  at  larger  -ky*  is 
attributed  to  data  scatter.  The  data  scatter  is  partly  due  to  a  lower 
signal-to-noise  ratio  in  the  backward-scatter  LDA  system  than  in  the 
forward-scatter  system.  Furthermore,  the  mechanically  generated  waves 
were  highly  distorted  at  such  high  wind  speeds;  a  small  error  in  obtain¬ 
ing  u  or  v  due  to  a  lower  signal-to-noise  ratio  can  contribute  to  a 
significant  error  in  r^j  . 

At  u,,,  ■  0.0  and  1.7  m/s,  the  phase  lag  8eems  to  stay  around 

180°  while  0*  and  9-  are  equal  to  0  near  the  surface  and  reach  180° 
fll  “22 

at  large  -ky*.  Similar  trends  are  also  observed  at  u.  «  2.5  m/s 
despite  large  amounts  of  data  scatter.  For  the  two  high-speed  experi¬ 
ments,  seem  to  start  from  360°  (or  0°)  for  small  -ky*  and  decrease 

gradually  to  around  90°  at  large  -ky*. 

The  comments  about  the  wave-induced  wave  stresses  uiuj  in 

A  A 

Sec.  4.2.7  are  also  valid  here.  In  general,  u^Uj/2  is  also  very  much 

A 

larger  than  r. 


.*■ » m  •  *.  «**  jy  >  *  . 


•.  s  *  *_  • . 


•-V-*  yv-'-  • 


Table  8 


Determination  of  Friction  Velocities 
for  Wave-Following  Frame  Experiments 


Method  1 

Method  2 

Method  3 

From  u+  *  -y*+ 

From  Hsu  et  al.  (1983) 

From  Fixed-Frame  Data 

(m/s)  u*  (mm/s) 

u*  (mm/s) 

u*  (mm/s) 

CHAPTER  5 


DISCUSSION 

The  representative  results  of  the  experiments  were  presented  in 
Chap.  4.  The  experiments  provide  a  comprehensive  data  set  for  the  flow 
field  directly  beneath  an  air-water  interface  in  a  well-defined  test 
facility  for  a  wide  range  of  wind  speeds,  surface  conditions,  and  refer¬ 
ence  coordinate  systems.  The  objectives  of  obtaining  velocity  data  have 
been  accomplished.  However,  it  is  crucial  to  sort  out  the  implications 
of  these  data  and  to  understand  the  physical  reasons  for  the  observed 
behavior  of  the  flow.  This  is  the  main  goal  of  this  chapter.  Whenever 
possible,  the  data  interpretation  will  be  compared  to  some  simple  analy¬ 
ses  of  the  governing  equations  or  well-established  results. 

5.1  CHARACTERISTIC  SCALES 

5.1.1  Case  I:  Wind  Wave  in  Eulerlan  Frame 

In  this  set  of  experiments,  the  air-water  interface  was  subjected 
to  the  shearing  action  of  the  turbulent  air  flow.  Thus,  shear  layers 
developed  adjacent  to  the  interface.  The  flow  was  turbulent  in  the 
water  boundary  layer  even  at  the  lowest  wind  speed.  The  definition  of 
turbulence  is  based  upon  the  highly  diffusive  nature  of  the  flow  (other 
than  due  to  molecular  diffusion).  This  was  observed  by  injecting  small 
amounts  of  dye  in  the  water  layer  during  some  preliminary  experiments. 
It  is  important  not  to  confuse  the  random  wind-generated  waves  with  tur¬ 
bulence.  The  random  waves  are  not  in  turbulent  motion  because  they  are 
essentially  irrotational  •  In  addition,  the  waves  are  mostly  nondissipa- 
tive  (though  they  often  are  dispersive),  while  turbulence  is  essentially 
dissipative  (Tennekes  and  Luraley,  1972).  This  is  the  main  reason  why  we 


can  decompose  the  flow  field  Into  three  parts:  the  mean,  Lhe  wave- 
induced  component,  and  the  fluctuating  or  turbulent  component.  Further¬ 
more,  it  may  be  possible  to  derive  some  scales  or  parameters  which  char¬ 
acterize  each  of  these  components. 

For  wind  speeds  less  than  3  m/s,  the  water  surface  remained  rela¬ 
tively  smooth  even  though  capillary-type  disturbances  were  present.  The 
wave-induced  quantities  were  expected  to  be  very  small.  Therefore,  the 
water  boundary  layer  developed  would  have  features  similar  to  those  of 
flows  over  flat  plates.  Of  course,  it  was  necessary  to  convert  the 
water  velocity  into  the  velocity  defect  form  (ug  -  u)  in  such  compari¬ 
sons.  When  the  mean  velocity  defect  profiles  are  plotted  in  wall  coor¬ 
dinates  (Fig.  12),  they  agree  with  profiles  of  flows  over  rigid  sur¬ 
faces.  The  profiles  are  logarithmic  which  is  usually  taken  as  evidence 
for  the  existence  of  the  inertial  sublayer.  In  addition,  the  profiles 
extend  to  follow  the  u+  ■  -y+  line,  which  suggests  the  existence  of  the 
viscous  sublayer.  Hence,  the  velocity  and  length  scales  appropriate  for 
flows  over  solid  surfaces  are  also  applicable  for  the  water  boundary 
layer  at  low  wind  speeds. 

In  the  inertial  sublayer,  the  characteristic  velocity  is  the 
friction  velocity  u*  and  the  characteristic  length  is  proportional  to 
the  distance  from  the  still  water  level  (i.e.,  -<y)  provided  that 
-yu  N  »  1  and  -y/<5  <<  1.  Of  course,  in  the  viscous  sublayer,  the 
viscous  shear  stress  dominates  and  the  contribution  of  the  turbulent 
shear  stress  becomes  insignificant.  So  in  the  viscous  sublayer,  the 
length  scale  becomes  v/u*,  while  u*  remains  as  the  correct  velocity 
scale.  In  the  outer  layer,  where  the  "mean"  flow  dominates  (i.e.,  u  -  0, 
Ug  -  u  -  Ug),  the  obvious  choice  of  the  velocity  scale  is  the  mean 
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Eulerian  surface  drift  velocity  Ug.  The  corresponding  length  scale  Is 
the  turbulent  boundary-layer  thickness  6.  In  addition,  the  fetch  L  can 
also  be  a  significant  length  scale. 

At  high  wind  speeds  (>  3  m/s),  waves  were  present  at  the  interface. 
The  mean  velocity  profiles  in  wall  coordinates  are  also  logarithmic. 
However,  the  slopes  of  some  of  the  profiles  are  different  from  those  at 
low  wind  speeds.  Furthermore,  there  is  an  apparent  change  of  slope  in 
some  of  the  profiles,  such  as  those  at  wind  speeds  between  3.2  and 
9.9  m/s.  There  is  no  change  of  slope  in  the  profile  at  the  highest  wind 
speed  but  the  slope  of  the  profile  is  clearly  different  from  the  low 
wind-speed  experiments.  The  possible  explanation  of  such  behavior  is 
discussed  in  Sec.  S.2.  Here,  we  focus  our  attention  on  finding  the 
relevant  parameters  characterizing  the  flows. 

The  inertial  sublayer  most  likely  persists  at  high  wind  speeds 
because  of  the  logarithmic  trend  of  the  velocity  profiles.  Conse¬ 
quently,  the  friction  velocity  u*  and  -icy  are  the  proper  scaling  velo¬ 
city  and  length  in  this  layer.  However,  we  must  keep  in  mind  that 
within  the  Inertial  sublayer,  the  wave-induced  quantities  may  affect  the 
flow  dynamics  in  the  sublayer.  In  the  outer  layer,  where  the  effects  of 
turbulence  and  waves  are  weak,  the  mean  Eulerian  surface  drift  us  and 
the  turbulent  boundary-layer  thickness  6  remain  the  appropriate  scales 
for  the  outer  mean  flow. 

The  turbulent  intensities  or  RMS  u'  and  RMS  v'  profiles  when  nor¬ 
malized  by  u*  and  6  coalesce  into  narrow  bands  (Figs.  16  and  17). 
Clearly,  u*  and  6  are  the  characteristic  scales  relating  turbulent 
intensities  and  their  decay  with  depth.  The  fact  that  6  is  the 
appropriate  length  scale  is  intuitively  obvious,  because  the  largest 


turbulent  eddy  is  probably  on  the  order  of  the  thickness  of  the  turbu¬ 
lent  boundary  layer.  Similarly,  because  u*  is  derived  from  the  turbulent 
shear  stress,  which  is  directly  responsible  for  turbulence  generation, 
the  turbulent  intensities  should  scale  with  u*  as  shown  here. 

So  far,  the  relevant  scales  for  the  mean  and  turbulent  parts  of  the 
flows  are  shown  to  be  similar  to  those  for  wall-bounded  shear  flows.  It 
is  well  known  that  wall-bounded  shear  flows  have  multiple  scales.  Here, 
the  wind-generated  waves  complicate  the  problem  further  by  introducing 
additional  velocity  and  length  scales.  The  wind-generated  waves  are 
usually  considered  irrotational  and  of  small  amplitude.  Thus,  some 
suitable  scales  for  the  wave  field  may  be  derived  from  the  small- 
amplitude  linear  water-wave  theory.  Some  obvious  velocity  scales  are 
the  phase  velocity  of  the  dominant  wave  or  the  RMS  wave  orbital  velocity 
at  the  surface.  The  phase  velocity  is  always  bigger  than  the  mean  sur¬ 
face  drift  velocity  in  the  experiments.  Hence,  there  is  no  counterpart 
of  the  critical  layer,  which  is  known  to  exist  on  the  air  side  of  the 
interface,  in  the  water.  In  addition,  the  phase  velocity  does  not 
relate  directly  with  the  measured  water  orbital  velocity.  Therefore, 
the  phase  velocity  is  not  an  important  velocity  scale  in  analyzing  the 
water  motion.  On  the  other  hand,  the  linear  water-wave  theory  can  be 

used  to  compute  the  orbital  velocity  at  the  surface  from  the  surface 

Z2  1/2 

displacement  (Eq.  4-4).  This  RMS  orbital  velocity  (uq)  is  closely 
related  to  the  measured  wave-induced  velocities.  Associated  with  this 
computed  RMS  orbital  velocity  is  the  integrated  wave  number  (1/y^,  see 
Eq.  4-5).  The  excellent  data  collapse  resulting  by  using  these  scales 
(Figs.  19  and  20)  suggests  that  they  are  the  characteristic  scales  for 
the  wave-induced  quantities. 
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Another  length  scale  of  Interest  is  the  RMS  wave  height  (n  ) 

This  length  scale  may  give  some  insight  into  the  surface  roughness  as 
defined  in  boundary- layer  type  flows  (Schlichting,  1979).  The  velocity 
profiles  in  wall  coordinates  show  a  downward  shift  with  increasing  wind 
speed  which  is  usually  associated  with  changes  in  surface  roughness. 
Kays  and  Crawford  (1980)  gave  a  law  of  the  wall,  which  incorporates  a 
roughness  Reynolds  number,  for  fully  rough  flows  with  no  pressure  gradi¬ 
ent  and  no  transpiration: 


u+  *  -  In  y+  +  B 
< 


(5-1) 


_  1  ,  32.6 

B  -  —  In  - - 

k  Re, 


Rev 


u*kS 


kg  ■  characteristic  length  of  the  roughness  element 


If  the  RMS  wave  height  is  the  only  parameter  in  determining  the  "surface 
roughness,"  then  ks  should  be  proportional  to  the  RMS  wave  height,  i.e., 

— j  1/2 

kS  "  c(n  )  for  some  constant  c.  Since  the  velocity  profile  for 

uw  ■  3.2  m/s  lies  on  the  line  corresponding  to  the  start  of  the  fully 
rough  flow  regime,  it  was  used  to  calibrate  the  simple  model.  The  con¬ 
stant  c  was  found  to  be  13.5.  Next,  the  various  values  of  B  were  calcu¬ 
lated  for  higher  wind  speeds.  The  relationships  given  by  Eq.  5-1  were 
compared  with  the  measured  profiles.  These  relationships  gave  only  very 
rough  correlation  with  the  data.  One  apparent  reason  is  that  measured 
profiles  have  a  different  slope.  More  importantly,  similar  to  the  air 
flow,  the  actual  roughness  most  likely  does  not  depend  on  the  entire 
wave  spectrum;  hence,  the  total  RMS  wave  height  should  not  be  included 
in  such  calculation  as  suggested  by  Kitalgorodskli  (1973),  Kondo  (1973), 
and  Phillips  (1977). 
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The  RMS  wave  height  actually  has  been  considered  as  a  length  scale 

1/2 

in  calculating  the  surface  RMS  orbital  velocity  (u  )  ,  because  it  is 
obtained  from  the  surface  elevation  spectrum.  It  may  seem  paradoxical 
that  the  RMS  wave  height,  which  is  a  wave-related  quantity,  would  have 
an  effect  on  the  surface  roughness,  which  is  related  to  turbulence  gen¬ 
eration.  On  close  examination,  this  is  not  at  all  unreasonable,  because 
the  waves  may  in  fact  have  some  effects  on  surface  roughness  and  give 
rise  to  wave-induced  turbulence. 

5.1.2  Case  II:  Mechanical  Wave  in  Eulerian  Frame 

The  main  distinction  between  this  case  and  the  wind-wave  case 
(Case  I)  is  the  presence  of  the  mechanically  generated  waves.  In 
Case  I,  the  waves  were  generated  solely  by  the  wind.  Here,  the  dominant 
waves  were  generated  by  mechanical  means.  The  mechanically  generated 
waves  were  then  acted  upon  by  the  turbulent  air  flow.  As  a  result, 
shear  layers  developed  on  both  sides  of  the  interface.  The  water  layer 
was  turbulent  for  any  non-zero  wind  speed  used  in  this  study 
(^1.7  m/s).  Small  wavelets,  hereafter  referred  to  as  ripples,  were 
also  generated  by  the  wind.  The  ripples  randomly  distorted  the  mechan¬ 
ically  generated  wave  profile  from  that  of  an  almost  pure  sinusoid. 

The  mean  velocity  profiles  are  logarithmic  (Fig.  27),  thus  implying 
the  existence  of  the  inertial  sublayer.  Although  the  slopes  of  the 
profiles  are  not  2.5  (■  1/0.4)  they  are  practically  identical  to  those 
in  the  high  wind-speed  wind-wave  experiments.  So  the  significant  veloc¬ 
ity  scale  is  u*  and  the  velocity  defect  distribution  varies  with  -icy. 
Of  course  the  value  of  k  may  not  necessarily  be  equal  to  that  usually 
taken  by  the  von  Karman  constant  (0.4)  because  of  a  different  slope  of 
the  velocity  profiles,  but  the  constant  should  remain  at  the  same  order 


83 


as  0.4  in  the  scaling  length.  At  u„  ■  1.7  and  2.S  m/a,  the  profiles 
deviate  from  the  logarithmic  regions  as  the  interface  is  approached 
(i.e.,  small  -y+) .  The  data  near  this  region  behave  as  if  the  profiles 
were  in  the  viscous  sublayer  but  at  a  higher  -y+  than  the  expected 
-y+  ~  11.  The  explanation  of  this  behavior  is  given  in  Sec.  5.2. 

The  friction  velocity  u*  is  the  characteristic  scale  for  the  RMS  u* 
and  v',  while  the  turbulent  boundary-layer  thickness  6  is  the  parameter 
that  governs  their  decay  with  depth  (see  Figs.  31  and  32).  Once  more, 
we  must  emphasize  the  difference  in  the  definition  of  the  fluctuating 
quantity  n'  between  Case  I  and  Case  II.  In  Case  I,  even  though  the  fre¬ 
quency  of  some  of  the  surface  waves  may  not  be  the  same  as  that  of  the 
dominant  wave,  they  are  still  classified  as  surface  waves.  In  the  pres¬ 
ent  case,  any  non-periodic  components,  or  any  components  which  do  not 
have  a  fixed-phase  relationship  with  the  mechanically  generated  wave  are 
considered  as  fluctuating  components.  For  example,  the  wind-generated 
ripples  are  considered  as  the  fluctuating  component  of  n.  The  ripples, 
of  course,  contribute  to  parts  of  the  RMS  u'  or  v'.  Although  n'  has  a 
"wave-like"  nature,  the  RMS  u*  and  v'  certainly  do  not  scale  with  any 
attributes  of  the  dominant  wave. 

It  is  relatively  straightforward  to  characterize  the  mechanically 
generated  wave  train  because  its  frequency  and  mean  amplitude  are  well- 
defined.  The  proper  length  scale  for  the  wave-induced  water  motion  is 
the  inverse  of  the  dominant  wave  number  1/k.  This  length  scale  is  a 
special  case  of  Eq.  4-5  for  a  monochromatic  wave,  that  is, 


where  6(f)  is  the  Dirac  delta  function,  fD  is  the  frequency  of  the 

wave,  and  its  amplitude.  The  wave  decay  depth  y^  then  becomes  1/k. 

Similarly,  using  the  above  expression  for  S^f)  and  substituting  into 

,71,1/2  .  _ 

Eq.  4-4,  the  velocity  scale  (uqJ  becomes  2irfjjng//2.  The  reason  for 
not  using  the  whole  n  spectrum  in  deriving  the  wave  velocity  scale  has 
already  been  given  in  Sec.  4.2.5  and  is  not  repeated  here. 

The  turbulence  level  is  low  at  a  depth  of  6  or  near  the  outer 
layer.  However,  the  wave  quantities  are  still  significant  in  this 
region  as  shown  in  Figs.  37  and  38.  This  is  quite  different  from  the 
wind-wave  case  where  bo th  the  turbulent  and  wave  quantities  are  small  in 
the  outer  layer. 

For  the  mean  flow,  the  clear  choice  of  the  velocity  scale  is  the 
mean  Eulerian  surface  drift  and  the  length  scale  is  6  where  there  is 
substantial  mean  velocity  gradient.  The  growth  of  the  boundary  layer  is 
related  to  another  length  scale — the  fetch  L. 

5.1.3  Case  III:  Mechanical  Wave  in  Wave-Following  Eulerian  Frame 

The  surface  conditions  in  this  case  were  identical  to  those  of 
fixed-frame  mechanical-wave  experiments  (Case  II).  The  major  difference 
between  Case  II  and  the  present  one  is  in  the  reference  coordinate. 
Here,  all  the  measurements  were  made  in  the  wave-following  frame  but  the 
obtained  data  are  still  Eulerian  in  nature.  Apparently,  the  relevant 
velocity  and  length  scales  in  Case  II  apply  here.  Consequently,  most 
comments  made  in  the  previous  section  also  follow. 

The  logarithmic  trends  of  the  mean  velocity  profiles  persist;  thus, 
u*  is  the  velocity  scale  and  the  velocity  defect  varies  with  -icy  in  the 
inertial  sublayer.  As  compared  to  Case  II,  the  velocity  probe  was  able 


to  reach  very  close  to  the  interface.  At  the  two  lowest  wind  speeds, 


Che  velocity  probe  was  able  Co  sample  at  depths  near  the  viscous 
sublayer  as  Indicated  by  data  points  which  follow  the  u+  ■  -y*+  line. 

Unfortunately,  the  friction  velocity  u*  and  the  turbulent  boundary- 
layer  thickness  6  cannot  be  found  from  the  measured  turbulent  shear 
stress.  However,  the  estimated  u*'s  and  5's  seem  to  collapse  the  turbu¬ 
lent  intensities  (RMS  u'  and  v’)  equally  well.  Therefore,  the  turbulent 
motion  has  a  length  scale  of  the  order  of  6  and  u*  is  of  the  correct 
magnitude  for  the  turbulent  fluctuations. 

The  wave  field  has  a  characteristic  velocity  of  2itfDns/^2  and  a 
length  scale  of  1/k.  The  RMS  u  and  v  collapse  into  narrow  bands  even 
for  | y |  <  | n | ,  where  ordinarily  u^(x,  |y|  <  |n|,  t)  are  not  well- 
defined.  The  wave-following  frame  does  not  seem  to  have  any  major 
effect  on  the  wave-related  velocities. 

Lastly,  the  mean  Eulerian  surface  drift  is  the  velocity  scale  and 
the  estimated  6  is  the  length  scale  for  the  mean  flow.  The  fetch  may 
also  be  a  relevant  length  scale  for  the  mean  flow. 

5.1. A  Summary 

In  Secs.  5.1.1  to  5.1.3,  we  clearly  showed  the  various  important 
length  and  velocity  scales  that  relate  to  different  components  or 
regions  of  the  flow  field.  Many  of  these  characteristic  scales  have 
been  suggested  by  other  investigators,  such  as  Navrotskii  (1967)  and 
Kitalgorod8kii  et  al.  (1982),  but  the  existence  of  these  scales  has 
never  been  experimentally  demonstrated.  Also,  instead  of  finding,  say, 
the  bulk  momentum  transfer  coefficient  at  the  air-water  interface,  we 
focus  our  inquiry  on  the  relevant  scales  of  the  fundamental  principles 
that  govern  the  transfer  process. 
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Some  of  the  characteristic  scales  are  summarized  in  Table  9.  Note 
that  not  all  the  scales  are  listed,  because  others  can  easily  be  derived 
from  those  tabulated.  The  absolute  magnitudes  of  these  scales  sometimes 
are  not  as  important  as  their  relative  magnitudes.  For  example,  in  the 
wind-wave  experiments,  is  always  much  smaller  than  6.  On  the  other 
hand,  6  is  always  less  than  the  corresponding  y^  in  all  mechanical-wave 
experiments . 

The  turbulent  boundary- layer  thickness  6  is  larger  in  the  wind-wave 
case  than  in  the  mechanical-wave  case  at  comparable  wind  speeds.  This 
may  be  due  to  the  suppression  of  the  wind-generated  ripples  by  the 
mechanical  waves  as  suggested  by  Bole  and  Hsu  (1969).  However,  the 
friction  velocities  remain  the  same  order  of  magnitude  at  equivalent 
wind  speeds  regardless  of  surface  conditions.  Of  course,  all  yn's  and 
6's  are  small  compared  to  the  fetch  of  13  m. 

5.2  MEAN  VELOCITY  PROFILES  IN  WALL  COORDINATES 

Some  mean  velocity  profiles  in  Cases  I  and  II,  where  waves  were 
dominant,  have  slopes  different  from  2.5  in  the  logarithmic  region. 
Furthermore,  there  is  a  change  of  slope  in  a  number  of  these  profiles. 
Such  behavior  is  not  found  in  typical  two-dimensional  boundary-layer 
type  flows,  even  under  mild  pressure  gradients.  Tennekes  and  Lumley 
(1972)  pointed  out  that  the  pressure  gradient  can  have  an  effect  on  the 
slope  of  the  velocity  profile,  but  we  will  show  that  the  pressure  gradi¬ 
ent  cannot  be  the  cause  of  the  change  in  slope  or  the  slopes  being  less 
than  2.5  in  some  of  the  mean  profiles.  At  low  wind  speeds,  say  u^  < 
3  m/s  in  both  Case  I  and  II,  the  pressure  gradient  due  to  water  set-up 
or  surface  tilt  was  negligible.  In  Case  I,  the  slopes  of  the  mean  pro¬ 
files  for  u.  <  3  m/s  follow  a  2.5  slope  while  the  profiles  in  Case  II 


have  a  smaller  slope.  The  primary  distinction  between  Cases  l  ami  11  1* 
the  presence  of  the  mechanically  generated  waves  in  Case  II.  Of  course 
in  Case  1,  the  wind  was  not  strong  enough  to  generate  any  waves  of  sig¬ 
nificant  magnitude.  Hence,  by  comparing  the  low  wind-speed  experiments 
in  Cases  I  and  II,  we  conclude  that  the  pressure  gradient  is  not  the 
primary  cause  for  the  anomalous  slopes  of  the  mean  velocity  profiles. 

Another  possible  cause  is  that  the  flow  might  not  be  two- 
dimensional.  The  law  of  the  wall  (Eq.  5-1)  does  not  hold  for  three- 
dimensional  flows  (White,  1972).  The  low  wind-speed  experiments  in 
Case  I  were  nearly  two-dimensional  near  the  channel  center-line  and  the 
mean  profiles  follow  the  expected  trend  of  such  flows.  In  Case  II,  the 
mechanically  generated  waves  were  two-dimensional  and  the  ripples  were 
even  smaller  than  those  in  Case  I  at  low  wind  speeds,  so  the  flow  field 
beneath  the  waves  should  also  be  nearly  two-dimensional.  Although  there 
are  data  scatter  in  the  mean  velocities,  any  unsteadiness  or  net  three- 
dimensional  effect  of  the  flow  cannot  give  such  consistent  trends  in  the 
mean  profiles.  Three-dimensionality  of  the  mean  flow  was  not  present 
and  therefore  cannot  be  the  cause  of  the  usual  slopes  of  the  mean  velo¬ 
city  profiles. 

The  law  of  the  wall  given  by  Eq.  5-1  can  be  rewritten  as: 

u+  ■  —  ln(— for  y  <  0  (5-3) 
k  y^ 

where  the  various  constants  are  absorbed  into  yQ.  In  Cases  I  and  II, 
measurements  of  the  mean  profiles  were  made  at  fixed  depths.  If  the 
mean  flow  with  velocity  distribution  given  by  Eq.  5-3  above  follows  the 
motion  of  the  dominant  waves,  then  Eq.  5-3  should  be  modified  as  fol- 


(5-4) 


■  ^  l ln(h  -  y )  -  In  yQ] 

where  h  ■  h(t)  *  n(t)  ■  surface  elevation.  The  difference  (Au+) 
between  Eq.  5-3  and  Eq.  5-4  gives  the  change  in  u+  due  to  the  logarith¬ 
mic  profile  moving  up  and  down  with  the  dominant  waves.  It  is  clear 
that  the  biggest  difference  should  occur  at  points  very  close  to  the 
surface  and  for  large  |  n |  -  Consider  a  typical  run  in  Case  II  where  |n| 
was  large,  with  the  sign  of  Au+  aside,  the  difference  is  at  most  1.3  u+ 
units.  This  difference  can  account  for  the  apparent  drop  in  u+  at  small 
-y+  ~  50  (see  Fig.  27)  for  low  wind-speed  experiments.  However,  away 
from  the  interface,  Au+  is  considerably  less  than  the  deviation  of  the 
measured  profile  from  the  relation  given  by  Eq.  5-4.  A  more  elaborate 
scheme  of  calculating  Au+  was  suggested  by  Phillips  (1977).  This  scheme 
was  incorporated  into  the  data  analysis  program  and  the  results  are 
tabulated  in  Appendix  D.  The  two  estimates  for  Au+  agree,  but  both  are 
too  small  to  account  for  the  observed  discrepency  at  large  -y+  values- 
If  indeed  the  anomalous  slopes  were  due  to  kinematics  of  the  waves  and 
measurements  being  made  in  a  fixed  frame,  then  the  anomaly  should  vanish 
in  the  wave-following  experiments  (Case  III).  The  mean  velocity  pro¬ 
files  in  Case  III  also  have  slopes  smaller  than  2.5  in  the  logarithmic 
region.  Thus,  the  change  in  slope  or  the  profiles  having  slopes  less 
than  2.5  cannot  be  the  effect  of  measuring  velocities  at  fixed  depths. 

Numerous  possible  explanations  for  the  observed  trends  of  the  mean 
profiles  were  rejected  based  on  the  above  discussions.  A  logical 
deduction  is  that  the  dynamics  of  the  waves  are  responsible  for  the 


unusual  mean  profiles. 


The  normalized  RMS  u'  and  RMS  u  in  Case  I  are 
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plotted  versus  -y+  In  Pigs.  59  and  60,  respectively.  Hence,  it  is  pos¬ 
sible  to  compare  the  change  in  relative  magnitudes  of  RMS  u'  and  RMS  u 
with  -y+.  For  u,,  <  3  m/s,  the  RMS  u'  is  always  larger  than  RMS  u  at 
all  depths.  At  higher  wind  speeds,  there  are  values  of  -y+  at  which 
RMS  u'  >  RMS  u  and  also  for  RMS  u'  <  RMS  u.  If  the  waves  do  affect  the 
mean  profiles,  then  in  the  region  where  RMS  u'  <  RMS  u,  the  mean  pro¬ 
files  may  deviate  from  the  2.5  slopes.  On  the  other  hand,  the  slopes 
should  approach  2.5  in  regions  where  RMS  u'  >  RMS  u.  For  Case  1,  it  is 
found  that  in  the  region  where  RMS  u'  >  2  (RMS  u),  the  mean  profiles 
have  slopes  nearly  equal  to  2.5;  and  where  RMS  u'  <  RMS  u,  the  slopes 
are  smaller  than  2.5.  Furthermore,  the  slopes  change  from  2.5  to  less 
than  2.5  at  depths  where  RMS  u'  a  2  (RMS  u),  such  as  those  for  3.2  m/s  < 
u„  <  9.9  m/s.  These  observations  were  derived  from  experiments  in 
Case  1;  if  these  observations  are  the  outcomes  of  some  basic  physical 
principles,  then  they  should  apply  to  Cases  II  and  III  also.  This  is 
indeed  the  case,  because  generally  RMS  u'  <  RMS  u  in  mechanical -wave 
experiments  and  the  mean  profiles  were  found  to  have  slopes  less  than 
2.5.  In  addition,  there  is  no  change  of  slope  (in  the  logarithmic 
region)  in  the  profiles  for  Cases  II  and  III,  because  there  does  not 
exist  any  region  where  RMS  u*  *  RMS  u-  So  in  the  light  of  the  experi¬ 
mental  evidence,  the  waves  seem  to  have  an  effect  on  the  mean  flow. 


5.3  POSSIBLE  EVIDENCE  OF  WAVE-GENERATED  TURBULENCE 

In  Case  I,  the  turbulent  Intensity  (RMS  u'  or  v’)  profiles  are  very 
similar  to  those  in  a  steady  pure  shear  flow.  That  is,  the  turbulence 
is  non-isotropic  with  RMS  u’  >  RMS  v'.  In  addition,  the  RMS  u  exceeds 
RMS  u',  but  at  most  by  a  factor  of  about  two.  Although  the  wave  fluctu¬ 
ation  is  strong,  it  remains  at  the  same  order  as  the  turbulent  motion. 
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Also,  the  wave  motion  decays  fairly  rapidly  with  depth.  On  the  other 
hand,  the  vertical  and  horizontal  turbulent  Intensities  in  Case  II  are 
almost  equal  in  magnitude,  i.e.,  the  turbulence  is  nearly  isotropic.  The 
RMS  u  now  exceeds  RMS  u'  by  more  than  a  factor  of  two.  This  factor 
attains  a  value  of  10  in  low  wind-speed  experiments  when  waves  were  pre¬ 
dominant.  The  decay  of  the  wave  motion  with  depth  is  not  as  rapid  as 
that  in  Case  I  because  the  mechanical  waves  have  a  longer  wavelength. 
Thus,  if  there  is  wave-generated  turbulence,  the  phenomenon  ought  to  be 
more  prominent  in  mechanical-wave  experiments.  The  apparent  deviation 
of  the  turbulent  intensities  in  Cases  11  and  III  from  those  of  the  wind- 
wave  experiments  may  indicate  such  a  phenomenon. 

If  the  waves  are  truly  irrotational,  the  mean  wave-induced  shear 
stress  (-u  v)  should  be  zero,  and  the  wave  would  have  no  interaction 
with  either  the  mean  or  the  turbulent  flow  fields.  It  is  clear  from  the 
data  that  the  waves  are  not  truly  irrotational,  thus,  the  mean  wave- 
induced  shear  stress  is  non-zero.  This  shear  stress  may  contribute  to 
turbulence  production  directly  or  indirectly.  The  mean  wave-induced 
shear  stress  is  much  larger  in  Cases  II  and  III  than  in  Case  I.  This 
again  implies  that  the  effect  of  the  mean  wave-induced  shear  stress  is 
stronger  in  mechanical-wave  experiments. 

The  experimental  data  support  the  postulate  that  the  waves  can 
affect  the  mean  flow  as  discussed  in  Sec.  5.2.  It  is  well-known  that 
turbulence  can  only  be  maintained  by  drawing  energy  from  the  mean  flow, 
or  as  suggested  above,  from  the  waves  also.  Hence,  the  mean,  wave,  and 
turbulent  fields  are  closely  coupled.  Although  the  experimental  evidence 
for  such  coupling  is  strong,  it  is  crucial  to  demonstrate,  from  a  theo¬ 
retical  viewpoint,  that  the  governing  equations  allow  such  interactions. 
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5.4  ANALYSIS  OF  THE  FIXED-FRAME  GOVERNING  EQUATIONS  PERTAINING 

TO  THE  EXPERIMENTAL  OBSERVATIONS 

The  boundary-layer  equatlona  for  two-dimensional  incompressible 
flows  can  be  derived  from  Che  Navlar-Stokes  equations  by  an  order-of- 
magnitude  analysis.  In  many  flow  situations,  such  as  those  over  flat 
plates,  the  characteristic  scales  are  well-known.  Therefore,  oue  can 
proceed  directly  to  the  boundary- layer  equations  without  rederiving 
them.  Here,  some  of  the  characteristic  scales  and  their  relative  magni¬ 
tudes  are  not  known  a  priori.  Hence,  an  order-of-magnitude  analysis  of 
the  governing  equations  is  impossible  unless  the  relevant  scales  are 
found  experimentally  or  otherwise.  The  characteristic  scales  for  dif¬ 
ferent  parts  of  the  flow  are  found  experimentally  in  Sec.  5.1;  thus,  an 
order-of-magnitude  analysis  can  now  be  performed  on  the  governing  equa¬ 
tions. 

5.4.1  Primary  Estimates 

To  do  an  order-of-magnitude  analysis  of  the  governing  equations, 
variables  will  be  replaced  by  their  "typical"  magnitudes.  The  changes 
in  the  variables  are  taken  to  be  of  the  same  order  as  the  typical 
values.  The  symbol  "A"  that  preceeds  a  variable  denotes  changes  in  that 
variable.  The  notation  denotes  a  typical  value  of  order-of- 

magnitude  accuracy  regardless  of  the  sign.  This  should  not  be  confused 
with  the  superposed  which  stands  for  wave  quantities.  In  the  fixed- 
frame  equations,  the  following  typical  values  are  used: 

Ax  -  L  for  mean  quantities 

-  cAt  **■  c/w  ~  y^  for  wave  quantities 

At  **•  1/ui 

Ay  ~  6  for  mean  quantities 

~  yn  ~  1/k  for  wave  quantities 
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yn  <  6  in  Case  I 
y^  —  6  in  Case  II 

U  ~  Ug  ~  Ug 

u,  Au  **•  aw 

v,  Av  ~  aw 

where 

L  3  length  scale  in  the  streamwise  direction 
w  *  circular  frequency  of  the  dominant  wave 
6  ■  turbulent  boundary-layer  thickness  (6  «  L) 
yn  ■  wave-decay  depth  (y^  «  L) 
us  =  mean  Eulerian  surface  drift  velocity 
a  -  amplitude  of  the  dominant  wave 
k  •  dominant  wave  number 
c  *  dominant  wave  speed  (■  w/k). 


It  is  assumed  that  the  mean  turbulent  quantities 

A(u'u' ) ,  A(u' v ' ) ,  A(vrv' )  ~  e^Ug  (5-5) 

where  e2  is  a  multiplicative  factor  (usually  «  1)  whose  magnitude  will 

be  determined  in  the  analysis.  This  order-of-magnitude  estimate  is 

arrived  at  by  the  following  reasoning.  The  mean  turbulent  quantities, 
2  2 

such  as  u'  ,  are  of  the  order  of  u*,  and  u*  is  found  to  be  roughly  pro¬ 
portional  to  Ug,  so  that  the  mean  turbulent  quantities  are  of  the  order 
of  e2u2. 

The  order-of-magnitude  changes  in  u'u',  u'v',  and  v'v*  are  derived 
as  follows.  The  phase-averaged  value  <u'u’>  is  assumed  to  be  given  by: 

<u'u’>  ~  e2  <u^>2 
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where  <ug>  is  the  phase-averaged  Euler ian  surface  drift  velocity.  This 
estimate  is  a  generalization  of  Eq.  5-5,  which  applies  to  the  mean  quan¬ 
tities.  Now,  if  <Ug>  is  given  by: 


<Ug>  ~  Ug  +  au>  cog  uit 


then 


2  2  2  2  2 
<Ug>  ~  ug  +  2ugaui  cos  wt  +  a'  w  cos 


ut 


■  Ug[l  +  2(~)  cos  ut  +  (^)  cos^  ut] 

s  s 


2  2  2 
•  u|[l  +  2a  cos  ut  +  a  cos  ut] 


where  a  ■  au/ug. 

From  the  experimental  results,  a  ]>  1  in  Cases  II  and  III,  so  u'u', 
u'v',  v'v*  ~  e^u|.  It  is  unnecessary  to  consider  Case  I  because  the 
wave-induced  turbulent  stresses  are  not  defined.  From  the  above  analy¬ 
sis,  u'u'  is  of  the  same  order  as  u'u',  and  similarly  for  the  other  com¬ 
ponents.  This  result  is  quite  different  from  those  of  Carr  ( 1981a ,b) 
and  Jayaraman  et  al.  (1982)  for  the  unsteady  boundary  layer.  They  found 
that  u'u'  is  much  smaller  than  u'u'. 

5.4.2  Continuity  Equations 

The  mean  continuity  equation  (Eq.  2-17a)  has  two  terms  and  their 
or der-of -magnitude  estimates  are  given  by: 


3u  US 
3x  ~  L 


.  3v  Av 

■nd  Jy  ~r 


I. 

L 


Each  term  must  be  of  the  same  order  for  the  equation  to  hold,  so  we  must 
have: 

-  6 


it  -  ■  -  ■ 
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that  is,  the  change  in  the  mean  vertical  velocity  is  very  small.  This 
agrees  with  the  experimental  observations  that  the  mean  vertical  velo¬ 
city  gradient  is  very  small  and  the  mean  vertical  velocity  is  very  close 
to  zero. 

Care  must  be  taken  in  analyzing  the  wave  continuity  equation 
(Eq.  2-17b)  in  order  to  avoid  some  pitfalls.  The  wave-induced  or 
periodic  velocities  are  known  to  be  of  the  same  order  of  magnitude,  and 
the  two  terms  in  the  equation  must  also  be  of  the  same  order.  Thus,  the 
length  scales  Ax  and  Ay  must  be  of  the  same  size.  However,  Ax  cannot  be 
the  same  as  the  primary  estimate  for  the  mean  variables  L;  otherwise  Ay 
would  be  L,  which  is  impossible.  Instead,  we  know  that 

3v  au 

ay  ~  yn 

u  ~  au 

k  ~  1/y^  and  c  *  u/k 


Since  we  require  a  balance  between  the  x  and  y  derivatives, 


and  because 


3u  3v  au  au 

M  '■  II  II  —  — 

3x  3y  yr,  c 


3u  __  au 
3x  ~  Ax  ’ 


it  follows  that  the  required  x  scale  is  Ax  ~  c/u  *  cAt;  this  gives  the 
typical  Ax  value  for  wave  quantities. 


5.4.3  Mean  Momentum  Equations 


There  are  nine  terms  in  the  mean  x-momentum  equation  (Eq.  2-18). 


The  order-of-magnitude  estimate  of  each  term  is  as  follows: 


Term  1 


Term  2 


Term  3 


Term  4 


Term  5 


Term  6 


Term  7 


Term  8 


Term  9 


(u  v  )  - 


i2  2  2 

S  a  in 


i  2 1 
r° 


Term  4  Is  negligible  compared  to  Term  5,  and  Term  6  is  negligible  com¬ 
pared  to  Term  7,  because  6  «  L.  As  yn  «  L,  Term  9  is  much  bigger 
than  Term  8.  The  estimate  in  Term  8  is  for  Case  I,  because  when  the 
waves  are  generated  by  the  wind,  the  wave  orbital  velocity  has  acquired 
a  magnitude  of  au)  at  the  fetch  L.  This  term  should  even  be  smaller 
compared  to  Term  9  for  Case  II,  because  there  is  little  growth  of  the 
mechanically  generated  waves  over  the  distance  L  for  the  wind  speeds 
used  in  the  experiments.  The  viscous  term  (Term  5)  is  important  when 


u  L  -1/2 

(-5-) 

v  v  ) 


or  in  regions  very  close  to  the  Interface.  The  turbulent  term  (Term  7) 
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Is  significant  when 


2  6 
£  ""  L 

2 

and  this  yields  the  order-of-magnitude  estimate  of  c  •  The  parameter  a 

is  always  on  the  order  of  one  or  larger,  except  for  low  wind-speed 

experiments  in  Case  I,  so  Term  9  is  always  significant  when  waves  are 

present-  Finally,  for  the  pressure  term  to  be  of  any  importance,  we  need 

to  have  _ 

—  LuS  2 

Ap  -  p  —  -  pus 

After  discarding  the  insignificant  terms,  the  mean  x-momentum  equa¬ 
tion  simplifies  to: 

_  2-  _ _ 

3-  (u  u)  +  _  (U  v)  -  -  -  +  v  ^  ^  (u  v  +  u  v  )  (5-6) 

This  equation  is  similar  to  the  steady  boundary  layer  x-momentum  equa¬ 

tion  except  for  the  extra  mean  wave-induced  shear  stress  term  -u  v.  The 
mean  wave-induced  shear  stress  and  its  vertical  gradient  are  not  small 
near  the  interface  when  waves  are  present,  especially  in  mechanical-wave 
experiments.  Hence,  the  waves,  which  lead  to  -u  v,  can  affect  the  mean 
flow.  In  Sec.  5.2,  it  was  pointed  out  that  the  mean  velocity  profiles 
in  wall  coordinates  change  slope  at  depths  where  RMS  u'  =  2(RMS  u).  The 
mean  wave-induced  shear  stress  is  related  to  the  wave-induced  velocity 
u,  which  of  course  is  represented  by  RMS  u.  On  close  examination  of  the 
data  (e.g.  Appendix  D),  the  depth  at  which  RMS  u'  =  2(RMS  u)  corresponds 
to  the  distance  over  which  there  is  substantial  change  in  the  mean  wave- 

induced  shear  stress.  In  Case  I,  the  low  wind-speed  experiments  have 

negligible  -u  v,  so  Eq.  5-6  is  identical  to  the  mean  x-raoraentum  equation 
for  a  steady  boundary  layer.  However,  this  asymptotic  result  does  not 


hold  for  low  wind-speed  experiments  In  Case  II,  because  -u  v  Is  Mil  11 


significant  due  to  the  mechanically  generated  waves.  It  is  clear  that 
the  experimental  observations,  such  as  those  in  Sec.  5.2,  agree  with  the 
implications  of  the  simplified  mean  x-momentum  equation  (Eq.  5-6). 

The  mean  y-momentum  equation  (Eq.  2-19)  also  has  nine  terms  as  well 
as  the  body-force  term.  Their  order-of -magnitude  estimates  are  as  fol- 
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It  is  readily  seen  that  Term  4  is  negligible  compared  to  Term  5,  Term  6 
is  much  smaller  than  Term  7,  and  Term  9  is  usually  much  larger  than 
Term  8.  So  in  this  first  step  of  simplification,  we  have: 


k  <*  u)  +  Ty  v) 


_  I  3P  + 


32v 


P#+v^-W(v,v,)-Fy(vv)-8 


2 

If  the  order-of-magnitude  estimate  for  e  ~  6/L  is  used  in  the  analy¬ 
sis,  then  Terms  1,  2,  and  5  are  also  negligible,  because  UgL/v  is  gener¬ 
ally  large.  Thus,  the  equation  is  further  simplified  to: 


1  3p 

p  "5y 


-g  -  ^  (v'v'  +  V  V  ) 


(5-6) 


In  experiments  where  waves  are  insignificant,  we  have  in  the  limit: 


v  v  ~  0  and  (v  v)  ~  0 


An  order-of-magnitude  estimate  for  the  pressure  term  (Term  3)  now 
becomes:  - 

V 

Ap  -  pg6  +  p  -j— 

~  Pg6 

That  is,  Term  7  is  small  compared  to  the  body  force.  So  in  this  special 
case,  the  pressure  term  is  just  equal  to  the  hydrostatic  pressure.  When 
waves  are  present,  we  have: 

u26  2  2 

Ap  ~  pg5  +  p  -j—  +  p  — — 

2  2 

f  ,  a  u)  . 

~  Pg6  +  P  — —  6 

yn 

2  2 

In  almost  all  experiments  a  w  /y  <<  g,  so  we  arrive  again  at: 

Ap  ~  pg6 
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So  the  hydrostatic  assumption  of  the  pressure  term  holds  in  all  cases, 
and  Eq.  5-6  simply  becomes: 


IE  .  . 


5.4.4  Wave-Perturbed  or  Periodic  Momentum  Equations 

The  periodic  momentum  equation  in  the  x  direction  (Eq.  2-20)  has  11 
terms.  Each  of  these  terms  has  the  following  order-of-magnitude  esti¬ 
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There  are  two  estimates  for  each  term  of  Terms  2  through  4  because  each 
of  them  has  two  characteristic  length  scales.  The  larger  of  the  two 
estimates  is  used  in  comparing  the  terms.  Terms  9  and  11  contain  turbu¬ 
lent  quantities,  but  because  the  phase  average  is  used,  the  wave  length 
scale  is  used  instead  of  the  turbulent  length  scale.  The  predominant 
term  clearly  is  Term  1.  Other  terms,  except  Term  5  are  roughly  of  the 
same  size.  Thus,  none  of  them  can  be  neglected  without  further  assump- 

tlons.  If  the  estimate  e  ~  6/L  is  used,  then  Terms  9  and  11  are 

2  o 

small  compared  to  Terms  8  and  10.  Furthermore,  if  ( v/ugL)(L//y^)  £  1, 
then  the  viscous  terms  (6  and  7)  are  negligible.  Hence,  we  have  the 
following  simplified  equation: 


Tt  +  2  k  (G  u)  +  h  (u  "  +  G  v) 


1  3p  3  3  -v 

p  Tx  '  -E  (u  u)  *  Fy  (u  v) 


Comparing  Term  1  with  the  remaining  terms,  we  see  that  only  the  pressure 
term  is  large  enough  to  balance  the  equation.  Furthermore,  the  wave- 
induced  waves  stresses  ~u^Uj  have  a  dominant  frequency  twice  that  of  u 
as  pointed  out  in  Sec.  4.2.7;  the  interaction  of  these  terms  with  Term  1 
must  be  nonlinear  and  is  generally  small.  Thus,  we  end  up  with: 


3u  _  _  1  3p 
3t  p  3x 


which  is  just  the  equation  of  motion  in  the  x  direction  for  small 
amplitude  gravity  waves  on  an  ideal  fluid  (Landau  and  Lifshitz,  1959). 
This  simple  analysis  shows  that  the  mean  and  the  turbulent  fields  have 
little  effect  on  the  wave  field.  However,  this  does  not  imply  the 
reverse,  that  is,  the  waves  may  have  some  significant  impact  on  the  mean 
and  turbulent  fields  due  to  the  nonlinear  terms. 
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The  periodic  y-moraentum  equation  (Eq.  2-21)  also  has  11  terms. 
They  are  listed  with  their  order-pf-magnitude  estimates  as  follows: 
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Similarly  to  the  periodic  x-momentum  equation,  Term  1  and  Term  5  are  the 
dominant  terms.  If  the  other  terms  are  neglected,  the  equation  of 
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motion  becomes  the  periodic  y-momentum  equation  for  small-amplitude 
gravity  waves  of  an  irrotational  fluid: 


3v  1  3p 

3t  p  3y 


The  body  force  due  to  gravity  is  not  included  in  the  equation  because  it 
is  accounted  for  in  the  mean  equation: 


1  0£  _  _ 
p  8y 
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Let  us  compare  the  above  analysis  with  the  experimental  data.  The 
RMS  wave-induced  velocity  was  computed  at  each  depth  according  to  the 
linear  water-wave  theory.  These  values  are  listed  in  Appendix  D  with 
the  measured  RMS  u  and  v.  The  RMS  v  are  within  5%  of  the  theoretical 
values.  The  RMS  u  are  always  less  than  the  theoretical  values,  and  less 
than  the  measured  RMS  v  also .  One  may  argue  that  in  the  wind-wave 
experiments,  the  waves  may  not  propagate  strictly  in  the  downstream 
direction  but  at  some  angle  to  it.  Hence,  only  a  portion  of  the  RMS  u 
was  measured.  However,  similar  results  are  found  in  Case  II  where  the 
mechanical  waves  are  two-dimensional.  One  explanation  is  that  Terms  2 
and  3  in  the  horizontal  wave  momentum  equation  are  bigger  than  Terms  2 
and  3  in  the  corresponding  equation  in  the  vertical  direction  (see 
Sec.  3.4.4).  That  is,  the  interaction  of  the  waves  with  the  mean  flow 
is  stronger  in  the  horizontal  direction  than  in  the  vertical  direction. 
Thus,  more  energy  is  "lost"  by  the  waves  in  the  horizontal  direction 
which  results  in  a  smaller  RMS  u  than  the  theoretical  value. 


5.4.5  Energy  Equations 

The  momentum  equations  presented  in  Sec.  5.4.4  show  that  coupling 
exists  among  the  mean,  wave,  and  turbulent  fields;  but  the  interaction 
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where 


D_  3-3 
Dt  *  3t  Uj  3x 

J 

is  the  derivative  following  a  fluid  particle  in  the  mean  flow. 

The  energy  equation  for  the  mean  flow  (Eq.  5-7)  does  not  contain 
any  more  information  than  the  momentum  equation  for  the  mean  flow 
because  the  former  is  obtained  from  the  latter  by  algebraic  manipula¬ 
tion.  However,  additional  insight  into  the  dynamics  of  the  flow  can  be 
gained  fronv__the  energy  equation.  The  terms  on  the  right-hand  side  of 
Eq.  5-8  that  have  the  form  3/3xj[  ]  represent  transport  or  redistribu¬ 
tion  of  energy  within  some  volume.  The  work  done  by  the  body  force  is 
represented  by  -F^u^.  The  last  terra  in  the  equation  stands  for  viscous 
dissipation.  The  terms  of  primary  interest  are  those  defined  as  I  and 
II.  Term  I,  when  it  is  positive,  represents  the  production  of  the  wave 
kinetic  energy  by  the  mean  wave-induced  stresses  -u^uj.  The  same  term 
appears  with  an  opposite  sign  in  the  equation  for  wave  kinetic  energy 
(Eq.  5-9).  Similarly,  Terra  II  represents  the  production  of  the  turbu¬ 
lent  kinetic  energy  by  the  mean  Reynolds  stresses  -u  ju  ’ .  This  term 
shows  up  as  an  energy  source  terra  in  the  equation  for  turbulent  kinetic 
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energy  (Eq.  5-10).  The  term  defined  as  III  In  Eqs.  5-9  and  5-10  repre¬ 
sents  the  turbulent  energy  production  by  the  waves  via  the  action  of  the 


phase-averaged  turbulent  Reynolds  stresses  -<ujuj>  ■  -u’Uj  -  r^j. 
This  term  can  also  be  rewritten  as  -r^ j(3u^/ 3x^) .  It  is  present  in  the 
energy  equations  for  both  the  wave  and  turbulence,  but  with  opposite 
signs . 

Terms  I,  II,  and  III  clearly  show  the  interaction  among  the  mean, 
wave,  and  turbulent  fields.  Term  I  denotes  a  summation  of  9  terms, 
among  which  -u  v  (3u/3y)  is  dominant.  To  find  the  mean  velocity  gradi¬ 
ent,  polynomials  were  least-square  fitted  to  the  mean  u  velocity  profiles 
and  their  gradients  were  computed.  Next,  the  product  -u  v  (3u/3y)  was 
calculated  (listed  in  Appendix  D).  The  mean  velocity  gradient  3u/3y  is 
positive  in  Cases  I  and  II  (for  u*,  j4  0),  and  -u  v  is  generally  negative. 
Thus,  the  wave  production  term  is  negative;  in  other  words,  we  have 
energy  transfer  from  the  wave  field  to  the  mean  field.  Hussain  (1983) 
suggested  that  negative  production  is  a  simple  consequence  of  coherent 
structures  in  turbulent  shear  flows.  Here,  whether  the  surface  waves 
can  be  classified  as  coherent  structures  is  unimportant;  we  merely  want 
to  point  out  that  negative  production  from  organized  flow  structures  is 
not  an  uncommon  phenomenon.  The  production  term  was  found  to  be  more 
negative  near  the  interface  and  decays  (less  negative)  rapidly  with 
depth  as  both  -u  v  and  3u/3y  diminish  away  from  the  interface.  Some 
typical  profiles  of  this  production  term  are  shown  in  Fig.  61  for  Case  I 
and  in  Fig.  62  for  Case  II. 

The  negative  production  of  the  waves  means  that  they  Increase  the 
kinetic  energy  of  the  mean  flow.  If  energy  is  drained  from  the  waves  to 
the  mean  flow,  then  how  do  the  waves  sustain  themselves?  Indeed,  in 
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some  flow  situations,  the  waves  grow  with  fetch.  Thus  far,  we  have  neg¬ 
lected  energy  transfer  from  the  wind  at  the  interface.  It  is  obvious 
that  there  exists  direct  energy  transfer  from  the  wind  to  the  mean  and 
wave  fields  in  the  water  across  the  interface.  The  energy  production  or 
dissipation  terms  in  the  equations  just  indicate  the  direction  of  trans¬ 
fer  within  the  water  layer.  Therefore,  the  mean  and  wave  fields  in  the 
water  sustain  themselves  by  drawing  energy  from  the  wind. 

Turbulence  seems  to  draw  its  energy  from  the  mean  field  through 
Term  II  or  from  the  wave  field  through  Term  III.  Note  that  the  phase 
average  is  not  defined  in  Case  I;  however,  the  phase  average  used  in  the 
wave  and  turbulent  kinetic  energy  equations  may  be  interpreted  as  some 
generalized  phase  average  based  on  the  dominant  wave.  Although  the  gen¬ 
eralized  phase  average  is  not  computed,  it  still  allows  us  to  gain  some 
insight  from  the  energy  equations. 

The  dominant  component  in  Term  II  is  -uT\Tr(  3u/3y) ,  and  it  is  always 
positive  in  the  experiments.  Thus,  energy  is  drained  from  the  mean 
field  to  the  turbulent  field.  The  dominant  turbulence  production  term 
was  computed  and  summarized  in  Appendix  D  for  Cases  I  and  II.  Some 
representative  profiles  for  the  dominant  turbulence  production  term  are 
shown  in  Figs.  63  and  64.  In  wind-wave  experiments,  the  magnitude  of 
the  turbulence  production  Term  II  is  higher  than  the  magnitude  of  the 
wave  production  Term  I  at  all  wind  speeds.  However,  both  terms  are 
about  the  same  order  of  magnitude  near  the  Interface  in  some  cases.  On 
the  other  hand,  the  wave-production  term  is  higher  in  magnitude  than  the 
turbulence-production  term  in  mechanical-wave  experiments.  Thus,  the 
sum  of  Terms  1  and  II  is  positive  In  the  mean  kinetic  energy  equation 
(Eq.  3-7).  That  is,  the  net  effect  of  the  two  production  terms  is 
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similar  to  the  turbulence  production  term  in  a  pure  shear  flow — a  net 
drain  of  energy  from  the  mean  flow.  This  may  be  the  reason  why  the 
turbulent  intensities  (RMS  u'  and  v')  have  trends  similar  to  those  for 
flows  over  flat  plates.  However,  the  net  production  is  negative  in 
mechanical -wave  experiments,  which  may  explain  the  unusual  trends  of  the 
turbulent  intensities  for  mechanical-wave  experiments. 

There  are  four  dominant  components  in  Term  III,  namely, 


3u 


3u  -  3v  -  3v 

“rll  3x*  “r12  3y *  ~r12  3x  *  and  _r22  3y 


If  we  assume  the  following 


uj_  *  cos(kx  -  tut  -  0  ) 

3uj 


JT  kui 


(5-11) 

(5-12) 


r  ■  r  cos(kx  -  cut  -  0_  ) 

ij  ij  fij 
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then  the  components  in  Term  III  can  easily  be  found  from  the  data  sum¬ 
mary  in  Appendix  D.  However,  in  fixed-frame  experiments,  there  is  much 


uncertainty  in  8j,^.  The  uncertainty  is  not  due  to  instrumentation  but 


to  kinematics.  All  6*  seem  to  have  a  value  about  180°  relative  to  n- 


The  rjj  are  expected  to  be  higher  near  the  interface.  Since  the  mea¬ 
surements  were  made  in  a  fixed-frame,  the  probe  was  closest  to  the 
interface  at  the  wave  trough,  which  is  the  180°  phase  angle  point  of  tv 
So  the  phase  angles  of  r^j  may  be  due  to  kinematics  of  the  waves  and  not 
due  to  dynamics.  The  effects  of  the  wavy  motion  may  not  be  as  severe  at 


deeper  depths  and  may  indeed  be  the  true  phase  angles,  but  Term  III 


at  such  depths  (~  1/k)  is  most  likely  negligible. 
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The  last  term  in  Eqs.  5-8  to  5-10  denotes  viscous  dissipation.  The 
dissipative  time  scale  is  usually  much  smaller  than  the  time  scale  of 
the  waves  and  the  large  eddies.  Thus,  the  viscous  dissipation  term  in 
the  mean  or  wave  kinetic  energy  equation  is  small.  Most  of  the  kinetic 
energy  is  dissipated  by  the  viscous  dissipation  term  in  the  turbulent 
kinetic  energy  equation. 

5.4.6  Summary 

In  Sec.  5.4,  the  characteristic  scales  obtained  in  the  experiments 
were  used  to  simplify  the  governing  equations.  Not  only  do  the  simpli¬ 
fied  equations  give  further  theoretical  insight  to  the  present  research, 
but  also  they  bear  out  many  of  the  experimental  observations.  The  mean, 
wave,  and  turbulent  flow  fields  are  shown  to  be  closely  coupled  in  the 
energy  equations.  Some  of  the  production  terms  are  estimated  using  the 
experimental  data  and  the  directions  of  transfers  are  identified.  It  is 
not  the  aim  of  this  study  to  solve  these  equations,  but  they  do  lend 
themselves  for  future  research  and  analysis. 


5.5  IMPLICATIONS  OF  THE  WAVE-FOLLOWING  FRAME  RESULTS 

The  wave-following  coordinate  system  is  the  most  logical  reference 
frame  for  the  mechanical-wave  experiments  from  either  an  experimental  or 
theoretical  viewpoint.  This  system  allows  velocity  measurements  above 
the  wave  trough  and  clearly  defines  the  interfacial  boundary  for  the 
governing  equations.  However,  this  reference  frame  has  some  minor 
shortcomings.  An  order-of-magnitude  analysis  was  performed  on  the 
fixed-frame  governing  equations  in  Sec.  5.4.  A  similar  analysis  can 
certainly  be  done  on  the  wave-following  frame  equations.  Unfortunately, 
such  a  task  is  quite  formidable,  because  the  coordinate  transform  intro¬ 
duces  many  extra  inhomogeneous  terms  into  the  equations.  A  complete 
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analysis  Is  beyond  Che  scope  of  this  study.  Instead,  analytical  rela¬ 
tionships  are  developed  In  Sec.  5.5.1  which  relate  the  measured  quanti¬ 
ties  in  (x,y,z,t)  and  (x*,y*,z*,t*).  Not  only  do  these  relationships 
allow  us  to  compare  some  of  the  measured  quantities  in  the  two  coordi¬ 
nate  systems,  but  they  also  demonstrate  the  significance  of  utilizing 
the  wave-following  frame  as  is  shown  in  Sec.  5.5.2. 

5.5.1  Relationships  Between  Measured  Quantities  in  (x,y,z,t)  and 
(x*,y*,z*, t*)  Coordinates  and  Their  Implications 
The  phase  average  is  well-defined  in  both  (x,y,z,t)  and  (x*,y*,z*,t*) 
coordinates  provided  -y  >  a.  If  x  ■  x*  and  z  “  z*,  as  in  the  exper¬ 
iments  (x,y  *  h  +  fn)  and  (x*,y*  ■  h)  with  -h  >  a  represent  the  same 
spatial  point.  Then,  as  shown  by  Hsu  et  al.  (1977),  the  phase  averages 
in  the  two  reference  frames  are  related  by: 

<g>(x,y,t)|  -  -  <g>(x*,y*,t*)|  (5-14) 

y*h+fn  y-h 

By  use  of  a  Taylor  series  expansion,  <g>(x,y,t)|  can  be  written  as 

y*h 

<g>(x,y,t)| 

y”h 

■  <g>(x,h  +  fq  -  fn,t) 

-  <g>(x,h  +  fn,t)  -  [|y  <g>(x,h  +  fn,t)]  fn 
2  2-2 

+  [i-y  <g>  (X,h  +  fn,t)J  i-S-  +  o(n3)  (5-15) 

3y 

The  above  expression  is  valid  for  all  values  of  fn  provided  <g>  is 
continuous  together  with  its  higher  derivatives  (say  up  to  second  or 
third  order)  on  an  interval  containing  h  and  h  +  fn.  A  rigorous  proof 
is  not  given  here  to  show  that  <g>  satisfies  these  conditions;  instead  a 
heuristic  argument  is  given  below.  Since  <g>  is  obtained  by  phase 
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averaging,  It  la  periodic  and  usually  sufficiently  smooth  to  satisfy  the 
above  conditions.  For  example,  u  and  v  are  almost  sinusoidal,  so  their 
derivatives  exist  almost  everywhere-  In  addition,  owing  to  their  sine¬ 
like  behavior,  convergence  of  the  series  expansion  is  implied.  Then, 
using  Eq.  5-14  and  the  chain  rules  in  evaluating  the  derivatives, 
Eq.  5-15  become 


<g>(x,y,t)| 


ysh 

<g>(x*,y*,t*) 


3<g> 

y*„h  fn  3y* 


y*"h 


.  -2  3  ,C2  3<«>'|  .  .-3 

+  11  TF  (-J  +  0<"  > 

Next,  we  decompose  <g>  into  g  and  g  and  obtain,  for  -h  >  a, 

^2  3 


(5-16) 


f2  it 


8<x’y)1y-h  *  U(x**y*>  “  **  ^  ly*-  (S-  JJr)]y*.h  +  0(^3)  (5'17) 

g(x,y,t)|  U  -  g(x*,y*,t*)|  -  fn  S*|  +  0(n2)  (5-18) 
y-h  y*-h  ly*»h 

These  two  equations  show  the  relationships  for  a  measured  quantity  in 
the  two  coordinate  systems  at  the  same  mean  depth.  The  difference 
between  g(x,y)  and  g(x*,y*)  is  of  second  order  in  the  wave  quantity. 
The  difference  in  g(x,y,t)  and  g(x*,y*,t*)  is  of  first  order.  Although 
these  equations  are  similar  to  those  of  Hsu  et  al.  (1977),  the  wave 
quantities  are  much  larger  in  this  study  and  the  equations  have  differ¬ 
ent  implications  here. 

If  g  represents  the  horizontal  velocity  u,  then  even  though  the 
difference  between  the  mean  velocity  is  of  second  order  in  the  wave 
quantity,  the  higher-order  term  may  not  be  small  because  u  is  of  the 
same  order  as  u*  The  magnitude  of  the  second  term  in  Eq.  5-17  is 
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(5-19) 


f  n  ~  a  —  ~  ka(au) 
iy*  yn 

Using  the  linear  water-wave  theory,  one  finds  that  this  term  Is  about 
5  mm/ s  near  the  surface  for  the  present  experiments*  The  measured  mean 
horizontal  velocity  Is  higher  in  the  wave-following  frame  experiments 
than  in  the  fixed-frame  experiments  at  equivalent  mean  depths  and  wind 
speeds.  The  difference  is  attributed  to  the  higher-order  terms  in 
Eq.  5-17,  especially  the  second  term  (Eq.  5-19).  They  are  on  the  order 
of  5-10  mm/s,  which  agree  with  the  observed  differences.  Away  from  the 
interface,  the  second-order  terms  are  negligible  due  to  the  decay  rf  the 
wave  and  the  function  f  with  depth.  Although,  the  higher-order  terms  in 
Eq.  5-17  may  be  large  compared  to  the  mean  velocity,  they  represent  only 
a  small  percentage  of  the  mean  velocity  defect  u<*  -  u«  So  the  moan 
velocity  defect  in  the  wave-following  frame  can  be  interpreted  as  the 
same  quantity  in  the  fixed  frame.  Similarly,  if  g  represents  the  verti¬ 
cal  velocity  v,  the  second  term  in  Eq.  5-17  is  small  even  though  v  is 
big,  because  n  and  v  are  almost  90°  out  of  phase;  the  third  term  is  also 
small  because  v  ~  constant  and  3v/3y  ~  0.  Thus,  the  mean  vertical 
velocity  can  be  taken  to  be  the  same  in  both  coordinates  also. 

The  mean  square  wave  quantities  g  in  the  two  coordinates  are  re¬ 
lated  by: 

?(**y>ly.h  ■  [?(**. y*)  “  2£n  If*  +  ]y*„h  +  0(n3)  (5_2°) 

If  g  represents  u,  we  have  the  following  order-of-magnitude  estimate  for 
Eq.  5-20: 

T?  2  2 
Term  1  u  -  a  m 
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Term  2 


Term  3 


un  ip  ~  i  ’  aa’l~T-> 

f2?(|^F)  -  because  f 


1  near  the  interface 


Term  2  is  small  because  a  <  6,  but  cannot  be  ignored  entirely,  espe¬ 
cially  near  the  interface  where  the  mean  velocity  gradient  3u/3y*  is  not 
negligible.  However,  in  the  region  -y*  ^  a,  3u/3y*  is  not  changing 
rapidly,  so  Term  2  remains  very  small.  Term  3  is  very  much  smaller  than 
Term  1  and  can  be  neglected.  When  g  represents  v,  Term  2  is  even 
smaller  than  that  for  u,  because  n  and  v  are  almost  90°  out  of  phase  and 
the  mean  velocity  gradient  3v/3y*  is  practically  zero.  Hence,  the  RMS  u 
or  v  can  be  considered  the  same  in  both  coordinates. 

The  above  discussion  explains  why  many  of  the  wave-following 
results,  such  as  RMS  u  and  the  velocity  defect  ug  -  u,  appear  the  same 
in  both  coordinates  for  -y  >  a.  Unfortunately,  it  is  not  possible  to 
establish  a  relationship  for  the  turbulent  component  g'  in  the  two  coor¬ 
dinate  systems,  because  higher-order  derivatives  may  not  exist  for  g'. 

There  is  much  difficulty  in  relating  the  r^j  measured  in  the  two 
reference  frames.  First,  these  quantities  are  generally  of  small  magni¬ 
tude  and  have  large  uncertainties,  so  the  series  expansion  technique  may 
not  be  valid.  Second,  the  measured  phase  angles  0-  in  the  fixed-frame 

iJ 

experiments  may  not  represent  the  actual  phase  angles  due  to  the  kine¬ 
matics  of  the  waves.  Third,  despite  the  data  scatter,  the  r^j  do  not 
exhibit  any  clear  trend  in  the  wave-following  frame.  Hence,  we  do  not 
attempt  to  develop  any  relationships  to  connect  the  r^  in  different 
frames. 


We  have  shown  that  many  measured  quantities  of  interest  are  coor¬ 
dinate  independent.  Thus,  we  can  compare  these  quantities  without 


resorting  to  complicated  equations,  but  the  equivalence  of  coordinate 
frames  for  -y  >  a  need  not  be  valid  for  higher-order  turbulent  quan¬ 
tities.  As  a  cautionary  note,  th  equivalence  should  not  be  taken  as 
a  proof  that  the  fixed-frame  equations  are  applicable  to  the  wave¬ 
following  frame. 

5.5.2  Significance  of  the  Wave-Following  Frame  Movements 

The  phase  average  <g>  clearly  is  not  defined  for  -y  <  a.  By  use 
of  a  traditional  perturbation  scheme,  g(x,y)  and  g(x,y,t)  at  -y  <  a 
can  be  calculated  from  g(x*,y*)  and  g(x*,y*,t*)  according  to  Eqs.  5-17 
and  5-18.  However,  interpretation  of  the  wave-following  measurements  in 
such  a  manner  would  seem  to  defeat  the  purpose  of  using  the  wave¬ 
following  coordinates.  Here,  we  want  to  point  out  what  we  have  gained 
in  using  the  transformed  coordinates  and  suggest  possible  improvements 
in  future  research. 

In  mechanical-wave  fixed-frame  experiments  (Case  II),  the  mean 
velocity  defect  profiles  in  wall  coordinates  clearly  illustrate  the 
shortcomings  of  the  fixed-frame  measurements.  The  decrease  of  u+  at  -y+ 
around  50,  as  pointed  out  in  Sec.  4.2.3  and  Sec.  5.2,  is  due  to  averag¬ 
ing  velocities  at  points  close  to  the  wave  troughs  and  points  some 
distance  (*  2a)  from  the  wave  peaks.  In  the  wave-following  frame 
(Case  III),  these  kinematic  effects  of  the  waves  are  totally  eliminated 
because  the  velocity  probe  is  essentially  at  some  constant  distance  from 
the  mobile  interface.  Furthermore,  the  velocity  profiles  at  low  wind 
speeds  clearly  show  the  existence  of  the  buffer  layer  and  the  viscous 
sublayer  which  cannot  be  revealed  in  the  fixed-frame  system.  The  exist¬ 
ence  of  these  layers  are  known  from  other  experiments,  such  as  those  of 


McLeish  and  Putland  (1975),  and  Wu  (1975).  However,  their  measurements 


were  Lagrangian  and  are  of  limited  value  in  comparing  with  Eulerlan 
data.  Here,  the  wave-following  frame  measurements  are  still  Eulerian 
quantities . 

In  the  wave-following  coordinates,  the  friction  velocities  for  the 
low  wind-speed  experiments  can  be  found  from  the  velocity  profiles  in 
the  viscous  sublayer  rather  than  relying  on  the  usual  velocity  profile 
curve-fitting  procedure.  The  friction  velocities  obtained  using  this 
technique  agree  well  with  the  direct  turbulent  shear  stress  measure¬ 
ments  . 

The  velocity  profile  curve-fitting  method  may  not  work  for  the  pro¬ 
files  in  this  study  (fixed  or  wave-following  frame)  because  the  waves 
affect  the  slopes  of  the  profiles;  thus,  the  slopes  need  not  be  constant 
in  the  logarithmic  region.  On  the  other  hand,  it  is  possible  to  obtain 
the  friction  velocity  without  directly  measuring  the  turbulent  shear 
stress  for  low  wind-speed  flows.  One  just  uses  the  wave-following 

frame  to  obtain  a  mean  velocity  down  or  near  to  the  viscous  sublayer 

—  « —  *1/2 

and  the  friction  velocity  is  given  by:  u*  ■  [v(us  -  u)/(-y  )]  '  as  in 
Sec.  4.3.4,  for  -y*  in  the  viscous  sublayer. 

In  Case  III,  the  mean  velocity  defect  profiles  in  wall  coordinates 
follow  the  logarithmic  trend  which  implies  that  the  mean  profiles  follow 
the  undulating  motion  of  the  mechanical  waves.  Thus,  the  dominant  waves 
do  not  contribute  to  the  "roughness”  of  the  flow.  The  ripples,  on  the 
other  hand,  seem  to  have  some  effect  on  the  flow  "roughness,"  as  the 
mean  profiles  shift  downward  slightly  with  increasing  wind  speed. 

The  turbulent  intensities  (RMS  u’  and  v’)  in  the  wave-following 
frame  appear  as  an  extension  of  the  fixed-frame  results  to  shallower 
depth  (-y*  <  a).  The  RMS  u’  and  v’  are  isotropic  and  increase  in 


magnitude  as  the  Interface  is  approached'  Similar  trends  are  observed 
in  the  fixed-frame  experiments,  but  the  behavior  of  the  turbulent  Inten¬ 
sities  near  the  interface  could  not  have  been  anticipated  from  the 
fixed-frame  measurements.  The  isotropy  of  the  turbulent  intensities  is 
unexpected,  but  may  be  explained  as  follows.  The  distribution  of 
ripples  appears  fairly  uniform,  so  the  random  ripple-induced  velocities 
are  expected  to  be  isotropic.  Furthermore,  the  ripple-induced  veloci¬ 
ties  are  strongest  near  the  interface.  In  both  Cases  II  and  III,  the 
ripple-induced  velocities  are  included  in  the  calculation  of  RMS  u'  and 
v*.  Thus,  the  turbulent  intensities  near  the  surface  are  mostly  domi¬ 
nated  by  the  ripple-induced  velocities.  Therefore,  the  RMS  u'  and  v' 
are  isotropic  due  to  isotropy  of  the  ripple-induced  velocities.  How¬ 
ever,  the  RMS  u'  and  v'  have  been  shown  to  scale  with  the  turbulent 
velocity  and  length  scales,  which  leads  us  to  believe  that  there  is  no 
ambiguity  in  defining  the  ripple-induced  velocities  as  turbulence.  Bole 
and  Hsu  (1969)  suggested  that  there  are  interactions  between  the  ripples 
and  the  dominant  waves,  and  since  the  ripple-induced  velocities  have 
turbulent  characteristics,  there  may  be  interactions  of  the  main  waves 
with  turbulence  through  the  action  of  the  ripples. 

The  linear  water-wave  theory  was  shown  to  give  fairly  good  results 
for  the  wave-induced  velocities  for  -y  >  a.  In  Sec.  5.5.1,  we  showed 
that  the  RMS  wave-induced  velocities  are  practically  coordinate  indepen¬ 
dent.  Here,  since  y  ■  h  and  y*  ■  h  +  fn  both  have  the  same  mean 
depth,  we  can  compare  the  calculated  wave-induced  velocities  given  by 
the  linear  theory  at  y  *  h  with  the  measured  values  at  y*  *  h  +  fn- 
The  linear  theory  also  gives  good  results  for  the  wave-induced  veloci¬ 
ties  for  -y*  >  a.  It  is  possible  to  calculate  the  wave-induced 


velocities  for  -y  <  a  though  these  velocities  may  not  he  defined  In 
the  fixed  frame*  Again,  we  can  compare  the  calculated  wave-induced 
velocities  for  -y  <  a  with  the  corresponding  measured  values  at 
-y*  <  a  because  y  m  y  *  y*.  The  calculated  velocities  agree  with  the 
measurements  fairly  well  at  u,,,  ■  0.0  and  1.7  m/s,  the  difference  Is 
about  102  for  the  RMS  u  and  only  about  62  for  the  RMS  v.  At  higher  wind 
speeds,  the  agreement  is  not  as  good  in  this  region  (-y*  <  a),  because 
the  calculated  values  include  the  ripple-induced  velocities  which  the 
measured  RMS  u  and  v  do  not.  This  exercise  shows  that  even  though  the 
calculated  velocities  are  from  a  linear  theory  in  the  fixed  frame,  the 
calculated  values  agree  very  well  with  the  measured  values  in  the  wave¬ 
following  frame  at  equivalent  mean  depths,  even  at  depths  less  than  the 
wave  amplitude. 

The  mean  wave-induced  shear  stress  -u  v  appears  to  have  trends  sim¬ 
ilar  to  those  in  the  fixed-frame  experiments,  except  for  the  case  of 
u„  =  2.5  m/s,  where  -u  v  is  positive  in  some  regions  of  the  flow.  This 
could  be  due  to  data  scatter,  but  the  distribution  of  -u  v  with  depth  is 
relatively  smooth  which  suggests  that  the  observed  trend  is  not  caused 
by  data  scatter.  Also,  at  about  this  wind  speed,  similar  changes 
in  -u  v  occur  on  the  air  side  of  the  interface  as  shown  by  Hsu  and  Hsu 
(1983).  However,  a  definitive  conclusion  cannot  be  made  at  this  stage, 
suggesting  further  experiments  around  this  wind  speed. 

The  turbulent  shear  stress  -u '  v  ’  measurements  are  of  the  correct 
order  of  magnitude  but  there  are  fair  amounts  of  data  scatter.  Consis¬ 
tent  friction  velocities  cannot  be  obtained  from  these  measurements. 
This,  however,  does  not  mean  that  the  wave-following  frame  is  unsuitable 
for  such  measurements,  but  just  that  the  limits  of  the  instruments  have 
been  reached. 
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Although  we  have  not  developed  the  energy  equations  in  the  wave- 
following  coordinates,  some  preliminary  analyses  show  that  many  of  the 
terms  are  very  similar  to  those  in  the  fixed-frame  equations.  There¬ 
fore,  we  made  calculations  on  some  of  the  dominant  kinetic  energy 
production  terms,  such  as  -u  v(3u/3y*)  and  -uTvT( 3u/3y*)  .  Their 
values  can  be  found  in  Appendix  D.  Some  representative  profiles  of 
these  terms  are  shown  in  Figs.  65  and  66.  Despite  the  data  scatter, 
the  term  -u  v(3u/3y*)  is  generally  negative,  but  it  is  always  positive 
very  close  to  the  interface  (the  first  or  second  data  points).  It  is 
unclear  whether  this  is  mere  coincidence  or  due  to  some  physical  phenom¬ 
ena;  only  further  experimentation  may  resolve  this  mystery.  If  the 
term  -u  v(3u/3y*)  is  indeed  positive  very  close  to  the  interface,  then 
this  implies  the  waves  draw  their  energy  from  the  mean  flow  in  that 
region.  This  result  is  not  unlike  what  Hsu  and  Hsu  (1983)  found  on  the 
air  side  of  the  interface. 

5.5.3  Summary 

This  first  attempt  to  extend  the  wave-following  coordinate  system 
from  the  air  side  of  the  interface  to  the  water  side  has  proven  to  be 
very  successful.  Not  only  is  the  wave-following  frame  able  to  extend 
the  fixed-frame  results,  but  also  interesting  features  which  could  not 
have  been  discovered  in  a  fixed  frame  are  revealed.  However,  refinement 
of  the  LDA  is  needed  in  order  to  resolve  higher-order  wave  or  turbulent 
terms . 
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CHAPTER  6 


CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  CONCLUSIONS 

The  main  objectives  of  this  study  have  been  accomplished.  The  fol¬ 
lowing  conclusions  can  be  drawn  from  the  results  and  discussions  pre¬ 
sented  in  Chaps.  4  and  5. 

1.  This  study  presents  a  comprehensive  data  set  on  momentum  transfer 
at  an  air-water  interface  under  a  wide  range  of  surface  conditions 
and  in  two  coordinate  systems.  The  data  set  should  prove  very  use¬ 
ful  in  evaluating  current  or  future  computational  schemes  as  well 
as  turbulence  models  for  air-sea  interaction. 

2.  A  nonintrusive  wave- following  velocity  probe  extends  the  wave¬ 

following  coordinate  system  (Hsu  et  al.,  1981)  from  the  air  side  of 
the  interface  to  the  water  side.  The  wave-following  velocity  probe 
can  reveal  features  in  the  flow  which  could  not  have  been  found 
from  a  fixed  frame. 

3.  It  is  possible  to  separate  the  flow  field  beneath  the  air-water 

Interface  into  three  constituent  fields:  the  mean,  the  wave,  and 
the  turbulent  fields.  Each  of  these  fields  can  be  characterized  by 
significant  velocity  and  length  scales. 

4.  A  simple  analysis  of  the  governing  equations  using  the  experiment¬ 
ally  obtained  characteristic  scales  show  that  there  are  interac¬ 
tions  among  the  mean,  wave,  and  turbulent  fields. 

5.  The  wave  field  affects  the  mean  flow  directly,  which  leads  to  a 

change  in  slope  in  the  steady  form  of  the  law  of  the  wall,  although 


the  logarithmic  region  persists  In  the  velocity  profiles.  Also, 
the  kinetic  energy  of  the  mean  flow  is  augmented  by  the  waves. 

6.  Small-amplitude  water-wave  theory  gives  good  predictions  for  the 

measured  wave  quantities. 

7.  Turbulence  draws  its  energy  mainly  from  the  mean  flow.  Turbulence 

can  be  augmented  by  the  waves  because  the  waves  affect  the  mean 

flow. 

8.  Turbulent  intensity  profiles  in  wind-wave  experiments  are  similar 
to  flows  over  flat  plates,  for  example  Klebanoff  (1955). 

9.  In  mechanical-wave  experiments,  the  turbulent  intensity  profiles 

have  very  different  trends  from  those  in  wind-wave  experiments. 
This  may  be  due  to  the  interaction  of  the  ripple-induced  velocities 
with  that  of  the  dominant  wave. 

10.  Direct  interactions  between  the  waves  and  turbulence,  such  as 

-r^ jOu^/BXj)  in  the  energy  equations,  cannot  be  quantified  satis¬ 
factorily.  These  interaction  terms  are  probably  small  compared  to 
other  types  of  interactions. 

6.2  RECOMMENDATIONS  FOR  FUTURE  WORK 

The  results  documented  here  represent  the  major  part  of  the  study. 

Further  analysis  of  the  data  may  include  the  following: 

1.  A  careful  study  of  the  velocity  spectra  may  reveal  how  the  various 
spectral  components  of  the  wave  and  turbulence  vary  with  depth  and 
wind  speed.  Then,  from  this  information,  we  can  examine  the  frozen 
turbulence  model  postulated  by  Lumley  and  Terray  (1982). 

2.  Velocity  profiles  at  fixed  phase  angles  along  the  waves  can  be 
obtained  for  mechanical-wave  experiments.  This  can  give  valuable 


information  as  to  how  the  wave  and  mean  velocities  change  with 
phase  (i.e.,  position)  along  the  main  wave. 

3.  The  mean  and  wave  momentum  equations  in  the  wave-following  coordi¬ 
nates  need  to  be  derived.  An  order-of-magnitude  analysis  can  then 
be  performed  on  these  equations  which  may  provide  some  further 
insights  into  the  wave-following  frame  results. 

4.  An  analysis  of  the  energy  equations  in  the  form  as  those  of  Hsu 
et  al.  (1982)  may  further  quantify  the  partition  of  energy  in  the 
vertical  and  horizontal  directions. 

The  merits  of  experimental  studies  about  the  air-water  interface 

are  clearly  illustrated  in  this  study.  The  following  are  some  recommen¬ 
dations  for  future  research: 

1.  The  wave-following  LDA  system  needs  to  be  improved  to  allow  more 
accurate  measurements  of  the  turbulent  shear  stresses. 

2.  The  present  experimental  program  can  be  run  at  different  fetches  to 
quantify  the  effects  of  wave  growth  on  water  velocity.  This 
enables  us  to  examine  the  gradual  development  of  the  interactions 
between  the  waves  and  the  mean  flow. 

3.  It  may  be  possible  to  develop  a  wave-following  coordinate  system 
for  the  wind-wave  experiments.  This  is  a  non-trivial  task  because 
the  random  nature  of  the  wind  waves  pose  the  biggest  challenge  in 
data  acquisition  and  reduction.  However,  a  great  deal  can  be 
learned  from  a  successful  wave-following  attempt  in  this  case. 

4.  A  complete  three-dimensional  map  of  the  velocity  field  in  the  water 
can  improve  on  the  present  data  set.  The  total  kinetic  energy  for 
the  mean,  wave,  and  turbulent  flow  fields  can  then  be  determined. 
This,  when  coupled  with  simultaneous  three-dimensional  velocity 


measurements  on  the  air  side,  may  very  well  be  the  ultimate  experi¬ 
ment  for  momentum  and  energy  transfers  at  an  air-water  interface. 

5.  A  study  of  the  total  energy  transfer  can  be  achieved  by  mapping  out 
the  temperature  field  in  some  of  the  above-mentioned  experiments. 
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APPENDIX  A 


WAVE- SEPARATION  METHOD 

The  theoretical  framework  for  the  wave-separation  method  is  out¬ 
lined  in  Sec*  2.2.2  and  Howe  et  al.  (1981).  Some  of  the  essential  steps 
are  given  below. 

A  A 

Let  x  and  y  be  some  measured  quantities  (such  as  u,  v)  at  a  fixed 
fetch.  Each  quantity  can  be  decomposed  into  three  parts,  i.e., 

a  ■  ~ 

X  ■  X  +  X  +  x' 

y  ■  y  +  y  +  y' 

Then,  after  defining 

x  ■  x  +  x' 

y  -  y  +  y\ 

a  linear  operator  L  is  defined  as  follows: 

x(t)  -  Jf  Lx(a)n(t  -  a)da  -  Lx(t)  *  n(t)  (A.l) 

and 

y(t)  -  _/”Ly(B)n<t  '  B)dB  "  Ly(t)  *  n(t>  (a. 2) 

where  n(t)  denotes  the  water  surface  elevation  at  the  same  fetch,  and  * 
denotes  the  convolution  process* 

The  cross-correlation  of  x  and  y  is  defined  as: 

x(t)y*(t  ~  t)  -  /  /"l  (a)L*(B)n(t  -  a)n*(t  -  t  -  gjdadB 

—oo  —oo  x  y 

where  *  denotes  the  complex  conjugate.  If  it  is  assumed  that  x  and  y 
are  the  outcomes  of  stationary  random  processes,  then  the  cross¬ 
correlation  function  is: 

CD  00  ^ 

Ro«(t)  -  /  /  L  (a)L  (e)R  (t  +  6  -  a)dadS 

y  -oo  -oo  x  y  hh 
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Fourier  transforming  the  above  equation  gives: 


APPENDIX  B 


ANALYSIS  OF  THE  OPTICAL  WAVE-FOLLOWER  SYSTEM 

In  this  appendix,  we  first  outline  the  procedures  used  in  calculat¬ 
ing  the  intersecting  point  of  the  three  LDA  laser  beans  when  they  are 
scanned  through  multiple  interfaces  (air,  glass,  and  water).  Then,  we 
show  that  the  Doppler  shifts  on  the  laser  beans  due  to  the  scanning 
mirror  are  identically  zero  when  the  wave-follower  is  used  in  a  homo¬ 
geneous  medium.  Next,  we  extend  the  results  for  the  single  medium  to 
the  air-glass-water  system  by  providing  some  heuristic  arguments.  A 
complete  theoretical  treatment  of  the  multiple  media  problem  is  too 
Involved  to  present  here.  Furthermore,  we  have  demonstrated  in  the 
wave-follower  qualification  tests  that  the  heuristic  arguments  are 
Indeed  correct. 

The  LDA  has  three  laser  beams  comprising  the  two-component,  two- 
color  system.  The  orientation  of  the  beam  is  shown  in  Fig.  B.l.  The 
blue  and  the  green  beams  lie  in  the  horizontal  plane  GOB  while  the  cyan 
beam  is  equidistant  from  both  the  blue  and  green  beams.  Also,  the  three 
beams  lie  on  the  surface  of  a  cone  with  its  vertex  at  the  focal  point  0 
and  the  base  at  the  front  lens  of  the  LDA  in  Fig.  B.l.  The  mirrors  in 
the  wave-follow  system  simply  fold  this  "cone"  at  different  sections. 
Therefore,  there  is  no  loss  of  generality  in  considering  the  section  of 
the  "cone”  which  begins  at  the  moving  mirror  instead  of  at  the  front 
lens  of  the  LDA.  Thus,  the  orientation  of  the  three  LDA  laser  beams 
when  they  traverse  through  the  air-glass-water  media  can  be  represented 
by  Figs.  B.2  and  B.3.  The  definitions  of  the  angles  and  distances  in 
the  figures  are  self-explanatory.  In  Fig.  B.2,  the  line  BOG  lies  along 


the  axis  of  rotation  of  the  scanning  mirror  and  BO'  and  GO'  represent 
the  paths  of  the  blue  and  green  laser  beams,  respectively.  Then  in  Fig. 
B.3,  CO'  represents  the  path  of  the  cyan  beam.  Note  that  00'  in  Fig. 
B. 3  represents  the  projection  of  the  surface  BO'G  in  Fig.  B.2  in  the 
plane  of  the  paper.  The  reader  is  cautioned  that  00'  does  not  represent 
any  laser  beam  in  that  plane  in  Fig.  B.3,  except  perhaps  at  the  point 
where  the  blue  and  the  green  beams  intersect.  Therefore,  the  angles 
along  00'  do  not  satisfy  Snell's  law  of  refraction.  In  is  clear  from 
Fig.  B.2  that  the  blue  and  green  beams  always  intersect.  However,  the 
intersection  point  of  the  blue  and  green  beams  may  not  intersect  exactly 
with  the  cyan  beam  for  non-zero  values  of  0g  because  the  sine  function 
in  Snell's  law  is  a  nonlinear  function. 

We  first  find  the  intersection  point  (xAB,  yAB)  of  the  blue  and 
green  beams.  This  point  is  taken  to  be  the  position  of  the  probe  vol¬ 
ume.  Using  Snell's  law  of  refraction,  some  geometry,  and  the  following 
nomenclature,  one  can  derive  the  key  relationships. 

The  nomenclature  is 

na  -  refractive  index  of  air  (1.00) 
n^  -  refractive  index  of  glass  (1.50) 
nw  -  refractive  index  of  water  (1.33) 

m  -  thickness  of  glass  (6.4  mm) 

£  -  distance  of  axis  of  rotation  of  scanner  mirror  to 

glass  (72  mm) 

f  -  distance  of  axis  of  rotation  of  scanner  mirror  to 
measuring  point  in  air  (411.5  mm) 

a  -  half  distance  between  the  blue  and  green  laser  beams 
at  the  optical  scanner  mirror  (25.6  mm) 


a 


half  intersecting  angle  of  LDA  beams  in  air  (2.5°) 


In  Fig.  B.2,  we  have  Chen 


na  sin  9S  -  ng  sin  0a  -  sin  0a 


a*  -  (l  -  , — L_)  a 
v  £  cos  6J 


a"  -  a'  -  a  sin  o'  tan  0' 

Q 


«•  -  tan-1 

vsin  QJ 


kAB 


£  +  m  +  a"( 


cot  0" 


8 


sin  a 


t) 


Yab  =  £■  tan  0  +  m  tan  0a  cos  o'  +  a"  cot  o'  *  f  sin  9 


In  Fig.  B.3,  we  have 
cos 

COS  ijl" 


cos  9' 

_ a _ 

sin  a'  sin  0' 
s 

c°»  e; _ 

sin  o'  sin  0" 
s 


Note  that  i|i'  *  0a  ,  41"  *  0"  in  general,  but  <ji  *  0. 

..  _  _ CM  ~  CP _ 

XPM  “  [tan  <JT  -  tan(0"  +  A0")] 


ypM  “  *pM(tan  ♦")  +  cp 
Cp  -  yp  -  xM(tan  f ) 

CM  "  "  xM^tan(0"  +  A9")l 

yp  *  yp  +  (tan  *')(xM  -  x^) 

yM  “  yM  +  l tan c e ’  +  A9’)](xM  -  Xjl) 
na  sin  6  -  ng  sin(9’  +  A0')  -  riy  sin(0"  +  A0") 

yM  “  yp  *  h  *  yp  ~  stn(9°+  °T  *  k  *  tan  a(f  “  1  sec  9) 


yp 


£  tan  9 


The  equation  for  the  line  c^O'  is  given  by 

~  ( tan( 0"  +  A0" ) ]  xM  -  c„  -  0  (B-3) 

Using  Eqs.  B-l  and  B-2,  one  can  calculate  the  location  of  the  point 

^xAB,yAB^*  T**e  distance  between  (xAB,yAB^  an<*  the  line  ^M0' 

(Eq.  B-3)  gives  the  limit  on  the  maximum  scan  angle  (0g)  for  a  given 
laser  beam  width  (at  the  probe  volume)  because  if  this  distance  exceeds 
the  beam  diameter  near  the  probe  volume,  the  cyan  beam  does  not  inter¬ 
sect  with  the  blue  and  green  beams  and  the  probe  volume  is  undefined. 
In  this  study,  the  beam  waist  was  about  0.15  mm  and  0g  was  limited  to 
about  7°  in  scanning  through  multiple  interfaces.  At  some  points  along 
the  scan  path,  the  three  laser  beams  may  not  intersect  exactly  at  the 
beam  waist,  leading  to  perhaps  additional  broadening  ot  the  LDA  signal 
as  pointed  out  by  Hanson  (1976).  However,  such  broadening  cannot  be  too 
severe  because  of  the  relatively  long  focal  length  of  the  front  lens. 
This  was  further  confirmed  by  the  results  of  the  qualifying  experiments. 

Figure  B.4  depicts  the  side  view  of  the  optical  wave-follower  when 
it  is  used  in  a  homogeneous  medium.  The  definitions  of  the  various 
angles  and  lengths  are  apparent  in  the  figure.  The  point  O',  which 
represents  the  LDA  probe  volume,  moves  with  some  velocity  Up  as  the 
scanning  mirror  rotates  with  an  angular  velocity  0.  In  addition,  we 
have  the  following  relationships: 

6  *  0max  8in  “ot 
:  d0 

0“-Tr“ui0  cos  iii  t 
dt  o  max  o 


where  fQ  is  the  frequency  of  oscillation  of  the  scanner  mirror  and  9MX 
is  the  maximum  scan  angle. 

In  Fig.  B.4,  CAO'  represents  the  path  of  the  cyan  laser  beam,  GOO* 
represents  the  projection  of  the  surface  where  the  blue  and  green  beams 
lie  onto  the  plane  of  the  figure,  and  MM*  represents  the  scanner  mirror. 

Because  the  direction  of  Up  is  perpendicular  to  00',  up  has  no 
velocity  component  in  the  plane  containing  the  blue  and  green  beams 
between  00'.  Furthermore,  the  blue  and  green  beams  reflect  off  the 
mirror  along  the  rotational  axis  where  the  velocity  of  the  mirror  is 
zero.  Thus,  the  blue  and  green  beam  frequencies  undergo  no  Doppler 
shift  due  to  the  scanning  mirror. 

When  the  mirror  ro‘  tes  through  an  angle  0,  the  point  O'  (see  Fig. 
B.4)  moves  through  an  arc  length  of  f ( 2 0 ) .  There  is  a  factor  of  2 
because  from  the  law  of  reflection,  a  reflected  ray  rotates  twice  the 
angle  0  as  the  mirror  rotates  through  an  angle  0.  Thus,  the  tangential 
velocity  of  the  point  O'  is  the  time  rate  of  change  of  the  arc  length, 
i.e.,  up  ■  f  ( 20) .  The  velocity  up  has  a  component  in  the  direction  of 
O' A,  so  we  have 

Up  cos(^  +  a]  *  -  Up  sin  a  ■  -  2f6  sin  a 

which  causes  a  frequency  shift  of  (-2f0  sin  a)/X  to  the  cyan  beam,  where 
X  is  the  wavelength  of  the  cyan  beam.  The  point  A  on  the  mirror  has  a 
tangential  velocity  u  ■  r9  which  has  components  in  the  directions  AC 
and  AO',  and  the  net  frequency  change  to  the  cyan  beam  due  to  u  is 
(2u/X)  cos[n/2  -  (8  -  B  +  a) ) .  Therefore,  the  net  change  of  frequency 
of  the  cyan  beam  as  it  travels  from  C  to  O'  is 


8in( 9  -  8  +  a) 


(B-4) 


because 

f  sin  a  *  r  sin(0  -  8  +  a)  =  d  (B-5) 

where  d  is  indicated  by  the  dashed  line  in  Fig.  B.4.  Equations  B-4  and 
B-5  are  very  general  equations,  the  only  constraint  is  that  the  light 
beams  originate  within  the  cone  of  light  bounded  by  the  LDA  lens  and  the 
focal  point.  Hence,  the  scanning  mirror  produces  no  Doppler  shift  to 
the  outgoing  laser  beams  and  the  incoming  backscattered  light  when  the 
wave-follower  is  used  in  a  homogeneous  medium. 

When  the  wave-follower  is  used  in  the  air-glass-water  system,  Eqs. 
B-4  and  B-5  are  still  valid  for  the  following  reason.  The  point  O'  in 
Fig.  B.3  indicates  the  actual  position  of  the  probe  volume  inside  the 
test  facility.  However,  to  an  observer  outside  the  facility,  the  appar¬ 
ent  intersection  point  of  the  laser  beams,  which  is  defined  by  the  beam 
angles  in  the  air,  is  at  0".  Thus,  as  far  as  the  receiving  optics  of 
the  LDA  is  concerned,  the  probe  volume  still  appears  as  if  it  were  in  a 
homogeneous  medium  and  the  scanning  mirror  produces  no  Doppler  shift  to 
the  laser  beams.  The  fact  that  the  beams  traverse  through  a  series  of 
media  and  back  only  constrains  the  maximum  scan  angle  but  does  not 
introduce  any  further  frequency  shift  to  the  beams. 

In  practice,  it  is  almost  impossible  to  have  the  aluminized  surface 
of  the  scanner  mirror  passed  through  the  axis  of  rotation  of  the  scan¬ 
ner.  The  effect  of  this  is  that  the  probe  volume  does  not  follow 
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exactly  an  arc  as  the  mirror  rotates.  This  minor  defect  can  be  mini¬ 
mized  by  aligning  the  blue  and  the  green  laser  beams  along  a  line  on  the 
scanner  mirror  closest  to  the  axis  of  rotation  and  parallel  to  it. 
Then,  the  maximum  deviation  of  the  probe  volume  path  from  an  arc  is 
0.1  mm  for  a  scan  angle  of  10 * . 


APPENDIX  C 


LEAST- SQUARES  COSINE-FITTING  TECHNIQUES 

To  fit  a  sequence  of  data  y^ti),  i  ■  1,  2,  n  with  a 

waveform  described  by: 

y  ■  A  cos(u)t  -  6)  , 

one  has  to  find  A  and  9  such  that 

S  ■  E(yi  -  A  cos(tot^  -  9))^ 

is  a  minimum.  The  conditions  for  the  existence  of  a  minimum  are: 


Thus , 


j)S  _  3S 
3A  *  39 


Z  cos(o)t^  -  8)  *  A  Z  cosz(wtj,  -  9) 


Z  y^  8in(iut^  -  8)  ■  A  I  sinCujt^  -  9)  cos(u)ti  -  9) 


Let 


Si  •  Z  yi  cosujti 


S2  ■  X  yi  sinmti 


S3  «  Z  cos^mti 


Z  sin  oitj 


55  =  Z  cosuiti  sinuiti 


then  Eqs.  C.l  and  C.2  can  be  written  as  follows: 


Si  cos9  +  S2  8in9  *  A(S3  cos^9  +  S^  sin^9  +  2Sj  cos9  sin9) 


S2  cos9  -  Si  sin9  «  A(Sj  cos^Q  -  S5  sin^9  +  (S4  -  SjJcosS  sin0) 


cosine 


(C.l) 

(C.2) 


(C.3) 

(C.4) 


Solving  for  9  and  A  from  Eqs.  C.3  and  C.4  gives: 
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Figure  63.  (a)  -u'v'  3u/3y  profiles  in  physical  units  (WW-EF) 
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Figure  65.  (b)  -uv  3u/3y*  profiles  in  dimensionless  units  (MW-WF-EF) 
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Figure  B.4.  Optical  wave-follower  in  a  homogeneous  medium. 
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ABSTRACT 


y  Because  the  global  weather  and  ocean  circulations  are  driven  by  transfers 
across  the  atmosphere-ocean  interface,  studies  of  It  remain  at  the  forefront 
of  research  refforts>in  atmospheric  sciences  and  oceanography.^  Many  questions 
concerning  the  momentum  and  energy  transfers  In  the  water  L  surface  layer 
beneath  the  air-water  Interface  remain  unanswered.  The  present ;l^kperlmental 
program  was  undertaken  with  the  goal  of  learning  how  the  energy,  which  Is 
transferred  from  the  wind,  is  distributed  among  the  mean,  wave,  and  turbulent 
flow  fields  in  the  water. ^  The  interactions  between  these  three  flow  fields 
were  also  examined. 

The  experiments  were  performed  in  the  Stanford  Wind,  Water  Wave  Research 
Facility.  A  two-component,  frequency-shifted,  laser  Doppler  anemometer  system 
was  used  with  frequency  trackers  to  measure  the  water  velocity.  Three  major 
sets  of  experiments  were  done:  i)  wind-generated  waves  in  Eulerlan  coordi¬ 
nates  with  wind  speeds  ranging  from  1.5  m/s  to  13.1  m/s j  ii)  wind-ruffled 
mechanically  generated  waves  (22  mm  in  amplitude  and  at  a  frequency  of  1  Hz) 
in  Eulerian  coordinates  with  wind  speeds  ranging  from  0.0  to  6.2  m/s,  and 
ill)  wind-ruffled  mechanically  generated  waves  in  wave-following  Eulerian 
coordinates  with  wind  speeds  ranging  from  0.0  to  6.1  m/s.  A  novel  optical 
wave-follower,  which  allowed  velocity  measurements  very  close  to  the  mobile 
surface,  was  constructed  and  deployed  in  the  wave-following  measurements.  The 
mean,  wave,  and  turbulent  fields  were  isolated,  and  from  them,  other  quanti¬ 
ties  of  interest  were  computed. 

1,i  ■  / 

The  results  show  that  the  waves  affect  the  mean  flow,  but  the  mean  velo¬ 
city  profiles  show  logarithmic  behavior.  The  wave  field  generally  agrees  with 
the  prediction  of  a  linear  water-wave  theory.  In  the  wind-wave  experiments, 
the  turbulent  quantities  behave  similar  to  those  in  flows  over  flat  plates. 
In  the  mechanical-wave  experiments,  the  turbulent  quantities  have  different 
trends  from  those  of  the  wind-wave  cases,  which  suggests  possible  weak  wave- 
related  turbulence.  The  mean  wave-induced  shear  stress  is  negative,  which 
implies  that  the  waves  augment  the  mean  flow.  Thus,  the  waves  can  also  trans¬ 
fer  energy  indirectly  to  turbulence  via  the  mean  flow.  The  turbulent  field 
draws  its  energy  from  the  mean  flow  through  the  mean  turbulent  shear  stress 
and  from  the  wave  field  via  the  wave-induced  turbulent  stresses.  ^Experimental 
observations  agree  with  simple  analysis  of  the  momentum  and  energy  equations. 
This  study  contributes  to  the  basic  understanding  of  the  fundamental  princi¬ 
ples  which  govern  transfers  at  an  air-water  interface. 


7 


r'. 

/, 


END 

FILMED 

7-85 

DTIC 


